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1.        INTRODUCTION 

In  1966  Calspan  Corporation  (formerly  Cornell  Aeronautical 
Laboratory,  Inc.)  began  development  of  a  general  mathematical  model  and 
computer  simulation  of  the  dynamic  responses  of  an  automobile  in  accident 
situations  under  Contract  CPR-11-3988  with  the  Bureau  of  Public  Roads. 

The  mathematical  model  of  vehicle  dynamics  developed  in  the  first 
year  of  that  effort  included  the  general  three-dimensional  motion  resulting 
from  vehicle  control  inputs,  traversal  of  irregular  terrain,  or  from  collisions 
with  simple  roadside  barriers.   The  model  was  subsequently  named  the 
Highway- Vehicle-Object  Simulation  Model  (HVOSM) <   Later,  the  model  was 
further  developed  and  a  comprehensive  validation  program  was  carried 
out  including  a  series  of  repeatable  full-scale  tests  with  an  instrumented 
vehicle  in  order  to  objectively  assess  the  degree  of  validity  of  the  vehicle 
model.   Extensive  measurements  of  the  vehicle  parameters  required  for  input 
to  the  HVOSM  were  made  under  a  subcontract  with  the  Ford  Motor  Company  as 
a  part  of  the  validation  procedure.  This  effort  was  reported  in  Reference  5 
and  the  model  as  described  therein  has  been  referred  to  as  the  V-3  version 
of  the  HVOSM. 

Modifications  were  subsequently  made  to  the  simulation  in  order  to 
study  the  effects  of  terrain  (specifically,  railroad  grade  crossings)  on 
vehicle  controllability.  The  impact  routines  were  removed  and  extended 
terrain  definition  capabilities  were  added  along  with  a  more  realistic 
model  of  suspension  properties.  This  program  version  (Reference  10)  has 
been  informally  referred  to  as  the  V-4  version  of  the  HVOSM  and  has  since 
been  used  extensively  for  study  of  roadway  and  roadside  geometries. 


Further  developments  of  HVOSM  aimed  at  providing  a  simulation 
model  more  suitable  for  the  study  of  the  complex  dynamics  resulting  from 
accident  avoidance  evasive  maneuvers  were  reported  in  Reference  11,   This 
version,  informally  called  the  V-7  version  of  the  HVOSM,  includes  a  detailed 
model  of  the  braking  and  engine-driveline  systems  and  an  empirically  based 
definition  of  the  relationships  between  longitudinal  and  lateral  tire  forces 
through  the  inclusion  of  rotational  degrees  of  freedom  of  the  four  vehicle 
wheels. 

During  development  of  the  HVOSM,  documentation  efforts  primarily 
fulfilled  the  objectives  of  maintaining  communication  within  the  program 
development  structure,  ensuring  quality  control  of  the  development  and 
providing  a  historical  reference.   It  was,  however,  recognized  early  in 
the  development  of  the  HVOSM,  that  this  state-of-the-art  advance  in  the 
modeling  of  a  vehicle  and  its  environment  could  be  put  to  best  use  through 
its  widespread  distribution  to  organizations  interested  in  its  application  to 
highway  safety.   As  a  result,  distribution  of  the  HVOSM  was  begun  before 
its  development  was  complete  and  before  instructional  documentation  could 
be  provided. 

Recognizing  the  need  to  bring  documentation  of  the  several  HVOSM 
versions  together  and  to  provide  the  highway  safety  community  with  an  effective 
description  of  the  programs  and  their  use,  the  Federal  Highway  Administration 
(FHWA)  awarded  Calspan  Corporation  contract  number  DOT-FH-11-8265  for  the 
purpose  of  providing  such  documentation  for  the  then  existing  versions  of 
the  HVOSM. 

Three  versions  of  the  HVOSM  were  covered  by  this  documentation. 
They  were,  the  HV0SM-SMI1  (Sprung  Mass  Impact)  version  (formerly  known  as 
the  V-3  version),  the  HV0SM-RD1  (Roadside  Design)  version  (formerly  known 
as  the  V-4- version)  and  the  HV0SM-VD1  (Vehicle  Dynamics)  version  (formerly 
known  as  the  V-7  version)..   Under  the  first  phase  of  that  effort,  only 
those  versions  as  developed  by  Calspan  were  covered  by  the  documentation. 


The  second  phase  of  contract  number  DOT-FH-1 1-8265  called  for 
extension  of  the  capabilities  of  the  HVOSM  by  adding  new  features,  including 
some  additional  modifications  made  by  other  research  organizations,  and 
providing  additional  ease  of  use  features. 

Accordingly,  Calspan  has: 

•  Generalized  the  basic  vehicle  model  to  include  the 
capability  for  simulating  an  independent  front  and 
rear  suspension  vehicle  and  a  vehicle  with  solid 
front  and  rear  axles. 

•  Generalized  the  tire  model  to  allow  specification 
of  up  to  four  different  tires  on  a  vehicle  and 
revised  the  friction  ellipse  tire  model. 

•  Combined  the  sprung  mass  impact  version  with  the 
roadside  design  version  resulting  in  only  two 
program  versions  at  the  end  of  the  second  phase. 

•  Incorporated  the  Preview-Predictor  Driver  Model 
described  in  Reference  15  into  the  vehicle  dynamics 
model. 

•  Incorporated  impact  forces  due  to  localized 
structural  hard  points  into  the  sprung  mass  impact 
algorithm.  This  modification  was  originally 
developed  by  the  Texas  Transportation  Institute 
(TTI)  and  was  added  as  reported  in  Reference  13. 


•  Extended  the  curb  impact  algorithm  to  allow  up  to 
six  planes  to  describe  a  curb.   This  modification  was 
also  developed  by  TTI  and  was  reported  in  Reference  16. 

•  Developed  a  road  roughness  algorithm  to  allow 
determination  of  the  effects  of  road  roughness  on 
vehicle  performance. 

•  Revised  input  and  output  format  to  provide  an  easy  to 
use,  more  flexible  data  interface, 

•  Developed  a  Pre-Processing  Program  to  calculate  a 
number  of  program  inputs  including  vehicle  and 
terrain  data  or  to  supply  input  cards  from  a  stored 
library  of  vehicle  data. 

The  documentation  provided  now  covers  the  two  program  versions: 
the  HV0SM-RD2  Version  (Roadside  Design)  and  the  HV0SM-VD2  Version  (Vehicle 
Dynamics) .   It  is  intended  to  be  a  base  to  which  further  developments  and 
modifications  to  the  HVOSM  can  be  added,  thus  providing  a  uniform  reporting 
format  and  centralized  source  of  information  for  the  many  HVOSM  users,   It 
consists  of  four  volumes,  each  describing  a  separate  aspect  of  the  HVOSM. 
Two  volumes  are  directed  toward  the  engineer/analyst  containing  the  analysis 
(derivation  of  governing  equations,  assumptions,  and  development  of  controlling 
logic)  and  experimental  validation.   Another  volume  is  directed  toward  the 
general  program  user  and  contains  analysis/program  symbology,  descriptions 
of  the  models  and  solution  procedures,  descriptions  of  input  requirements 
and  program  output,  and  a  number  of  program  application  examples.   The 
remaining  volume  of  documentation  is  intended  for  use  by  those  interested 
in  the  detailed  computer  programs.   This  volume  contains  descriptions  of 
the  computer  code  including  a  discussion  of  subroutine  functions,  annotated 
flowcharts  and  program  listings.   Also  included  are  a  list  of  program  changes, 
a  description  of  program  stops  and  messages,  and  computer  system  requirements 
necessary  to  run  the  programs. 
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The  HVOSM  analysis  is  covered  in  this  report.  Section  2  contains 
a  list  of  symbols  and  their  definitions  and  the  notation  used  in  the  analysis 
The  analytical  derivations  of  the  HVOSM  mathematical  model  are  presented 
in  Section  3,  and  references  are  listed  in  Section  4. 


2.        SYMBOLOGY  AND  NOTATION 

2. 1  Symbology 

The  HVOSM  symbology  is  presented  in  this  section.   The  listing 
of  symbols  is  ordered  with' respect  to  analytical  symbol  and  includes  a 
corresponding  program  symbol,  a  brief  definition  and  an  equation  number 
referencing  the  calculation  of  the  variable.   Input  variables  are  indicated 
by  an  I  in  the  equation  number  column. 

2.2  Notation 

The  time  derivative  of  a  variable  is  indicated  by  a  dot  over  the 

2 
symbol  for  the  variable  (i.e.,  a  =  da  ,  a  =  d  a  ). 

dt        at2 

The  following  subscript  notation  is  employed: 

0  =  initial  value  of  a  variable  at  zero  time. 

t  =  value  of  a  variable  at  the  end  of  the  current  time 
increment . 

t-1  =  value  of  a  variable  at  the  end  of  the  previous 
time  increment. 

F  =  front . 

R  =  rear,  or  rear  axle. 

i  =  wheel  identification  --  1,  2,  3,  4  =  RF, 
LF,  RR,  LR,  respectively. 

j   =  identification  of  vehicle  end.   j  =  F,  R 
for  front  end  and  rear  end,  respectively. 

s  =  sprung  mass. 

u  =  unsprung  mass. 


ANALYTICAL 
SYMBOL 


Wci 


APD 

APDHAX 


as:A.,c5. 


d    ,b  .c 


w% 


[A] 


[Aj] 


AMU 

AMUG 
(AP)F 

(AP)0 


PROGRAM 
SYMBOL 


APD 
APDMAX 


AS(4) 
BS(4) 
CS(4) 


AX(4) 
BX(4) 
CX(4) 

AY  (4) 
BY(4) 
CY(4) 

AHTX(3,3) 


AINTI 
AJMTX(3,3) 

AMU 

AMU6(5) 
APF(21) 

APR(21) 


155 


39 


iot 


53 


m 


Distance  along  vehicle  fixed  x 
axis  from  the  sprung  mass  center 
of  gravity  to  the  center  line  of 
the  front  wheels 

Directional  components  of  a  Hne 
perpendicular  to  both  the  normal 
to  the  wheel  plane  and  the  radial 
tire  force,  F„ 
K1 

Accelerator  pedal  deflection 
and  maximum  accelerator 
pedal  deflection 

Directional   components  of  a  line 
perpendicular  to  both  a  normal  to 
the  tire-terrain  contact  plane 
and  the  line  of  intersection  of 
the  wheel  and  ground  planes 

Direction  components  of  a  line 
perpendicular  to  both  a  normal 
to  the  tire-terrain  contact  plane 
and  the  vehicle  fixed  y  axis 

Direction  component  of  a  Hne  per- 
pendicular to  both  a  normal  to  the 
tire-terrain  contact  plane  and  the 
vehicle  fixed  x  axis 

Transformation  matrix  from  vehicle 
fixed  to  space  fixed  coordinate 
systems 

Intersection  area  of  cutting  plane 
1  with  the  sprung  mass 

Transformation  matrix  from  wheel 
fixed  to  vehicle  fixed  coordinate 
systems 

Tire-terrain  friction  coefficient 
at  zero  speed  and  nominal   tire 
loading 

Tire-terrain  friction  coefficient 
factor  for  5  terrain  tables 

Anti -pitch  coefficients  for  front 
suspension  positive  for  anti- 
pltch  for  forward  braking 

Anti-pitch  coefficients  for  rear 
suspension,  effective  at  the 
wheels;  positive  for  anti -pitch 
effect  for  forward  braking 


in. 


in 


in' 


lb/lb-ft 


Ib/lb-ft 
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A3.A4 
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,FP1 
*FP2 


r<y 


CO 

cF.cR 


c    c 


[C,] 


CONS 


^ 
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ARBRF 
ARBRR 


A0.A1.A2 


A3.A4 


BFP1 
BFP2 


BNMTX(3,3) 


'"So 
fCTJ 


CFP 
CRP 


CMTX(3.4) 

CONS 

CNMTX(3,3) 

RRMC(4) 
TCT(12) 


Drive  axle  ratio  (propeller  shaft 
speed/wheel  speed).  Default  of 
1.0 

[     Constant  coefficients  for  tire 
side  force  due  to  slip  angle 

[     Constant  coefficients  for  tire 
side  force  due  to  camber  angle 

Distance  along  the  vehicle  fixed 
x  axis  from  the  sprung  mass  cen- 
ter of  gravity  to  the  centerline 
of  the  rear  wheels   (entered  posi- 
tive) 

BMTX(3,3)  I'M  Transformation  matrix  from  wheel 
fixed  to  space  fixed  coordinate 
systems 

I     First  and  second  order 

coefficients  for  relationship 
between  brake  pedal   force 
and  brake  system  pressure 
Transformation  matrix  from 
orientation  of  vehicle  axes  at 
indexing  to  space  fixed  axes 
(Euler  angles- %', 6'n,Q^      ) 

215!  Small  angle  camber  stiffness 


no 


Front  and  rear  viscous  damping 
coefficient  for  a  single  wheel, 
effective  at  the  wheel  for  the 
front  and  at  the  spring  at  the 
rear 

Front  and  rear  coulomb  damping 
for  a  single  wheel,  effective 
at  the  wheel   for  the  front  and 
at  spring  for  the  rear 

Coefficient  matrix  for  simul- 
taneous solution  of  the  ground 
contact  point 

Ratio  of  conserved  energy  to 
total  energy  absorbed  by  the 
sprung  mass 

Transformation  matrix  from 

vehicle  fixed  axes  to  most 

recently  indexed  axes  (Euler 
angles**  Wt-et  >  0c     ) 

Rolling  resistance  moment 
coefficient 

Small  angle  cornering  stiffness 

Closed  throttle  engine  torque 


psi/lb, 
psi/lb^ 


lb/rad 
Ib-sec/in 


lb 


lb-in/lb 

lb/rad 
lb-ft 
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![D] 


Dax 

DELB 
DELE 
DDEL 

DIST 
DRWHJ 


WiAi 


EN 

FAPi 
FARi 


FBRK 

F 

Fcxui 
Fcyui 
Fczui 


rJFi 
(FnL 
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CT(4) 
CPSP 


CONE 
CTWO 
CTHREE 

DMTX(10,11 


DELTAX 

DELB 
DELE 

DDEL 

DIST 


DRWHJ 


Dl(4) 

D2(4) 
D3(4) 

EI 


EN 
APITCH 

FB 

FBRK 
FC(4) 

FCXU(4) 
FCYU(4) 
FCZU(4) 

FJP(35) 

FJF(4) 
FN 


40 


m 


Circumferential   tire  force 
stiffness 

Coulomb  resistance  torque  in 
the  steering  system  effective 
at  the  wheels 

Coefficients  in  relationship 
approximating  aerodynamic  and 
rolling  resistance 

Mass  matrix  of  coupled  second 
order  differential   equations. 
Column  11   contains  the  forcing 
functions 

Desired  vehicle  acceleration 

Beginning,  end,  and  incremental 
wheel   deflection  for  entered 
front  wheel   camber  table 

Desired  speed  differential 
nulling  distance 

Incremental   tire  deflection  for 
calculation  of  the  equivalent 
tire  force-deflection  character- 
istic in  the  radial  mode 

Direction  components  of  a  line 
perpendicular  to  the  normals  of 
both  the  wheel   plane  and  the  tire- 
terrain  contact  plane 

Error  between  predicted  and 
desired  path  at  the  ith 
viewing  position 

Number  of  points  at  which  ei 
is  determined 

Anti-pitch  force  at  wheel   i 

Force  at  wheel   i   due  to 
auxiliary  roll  stiffness 

Resistance  force  normal   to  the 
contact  surface  of  a  deformable 
barrier 

Brake  pedal   force 

Circumferential   tire  force 

Components  of  the  circumferential 
tire  force  along  the  x,y,  and  z 
axes     • 

Table  of  equivalent  radial  spring 
forces  as  a  function  of  deflection 

Jacking  force  at  wheel  i 

Vehicle  force  produced  by  deforma- 
tion of  the  vehicle  structure  nor- 
mal  to  the  contacted  surface 


lb 
lb-in 


in/sec' 


in 


lb 
lb 

lb 

lb 
lb 
lb 

lb 

lb 
lb 
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NST1 


FRi 


(FRICT) 


pRxui 

Ryui 

FRiui 


L  F 


Rx  i 
Ry  i 
Rz'i 


*  si  max 
IF  . 


Fsxui 
Fsyui 
Fszui 

Fxui 
Fyui 
Fzui 

Z  Fxu 

t  Fyu 

ZFys 
ZFzs 


Z  F 


Zl 


1F1 


IRi 
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FNSTI(3) 
FR(4) 

FRCP(4) 

FRICT 

FRXU(4) 
FRYU(4) 
FRZU(4) 

SFRX(4) 
SFRY(4) 
SFRZ(4) 

FS(4) 


SFXS 


FSXU(4) 
FSYU(4) 
FSZU(4) 

FXU(4) 
FYU(4) 
FZU(4) 

SFXU 
SFYU 

SFYS 

SFZS 

SFZ1 


F1FI(Z) 

F1RI(2) 


11D 


B> 


3H 


Structural  hard  point  force 

Radial  tire  force  in  the  plane  of 
the  wheel 

Tire  force  perpendicular  to  the 
tire-terrain  contact  plane 

Friction  force  acting  between  the 
vehicle  sprung  mass  and  barrier 


Components  of    F' . 

sprung  mass  axes     for  wheel   i 


r^    along  the 


Summation  of  the  components  of 
radial  spring  mode  forces  over 
tire  i,  with  respect  to  space 


Tire  side  force  in  the  plane  of  the 
tire-terrain  contact  patch  perpen 
dicular  to  the  line  of  intersection 
of  the  wheel   plane  and  ground  plane 

Resultant  side  force  correspond- 
ing to  small   angle  properties  for 
slip  and  camber  angles 

Maximum  achievable  side  force  as 
limited  by  the  available  friction 

Sprung  mass  Impact  force  or  com- 
bination of  rolling  resistance 
and  aerodynamic  drag  acting  along 
the  vehicle  x  axis 

Components  of  tire  side  force, 
F  .   along  the  sprung  mass  axes 

Total   tire  force  components  along 
the  vehicle  axes 


2*1 

3S3|  Resultant  forces  acting  on  the 
351)  vehicle  through  the  unsprung 
masses  in  the  x  and  y  directions 


Sprung  mass  impact  force  acting 
along  the  vehicle  y  axis 

Sprung  mass  impact  force  acting 
along  the  vehicle  z  axis 

Resultant  force  transmitted 
through  the  suspensions  in  the 
z  direction 

Front  and  rear  suspension  coulomb 
damping  forces  for  a  wheel,  effec- 
tive at  the  wheel   for  the  front 
and  at  the  spring  for  the  rear 


lb 

lb 
lb 
lb 

lb 
lb 

lb 

lb 

lb 
lb 

lb 

lb 

lb 

lb 
lb 
lb 

lb 
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DEFINITION 

UNITS 

r2Ri 

F2FI(2) 

m 

Front  and  rear  suspension  spring 

lb 

h 

XIR 

I 

Rear  unsprung  mass  moment  of 

2 

lb-sec  -in 

F2RI(2) 

\&> 

and  bumper  forces  for  a  wheel , 

inertia  about  a  line  through 

effective  at  the  wheel   for  the 

its  center  of  gravity  and 

front  and  at  the  spring  for  the 

parallel  to  the  vehicle  x  axis 

2 

lb- sec  -in 

rear 

L-i 

FIWJ(4) 

I 

Rotational   inertia  of  an  indi- 

9 

G 

i 

Acceleration  due  to  gravity 

2 

in/sec 

WJ 

vidual  wheel   at  the  front  or  rear 

2 
lb-sec  -in 

GEAR, 

GEAR1 

I 

Transmission  gear  ratios 

— 

Ix.Iy.Iz 

XIX 

I 

Spring  mass  moments  of  inertia 

gear' 

GEAR, 
GEAR^ 

GEAR2 

XIY 

about  the  vehicle  axes 

GEAR  3 

XIZ 

GEAR4 

Ixz 

XIXZ 

I 

Spring  mass  roll-yaw  product  of 

lb-sec  -in 

Gli 

GN(1,J) 

i 

Lever  arm  lengths  in  brake  types 

in 

inertia 

lj 

1,2  and  3 

(I'x)t 

+7 

Effective  inertial  term    due  to 

G,4 

GN(2,J) 

I 

Brake  actuation  constant,  assumed 

time  varying  positions  of  the 

2j 

to  be  equal   for  both  shoes  of 
brake  types  1   and  2 

(I'z)t 

47 

unspring  masses 

Effective  inertial   term    due  to 

G3j 

GN(3,J) 

i 

Effective  lining-to-drum  or  lining 

time  varying  positions  of  the 

to-disk  friction  coefficient  at 

unspring  masses 

design  temperature  for  all  shoes 

(I'xz)t 

47 

Effective  inertial   term    due  to 

or  disks  in  types  1,2  and  4  and 

time  varying  positions  of  the 

for  the  primary  shoe  of  type  3 

in2            • 

unspring  masses 

G4j 

GN(4.J) 

i 

Cylinder  area  for  actuation  of 

(I'yz)t 

m 

Effective  inertia     term    due  to 

leading  shoe  of  brake  type  1,  or 

time  varying  positions  of  the 

for  each  shoe  in  types  2  and  3. 

unspring  masses 

Also  used  for  total   cylinder  area 

ik        2   . 

per  side  of  disk  in  type  4 

V 

XI  PS 

I 

Moment  of  inertia     of  the  steering 

lb-sec  -in 

system  effective  at  the  front 

in2 

wheels  (Includes  both  wheels) 

G5j 

GN(5,J) 

I 

Cylinder  area  for  actuation  of 

trailing  shoe  of  brake  Type  1 

Kd 

FKD 

I 

Performance  parameter  charac- 

sec /in 

Vgiij 

GN(6,J)- 

i 

Brake  dimensions  for  type  3. 

in 

terizing  understeer/oversteer 

GN(ll.J) 

properties  of  the  vehicle 

612j 

GN(12,J) 

I 

Effective  lining  to  drum  friction 

KF,KR 

AKF 

I 

Front  and  rear  suspension  load 

lb/in 

coefficient  for  secondary  shoe  of 

AKR 

deflection  rate  in  the  quasi  - 

brake  type  3 

linear  range  about  the  design 

G130 

GN(13,J) 

I 

Mean  lining  radius  for  brake  type 
4 

in 

position  effective  at  the  front 
wheels  and  the  rear  springs 

G14j 

GN(14,J) 

I 

Coefficient  of  heat  transfer  for 
convectlve  losses 

KFC,KRC 

AKFC 
AKRC 

I 

Coefficients  for  the  compression 
bumpers  of  the  front  and  rear 
suspension  effective  at  the 

G15j 

GN(15,J) 

i 

Specific  heat  of  brake  assembly 

BTU/lb/°F 

front  wheels  and  rear  springs 

G16j 

GN(16,J) 

i 

Effective  weight  of  brake  assem- 

lb 

1                ) 

^C.^C 

AKFCP 

i 

Coefficients  for  the  cubic  terms 

bly  for  heat  absorption 

AKRCP 

of  the  suspension  compression 

hi 

HI  (4) 

Tire  rolling  radius 

1n 

bumpers 

!W 

FIDJ(2) 

i 

Drivel ine  Inertia  for  front  or 
rear  (Note  that  a  value  of  zero 
is  entered  at  the  non-driving 
end  of  the  vehicle) 

lb-sec2-in 

KFE,KRE 

AKFF 
AKRE 

I 

Coefficients  for  the  extension 
bumpers  of  the  front  and  rear 
suspension  effective  at  the 
front  wheels  and  rear  springs 
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"fe.Se 

AKFEP 

AKREP 

i 

Coefficients  for  the  cubic  terms 
of  the  suspension  extension 
bumpers 

Kp 

FKP 

M8 

Driver  steer  control   gain 

KRS 

AKRS 

1 

Rear  axle  roll-steer  coefficient, 
positive  for  roll  understeer 

KS1,KS2 

FKS1 
FKS2 

1 

Drivers  estimate  of  vehicle 
braking  and  accelerating  gains 

^Tl 

AKST(3) 

I 

Structural   hard  point 
spring  rates 

lb/in 

KT 

AKT 

I 

Radial   tire  rate  in  the  quasi- 
linear  range 

lb/in 

Kv 

AKV 

I 

Load -deflect ion  characteristic 
of  the  vehicle  structure 

lb/1n3 

Kfs,K*si. 

KiS2,^S3 

AKDS 
AKDS1 
AKDS2 
AKDS3 

I 

Coefficients  of  the  cubic 
representation  of  rear  wheel 
steer  as  a  function  of  de- 
flection for  independent  rear 
suspension 

V 

AKPS 

I 

Load-deflection  rate  for  the 
linear  steering  stop,  effective 
at  the  wheels 

lb-1n/rad 

Kl 

AKT 

I 

Slope  of  PR  vs  Pp  for  values  of 
Pp  between  P,   and  P2 

h 

AK2 

I 

Slope  of  PR  vs  Pp  for  values  of 
Pp  greater  than  P_ 

(LF)i 

FLF 

I 

Fade  coefficient  for  brake  at 
wheel   i 

MS 

XHS 

I 

Sprung  mass 

lbsec2/in 

"lIF'^JR 

XMUF 
XMUF 

I 

Front  (both  sides)  and  rear 
unsprung  masses.     Note 

Ml=M2=MUF/2>  M3=MUR 

lbsec2/1n 

MrM2 

XMUF 

Right  and  left  front  unsomnq 
masses 

lbsec  /in 

H3 

XMUR 

I 

Rear  unsprung  mass 

lbsec2/1n 

NBX 

NBX(5) 

I 

Number  of  x'   boundaries  supplied 
for  5  terrain  tables 

NBY 

NBY(5) 

I 

Number  of  y'   boundaries  supplied 
for  5  terrain  tables 

NDEL 

NDEL 

Number  of  entries  1n  the  front 
wheel  camber  table 

NX 

NX(5) 

Number  of  x'   grid  points  1n  5 
terrain  table* 

NY 

NY(5) 

Number  of  y'  grid  points  1n  5 
terrain  tables 

NZTAB 

NZTAB 

! 

Number  of  terrain  tables  entered 

ANALYTICAL 
SYMBOL 
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o 

z 

1 

DEFINITION 
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PrP2 

PONE 
PTWO 

I 

"Break"  pressures  for  brake 
system  proportioning  valve 

psig 

(RATI0)i 

Factor  used  to  modify  the 
nominal   tire-terrain  friction 
coefficient  at  wheel   1  to 
reflect  the  effects  of  vehicle 
speed  and  tire  loading 

RBB 

RBB 

2&> 

Constant  for  barrier  bottom 
plane 

1n 

RBi 

RBI 

2W 

Constant  for  barrier  face  plane 

1n 

RBT 

RBT 

20i 

Constant  for  barrier  top  plane 

1n 

RBl 

RBI 

m 

Constant  for  the  plane  perpen- 
dicular to  the  barrier  face 
plane  and  containing  the  axis 
of  rotation 

in 

NZ5 

NZ5 

i 

Flag  to  indicate  whether  the 
variable  increment  terrain  table  is 
supplied,  =  0,no,r10,  yes 

Z% 

SNPF 

347 

Roll  moment  acting  on  the 
front  axle 

lb-in 

2V 

SNPR 

bbO 

Roll  moment  acting  on  the  rear  axle 

lb-in 

*N«S 

SNPS 

»6 

Roll  moment  on  the  sprung  mass 
resulting  from  sprung  mass  impact 
forces 

lb-in 

ZNeS 

SNTS 

W 

Pitch  moment  on  the  sprung  mass 
resulting  from  sprung  mass  impact 
forces 

lb-in 

z"n 

SNPSS 

310 

Yaw  moment  on  the  sprung  mass 
resulting  from  sprung  mass  impact 
forces 

lb-in 

> 

8(8  • 

SNPU 
SNTU 
SNPSU 

457 
359 

Moments  acting  on  the  sprung  mass 
produced  by  forces  acting  on  the 
unsprung  masses 

lb-1n 

P.Q.R 

P.Q.R 

+8 

Scalar  components  of  the  sprung 
mass  angular  velocity  along  the 
vehicle  x,y  and  z  axes 

rad/sec 

PC 

PC 

I 

Hydraulic  pressure  in  brake  system 
master  cylinder 

psig 

PJ 

PP(2) 

197 

Hydraulic  pressure  In  brake 
cylinders  at  front  or  rear  brakes 

psig 

(PS) 

Prop  shaft  speed 

rpm 

10 
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(PT) 

XPS 

I 

Pneumatic  trail  of  front  tires 

In 

(TQ)E 

TQE 

110 

Engine  torque 

lb-ft 

RF,RR 

RF.RR 

I 

Auxiliary  roll   stiffness  of 
the  front  and  rear  suspensions 

lb/1n/rad 

;to)f(to)r 

TQF(50) 
TQR(50) 

I 

Front  and  rear  torque  tables  for 
a  single  wheel  and  effective  at 

lb-ft 

(RPHE) 

RPME 

in 

Engine  speed 

rpm 

the  wheel   (positive  for  traction, 

(RPS)1 

RPSI(4) 

44 

Rotational   velocity  of  wheel 

rad/sec 

negative  for  braking) 

i,  positive  for  forward  motion 

(TR) 

TTTR 

1 

Transmission  ratio  (speed  ratio 

of  the  vehicle 

of  engine  to  prop  shaft) 

RRM1 

RRH(4) 

352 

Rolling  resistance  moment 
acting  on  wheel  1 

lb-1n 

TR  «TR 
Kl     K2 

TESTRl 
TESTR2 

I 

Lower  and  upper  skid 
thresholds 

«W 

RW 

I 

Undeflected  tire  radius 

in 

TS 

TS 

I 

Distance  between  spring  mounts 

In 

RWHJB 

RWHJB 

I 

Beginning  and  ending  radii  for 

in 

for  a  solid  rear  axle 

RWHJE 

RWHJE 

calculation  of  the  radial   tire 

force-deflection  characteristic 

TSF 

TSF 

1 

Distance  between  spring 

in 

used  in  the  radial   tire  mode 

mounts  for  a  solid  front  axle 

SET 

SET 

I 

Ratio  of  permanent  deflection  to 
maximum  deflection  of  deformable 
barrier 

(TS) 

TTTS 

I 

Throttle  setting  expressed  as 
the  decimal  portion  of  wide  open 
throttle 

Si 

51(4) 

173 

183 

Total   suspension  force  for  a 
wheel,  acting  at  the  front  wheels 
and  rear  springs 

lb 

TsrTs2 

TESTSl 
TESTS2 

I 

Driver  threshold/indifference 
levels  for  positive  and 
negative  speed  errors 

in/sec 

(SLIP), 

SLIP(I) 

241 

The  amount  by  which  the  rotational 

(TYPE) 

NBTYPE 

I 

Brake  type  indicator 

speed  of  wheel   1   is  less  than  that 

of  free  rotation  expressed  as  a 
decimal  portion  of  the  speed  of 
free  rotation 

Tl? 

TIPS  I 

a 

Coulomb  friction  torque  in  steer- 
ing system  effective  at  the 
wheel 

lb-1n 

(SLIP^ 

SLIPP 

m 

The  value  of  (SLIP)1,  at  a  given 
wheel  center  speed  OV.,  for  which 

the  value  of/^.   Is  a  maximum 

T*» 

T2PSI 

it. 

Resistance  torque  produced  by 
the  front  wheel  steer  stops, 
effective  at  the  wheel 

lb-1n 

SPn 

ST(5,2) 

I 

Coefficients  for  straight  line 
segments  defining  the  desired 
path 

u.v.w 

u.v.w 

"3 

Scalar  components  of  linear 
velocity  of  the  sprung  mass  along 
the  sprung  mass  x.y  and  z  axes 

in/sec 

u',v' ,w' 

DXCP 

^ 

Scalar  components  of  linear 

in/ sec 

©Ji.CSJt, 

SI  I 

m 

Characteristic  lengths  of  inter- 

1n 

DYCP 

velocity  of  the  sprung  mass  along 

S2I 

Z85 

section  area  between  the  sprung 

DZCP 

the  space  fixed  x'.y'  and  z' 

(5,); 

S3I 

2fe 

mass  and  barrier 

axes 

t 

T 

Time 

sec 

ui,vi,  n 

UI(4) 

90- 
98 

Scalar  components  of  the  tire 

1n/sec    . 

T. 

TESTB 

I 

Braking  Indifference  level 

1n/s'ec 

VI  4) 

contact  points  linear  velocity 

b 

WI(4) 

along  the  vehicle  axes 

VTE 

TB.TE 

I 

Beginning,  ending  and  Incremen- 

sec 

UGi 

UG(4) 

lOi 

Wheel  center  forward  velocity  in 

in/sec 

TINCR 

TINCR 

tal  times  for  entry  of  control 

direction  parallel  to  the  tlre- 

tables   (TQ)F,   (TQ)R  andfp 

terraln  contact  plane 

VTR 

TF.TR 

I 

Front  and  rear  track 

in 

UGwi 

UGW(4) 

19? 

Ground  contact  point  velocity 

1n/sec 

T1 

TI(4) 

m 

Circumferential  tire  force 

lb 

along  the  circumferential  direct- 
ion of  the  wheel 

resulting  from  applied  torque 

U',n 

UNP(17) 

1SI 

Components  of  the  velocity  of  the 

1n/sec 

VTL       - 

TIL 

i 

Driver  steering  model  lag 

sec 

VNP  17) 

three  or  four  points  that  define 

TL 

and  lead  times 

•  ,n 

WNP(17) 

the  intersection  area  of  the 

(tqYbT 

(TQ)Dj 

TQB(4) 
TQ0(4) 

2H 

Brake  torque  at  wheel  1 

Drive  Hne  torque  at  prop  shaft 

lb-ft 
lb-ft 

n 

barrier  and  vehicle  along  the  spec 
fixed  axes 

!- 

at  vehicle  end  j 

11 
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UV!ST1 
U.ST1 
M  ST1 

UT 

vGi 


VGR 
VGR 
VGR 
VGR 
VGR 
VGR 


"43 
(VTAN) 

ME1 


(WOT) 

xB,xE 

XINCR 

xBB,yBB 
ZBB 


"BDRY 


BT 
% 
ZBT 


cpn 
ycpn 
zcpn 
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SYMBOL 


cpn 
cpn 


cpn 

^opf  *spt> 
^'gp; 


URP 
VRP 
WRP 

UPT(4) 
VPT(4) 
WPT(4) 

UT 

VG(4) 


VGR  12 
VGR21 
VGR23 
VGR32 
VGR  34 
VGR43 

VTAN 

WEIGHT(I) 

XIMPOR(I) 

TWOT 

XB(5) 
XE(5) 

XINCR(5) 

XBB 
YBB 
ZBB 

XBDRY(4,5) 

XBI 
YBI 
ZBI 


XBT 
YBT 
ZBT 


XCP 
YCP 
ZCP 

XCPN(3) 
YCPN(3) 
ZCPN(3) 

XCPNP(3) 
YCPNP(3) 
ZCPNP(3) 

XGPP(4) 
YGPP(4) 
ZGPP(4) 

XP'M] 

YP(4) 
ZP(4) 


Components  of  the  velocity  of  the    in/sec 
point  of  application  of  the  sprung 
mass  Impact  force  along  the  space- 
fixed  axes 

Components' of  the  velocity  of  the     1n/sec 
deformed  structural  hard  points 
along  the  space  fixed  axes 

Total  vehicle  velocity  in/sec 

Contact  point  lateral .velocity  In/sec 

in  the  direction  parallel  to  the 
tire-terrain  contact  plane 

Vehicle  speed  at  which  trans-  mph 

mission  upshifts  and  downshifts 
occur 


m 


iki 


n& 


m 


IV, 


Tangential  velocity  between  the 
vehicle  and  barrier 

Driver  steering  error  weighting 
function 

Driver  steering  error  importance 
weighting  function 

Wide  open  throttle  torque 

Beginning,  ending  and  incre- 
mental x1  for  terrain  tables 


Coordinates  of  the  intersection 
of  the  z'  axis  with  the  barrier 
bottom  plane  1n  the  vehicle  axes 

x'  intercept  for  angled  bound- 
aries within  terrain  tables 

Coordinates  of  the  intersection 
of  the  y'  axis  with  cutting 
plane  i,  in  the  vehicle  axes 


Coordinates  of  the  intersection 
with  the  barrier  top  plane  in 
the  vehicle  axes 


Coordinates  of  the  origin  of  the 
vehicle  axes  (sprung  mass  center  of 

gravity)  with  respect  to  the 

space  fixed  axes 

Coordinates  of  the  vehicle 
corner  n  1n  the  vehicle  axes 

Coordinates  of  the  vehicle 
corner  n  1n  the  space-fixed  axes 

Coordinates  of  the  ground 
contact  points  with  respect  to 
the  space-fixed  axes 

Coordinates  of  the  wheel  centers 
with  respect  to  the  space  fixed 
axes 


in/sec 


lb-ft 
In 


in 


in 


in 


in 


in 
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n_ 

X'pi 

y'pi 

XR1 


t£xR) 

(£y) 

(Sl5)J 

JSTI 
ySTi 
ZSTi 

XSTi 

ZST1° 


*VF 
*VR 

x'vp1 
y'vPi 

(y) 
{*) 

Ve 

YINCR 

YBDRY 
y'B 

y'cry'c2 

y,C3,y  C4 
y,C5,y'c6 

c 

z'bb 
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SYMBOL 


XNN(17) 

YNN(17 
ZNN(17) 

X 

Y 

XRI 
YRI 
ZRI 


SXR 
SYR 
SZR 

XSTI (3) 
YSTI(3) 
ZSTI (3) 

XSTI0(3) 
YSTI0(3) 
ZSTI0(3) 

XVF 


XVR 

XVP 

YVP 


X1.Y1.Z1 
X2.Y2.Z2 


VAR 
DER 


YB(5) 
YE(5) 
YINCR(5) 

YBDRY(4,5) 


YBP 


YC1P 
YC2P 
YC3P 
YC4P 
YC5P 
YC6P 


1% 


H& 


m 


DEFINITION 


YV 


ZBBP 


Coordinates  of  intercept  points   in 
between  the  barrier  and  sprung 
mass  in  the  vehicle  axes 

Coordinates  of  the  location  on 
the  desired  path  at  which  the 
1th  error  is  determined 

Coordinates  of  the  centrold  of    in 
the  intersection  area  on  cutting 
plane  i,  projected  on  to  the 
actual  vehicle  barrier  interface 
of  the  previous  time  increment 

Coordinates  of  the  point  of  ap-    in 
plication  of  the  sprung  mass 
impact  force 

Coordinates  of  the  deformed 
structural  hard  points  In 
the  vehicle  axes 

Coordinates  of  the  underformed 
structural  hard  points  1n  the. 
vehicle  axes 

Distance  from  the  sprung  mass     in 
e.g.  to  the  vehicle  front  along 
the  x  axis 

Distance  from  the  sprung  mass 
e.g.  to  the  vehicle  rear  along 
the  x  axis 

Driver  prediction  of  vehicle      in 
location  at  the  1th  sample 
increment  in  the  future 

Coordinates  of  accelerometer      in 
positions  with  respect  to  the 
vehicle  axes  for  which  acceler- 
ation components  are  output 

System  dependent  variable, 
integral  of  fy} 

First  derivatives  with  respect 
to  time  of  the  system  dependent 
variables 

Beginning,  ending  and  incremental  1n 
y'  for  terrain  tables 


Lateral  position  of  y'  terrain    1n 
boundaries  with  respect  to  space 

Lateral  position  of  the  barrier 
face  plane  with  respect  to 
space 

Lateral  positions  of  slope       fin 
changes  defining  a  curb 


Distance  from  the  sprung  mass 
e.g.  to  the  vehicle  side 

Elevation  of  the  bottom  barrier 
plane  1n  space 
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z'bt 

ZBTP 

I 

Elevation  of  the  top  barrier 

in 

<L«Pr*X- 

255 

Direction  angles  of  a  line  per- 

rad 

plane  in  space 

C<   '  CI'    C! 

pendicular  to  the  normals  of  both 

the  wheel   plane  and  tire-terrain 

Z'C2,Z'C3 

ZC2P 

I 

Elevation  of  curb  at  slopeC2 

1n 

contact  plane  with  respect  to 

Z'C4,Z'C5 
Z'C6 

ZC3P 

I 

Change  lateral  positions 

space 

ZC4P 
ZC5P 

V^i 

84 

Direction  angles  of  a  normal   to 
the  tire-terrain  contact  plane 

rad 

ZC6P 

_  . 

at  wheel   i  with  respect  to  space 

ZF 

ZF 

I 

Static  distance  along  z  axis 

in 

%A% 

ffe 

Direction  angles  of  the  resultant 

rad 

between  the  sprung  mass  center 

radial   force  on  wheel   i  with 

of  gravity  and  the  center  of 

respect  to  the  vehicle  axes 

gravity  of  the  front  unsprung 

masses 

iM 

145 

Direction  angles  of  a  line  from 

wheel   center  i   to  the  ground 

Z'              ZGPi    ; 

I 

Input  elevations  of  the  terrain 

In 

contact  point  of  tire  radial 

table  grid  points 

spring  j  with  respect  to  the 

Z'Gi 

ZPGI(4) 

126 

Ground  elevation  with  respect  to 

in 

vehicle  axes 

the  space  axes  of  the  point 

%-AiTii 

148 

Direction  angles  of  the  resultant 

rad 

beneath  the  wheel  centers 

radial   force  on  wheel   i  with 

Z'G1 

A  vector  through  the  ground  con- 

respect to  the  space  axes 

tact  point  normal  to  the  actual 

°iA>r>t 

257 

Direction  angles  of  a  line  per- 

rad 

or  equivalent  ground  contact 

pendicular  to  both  a  normal   to 

plane 

the  tire-terrain  contact  plane 

ZR 

ZR 

I 

Static  distance  along  the  z  axis 

in 

and  the  wheel   axis  with  respect 

between  the  sprung  mass  center 

to  space 

of  gravity  and  the  rear  axle  roll 

a  ,  A, rr* 

101 

Direction  angles  of  the  x  axis 

center 

with  respect  to  space 

ZVB 

ZBV 

I 

Distance  from  the  sprung  mass 

1n 

vA ,7y 

w 

Direction  angles  of  the  y  axis 

e.g.  to  the  plane  defining  the 

with  respect  to  space 

bottom  of  the  vehicle  along  the 

z  axis 

y^rW'j"'. 

95 

Direction  angles  of  a  normal   to 

the  wheel   1  with  respect  to  space 

Zirr 

ZVT 

I 

Distance  from  the  sprung  mass 

1n 

VT 

e.g.  to  the  plane  defining  the 

VA 

80 

Direction  angles  of  kingpin  axis 

top  of  the  vehicle,  along  the 

of  wheel   1 

z  axis 

>1 

Slip  angle  at  wheel   1 

rad 

at.    6   T 

B*    B*  B 

2W> 

Direction  angles  of  a  normal  to 

p'i  . 

BETP(4) 

m 

Equivalent  slip  angle  produced  by 

rad 

the  barrier  face  plane  in  the 

camber  of  wheel   1 

vehicle  axes 

01 

BETBR(4) 

113 

Non  dimensional   slip  angle 

OLjd     7 

277 

Direction  angles  of  a  normal  to 

variable  for  wheel   1 

BT'BT'BT 

the  barrier  top  plane  in  the 

(T2)t 

vehicle  axes 

GAM1 

47 

1 

<X  jO    r 

273 

Direction  angles  of  a  normal  to 

GAM2 

llnertial  expressions 

BTBl'lBl 

the  plane  perpendicular  to  the 

<Vt 

GAM3 

( 

barrier  face  plane  and  containing 

the  axis  of  rotation 
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(r4)t 
(r5)t 

<r6>t 

GAM4 
GAM5 

GAM6 

n 

\ Inertia!  expressions 

(r7)t 

GAM7 

( 

<rfi>t 

GAM8 

J 

(7-9)t 

GAM9 

DELBB 

Barrier  deflection 

in 

4 

DELI 

Right  front  suspension  deflec- 
tion for  independent  front 
suspension  or  front  axle  roll 
center  deflection  relative  to 
the  vehicle  from  position  of 
statis  equilibrium 

in 

4 

DEL2 

Left  front  suspension  deflec- 
tion relative  to  the  vehicle 
from  static  equilibrium  position 

in 

J3 

DEL3 

Right  rear  suspension  deflec- 
tion for  Independent  rear  sus- 
pension or  rear  axle  roll  center 
deflection  relative  to  the 
vehicle  from  static  equilibrium 
position 

in 

* 

DEL4 

Left  rear  suspension  deflection 
relative  to  the  vehicle  from 
static  equilibrium  position 

in 

A* 

DELG 

i 

Distance  between  road  roughness 
input  points 

in 

A< 

DELTA (4) 

m 

Distance  between  the  wheel   cen- 
ter and  ground  contact  point 

1n 

Af 

DT 

i 

Numerical   Integration  step 
interval 

sec 

&t3 

DELTB 

I 

Time  increment  for  use  during 
barrier  impacts 

sec 

Ai1 

DELTC 

I 

Numerical   integration  step  size 
for  curb  impact  option 

sec 

Atn 

DTR 

Integration  step  size  for  use 
with  wheel  spin  equations  of 
motion 

sec 

ATHF1 

ATHR1 

DTHF1 
DTHF2 
DTHR3 
DTHR4 

I 

Front  and  rear  half-track 
changes  with  suspension 
deflection 

In 

Ay'a 

DELYBP 

I 

Incremental   deflection  of  the 
barrier  position 

1n 

Av,j 

DPSILF 

H8 

Ideal  steer  angle  change 

rad 

-  -fa 

EPSB 

I 

Acceptable  error  1n  the  force 
balance  between  the  vehicle 
structure  and  barrier 

lbs 
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SYMBOL 
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«>•*« 

EPSF 
EPSR 

I 

Friction  lag  1n  front  and  rear 
suspensions 

1n/sec 

?n 

EPSL 

Permanent  set  of  the  barrier 
for  secondary  impacts 

1n 

Ev 

EPSV 

I 

Friction  lag  in  the  vehicle- 
barrier  friction  force 

1n/sec 

e* 

EPSPS 

I 

Friction  lag  in  steering  system 

deg/sec 

?e 

ZETAB 

I 

Threshold  value  of  wheel  rotation- 
al  velocity  below  which  logic  is 
applied  to  limit  brake  torques 

rad/sec 

Si 

J/1 

Suspension  displacement  of  the 
relative  to  the  vehicle  from  the 
position  of  static  equilibrium 

in 

Ccjn 

CDO 
CD1 
CD2 

Coefficients  for  unloading  force 
deflection  characteristic  of  the 
barrier 

% 

THESKD 

lib 

Vehicle  slip  angle 

rad 

Bei 

THGI(4) 

124 

Pitch  angle  of  terrain  under  wheel 
i  relative  to  the  space  axes 

rad 

®n 

THETN 

Value  of  0  at  t=0  or  at  the  nth 
indexing  of  the  axes 

rad 

fit 

THETT 

67 

integrated  value  of  e  from  t=0 
or  the  nth  Indexing  of  the  axes_ 

e;s,. 

k>t 

Angle  between  the  x  axis  and  the 
tire-terrain  contact  plane  at 
wheel   i 

rad 

!    xe 

TLAMB 

m 

Coefficient  for  inertlal  coupling 
terms  in  relationships  for 
driving  end  of  vehicle 

**** 

XLAMF 
XL  AMR 

i 

Ratio  of  conserved  to  absorbed 
energy  1n  the  front  and  rear 
suspension  bumpers  or  multiple 
of  Kp,K„  for  use  1n  simulating 

suspension  bumpers 

*r 

XLAMT 

i 

Multiple  of  K_  for  use  in  non- 
linear range  of  tire  deflection 

Mb 

XLM1(4) 
XLM2(4) 
XLM3(4) 
AMUB 

10  1 

<O0 
109 

1 

Contants  for  simultaneous  solution 
of  the  ground  contact  point 

Effective  coefficient  of  friction 
between  the  vehicle  sprung  mass  an 
barrier 

1 

Ma 

AMUC 

I 

Tire-curb  friction  factor 

14 
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/<*&[ 

XMUGI(4) 

i 

Nominal  coefficient  of  friction 
be.tween  tire  1  and  ground 

/*l 

XMUK4) 

T 

Peak  value  of  friction  coefficient 
for  side  forces  for  prevailing 
conditions  of  speed  and  load  at 
wheel  1 

MmL 

XMUM(4) 

I 

Nominal  test  surface  friction 
coefficient  on  which  tire 
properties  were  measured 

*-*-l 

XMUX(4) 

20 

Effective  friction  coefficient 
between  tire  and  terrain  at 
wheel   i  in  the  direction  along 
the  tire  circumference 

"*PL 

XMUXP(4) 

I 

Peak  circumferential  friction 
coefficient  for  tire  1 

M*5. 

XMUXS(4) 

I 

Sliding  circumferential   friction 
coefficient  for  tire  1 

K 

PI 

3.14159... 

s 

RHO 

I 

Distance  between  rear  axle  center 
of  gravity  and  roll  center, 
positive  for  roll   center  above 
e.g. 

Distance  between  front  axle 
center  of  gravity  and  roll   center, 
positive  for  roll  center  above  eg 

in 

Se 

RHOF 

I 

In 

S*i 

RHOS(I) 

1% 

Ratio  of  circumferential  to  peak 
side  force  friction  coefficients 
for  prevailing  conditions  of  speed 
and  load 

y^max 

RHOMAX 

m 

Maximum  value  of    si  at  the  exis- 
ting forward  velocity  of  wheel  1 

^P 

SIGR 

I 

Coefficients  for  the  polynomial 
form  of  barrier  load  deflection 
characteristic 

°T 

SIGT 

I 

Maximum  radial  tire  deflection 
for  quasl-Hnear  load-deflection 
characteristic 

in 

rA 

TAUA 

I 

Ambient  temperature 

°F 

*i 

TAU(4) 

Temperature  of  brake  assembly 

°F 

rCo 

TAU0(4) 

I 

Initial  temperature  of  brake 
assembly 

°F 

0Af 

PHIT 

THETT 

PSIT 

Euler  angles  of  sprung  mass 
axes  relative  to  Inertial  axes 

rad 

PMC  (50) 
PHIRC(50) 

I 

Table  of  front  and  rear  wheel 
camber  as  a  function  of  de- 
flection 

deg 

*cH 

PHICI(4) 

8b 

Camber  angles  of  wheels  relative 
to  the  normal  to  t1re-terra1n 
contact  plane 

rad 

PHICT,  PHIC2 
PHIC3,  PHIC4 
PHIC5,  PHIC6 

I 

Curb  slope  angles 

rad 
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<Pf 

PHIF 

Angular  displacement  of  front 
axle  relative  to  the  vehicle 
about  a  line  parallel  to  the 
x-axis  through  the  front  roll 
center 

rad 

% 

PHGI(4) 

85 

Camber  angle  of  terrain  under 
wheel  i 

rad 

0  i 

PHII(4) 

Camber  angles  of  four  wheels 
relative  to  vehicle 

rad 

*P'n 

PHIN 

Value  of        at  t=0  or  at  the 
nth  indexing  of  the  axes 

rad 

0* 

PHIR 

Angular  displacement  of  rear 
axle  relative  to  the  vehicle 
about  a  line  parallel   to  the  x  - 
axis  through  the  rear  roll  center 

rad 

<Pr 

PHIT 

oB 

Integrated  value  of  <<>    from  t=0 
or  the  nth  axis  indexing 

rad 

0yai 

Vi 

Angle  between  y  axis  and  tire- 
terrain  contact  plane 

rad 

^BDRY 

PSBDRY(4,5) 

I 

Angle  of  interpolation  boundaries 
in  terrain  tables,  measured  from 
the  x'axis 

rad 

Wf 

PSIF(50) 

I 

Table  of  front  wheel   steer  angle 
vs  time 

rad 

wt 

PSII(4) 

Steer  angles  of  wheels  relative 
to  vehicle  (positive-clockwise 
as  viewed  from  above) 

rad 

Vi 

PSIIP(4) 

m 

Steer  angles  of  wheels  in  tire- 
terrain  contact  plane 

rad 

K 

PSIN 

Value  of  f  at  t=0  or  the  nth 
Indexing  of  the  axes 

rad 

Vt 

PSIT 

v> 

Integrated  value  of  y  from  t=0 
or  the  nth  axis  indexing 

rad 

QF 

or 

OMEGF 
OMEGR 

I 

Maximum  suspension  deflections 
from  the  equilibrium  position 
for  linear  load -deflection 
characteristic  of  the  springs 

in 

nFC 
nRC 

OMEGFC 
OMEGRC 

I 

Front  and  rear  suspension 
deflections- at  which  the 
compression  bumpers  are 
contacted,  measured  at  the 
front  wheels  and  the  rear 
springs 

1n 

°FE 

OMEGFE 
OMEGRE 

I 

Front  and  rear  suspension 
deflections  at  which  the 
extension  bumpers  are  con- 
tacted, .measured  at  the  front 
wheels  and  rear  springs 

in 

°T 

OMEGT 

I 

Multiple  of  A2  at  which  the 
assumed  parabolic  variations 
of  small  angle  cornering 
and  camber  stiffnesses  with 
tire  loading  are  abandoned 

a 

OMGPS 

I 

Front  wheel   steering  angle  at 
which  the  linear  steering 
stops  are  engaged 

rad 
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3.         HVOSM  MATHEMATICAL  MODEL  ANALYSIS 

3. 1       Derivation  of  Equations  of  Motion 

The  deviation  of  the  equations  that  govern  the  general,  unsteady 
motion  of  the  vehicle  is  based  on  Euler's  equations  of  motion  written  with 
respect  to  a  moving  axis  system  in  order  to  avoid  the  appearance  of  the 
derivatives  of  the  moments  and  products  of  inertia  of  the  vehicle  sprung 
mass  in  the  equations.   These  equations  take  the  general  form  (see  for  example, 
Reference  1): 


w(w'  +  vP-ilQ)  =  Ft  (1) 

where  ( U  ,  V  >  W ) ,    (  P ",  Q_  >  R.   ) ,  (  L%   >  Ly  >  L£    ) ,  (  fy   ,    /^  ,  /^  ) ,  and 
(  A^x  »  kJu  '  Mi.   }   are  tne  scalar  components  along  the  vehicle  axes  of  the  linear 
velocity,  the  angular  velocity,  the  angular  momentum,  the  external  forces, 
and  the  external  moments,  respectively. 

Applying  the  assumption  that  symmetry  exists  in  the  Y-ifi  and  ^2. 
planes,  i.e.,  the  Y-ty  and  HZ.  products  of  inertia  are  zero,  (Z)o/=  £yt  =  0) 
the  angular  equations  become: 


(2) 


ly  a  -  Czx:P)iM  -  £?(!*  -jx)  =  sly  (5) 

(4) 

-ly.tP  +  ufc-paCTt-iyi+azitz  =  ^2 
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3.1.1      Sprung  Mass 

In  the  following  derivations,  inertial  forces  produced  by 
accelerations  of  the  unsprung  masses  are  treated  as  external  forces  acting 
on  the  sprung  mass.  That  is,  the  Euler  equations  of  motion  for  the  sprung 
mass  are  modified  to  include  inertial  effects  of  the  unsprung  masses. 

For  the  case  of  sprung  mass  accelerations  in  the  X-  direction,  each 
of  the  unsprung  masses  is  assumed  to  act  as  a  point  mass,  With  this  assumption, 
the  equating  of  inertial  and  applied  forces  in  the  X  direction  is 

M5r#-Tr£*.*r<2)  -  Zf**  '+  2fK6  -  (2M)gsin0  -2Fxxi  (5) 

where  %iiy.L   are  the  unsprung  mass  inertial  forces  in  the  X   direction. 

The  scalar  component  in  the  X  direction  of  the  acceleration  of  a 
particle  in  a  moving  frame  of  reference  (e.g.,  see  Reference  1)  is 

a?;   =  U  -  vfi  +■  urQ   +    X  +■  2  Qj    -  Z  Ry 

-%(4*+X*)+y(PQ-R)   +    }(e*  +  Q)  (6) 

In  the  vehicle  representations  depicted  in  Figure  3.1,  constraints 
on  the  suspension  system  preclude  motions  of  the  unsprung  masses,  relative 
to  the  sprung  mass,  in  the  y. .  direction.  Therefore,  in  the  present  application 

00 

of  Equation  (6),  )C    =  0.  Also,  dynamic  effects  of  motion  of  independently 
suspended  unsprung  masses,  relative  to  the  sprung  mass,  in  the  M   direction 
are  assumed  to  be  negligible.  Application  of  (6)  to  the  cases  of  the 
independently  suspended  unsprung  mass  centers  of  gravity  yields  the  following 
expression  for  inertial  forces: 

Fixi  =  Mi#*i~  MLl\X-vg+ Ufa  +<2jgl$l  -  x(aL+e2)+y(Pa-£)<-i(p£+&)l  a) 
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18 


where 


X-,    =  *- 

y^TF/z 

£,=    2f^/ 

)CZ  =■  &, 

i/^-V^ 

4  -  zf  '4 

^3  =  -l? 

#a  s  r*  A 

£3  =  2z+J3 

Af  -.-£ 

^f  =  -fc/Z 

fy    '    ?£   ^c£f 

and  ^L     =^VF/Z   for  I  =  1,   2  or  Mi.    =  ^ud/l    for   ^  =  3,  4. 

Similar  application  of  (6)  to  the  beam  axle  centers  of  gravity 
yield: 

^l  =  M-LaH  *  M£yL-v-R.twCi  +  ZQiL+iiijL  -  U&+&2)  +y(pci-i)+z(F&-ci)l  m 

where 

*,  -  cl.  y,=  -pr<f>F  £,=  2f+S,t-/>F 

X3  =  ^  2/3  -  -/^  £3  =  %+4  V 

and  Alt   =  A/fi//=  for   i.  =  1  and  ^  =  Mt^e   for  i   =  3. 

Substitution  of  (7)  and/or  (8)  into  (5)  yields  the  equation  of 
motion  for  acceleration  of  the  sprung  mass  in  the  /.  direction  as  shown  in 
the  first  row  of  the  l\0\\      and  \\E\[     matrices  in  Section  3.1.5. 

The  equation  of  motion  for  acceleration  of  the  sprung  mass  in  the 
M      direction  is: 

M*tir+uZ-vrP)  =  £FyU.i-EFy«f  +  (2M)ga>s0si*</>-2fX;fL      (9) 
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The  scalar  component  in  the  u      direction  of  the  acceleration 
of  a  particle  in  a  moving  frame  of  reference  can  be  expressed  as: 


a^=  ir  +VLR-TArP+ij  +  2PJL -2Pt  +  x(P$  +  P) -y(P*+Pz)±z(QR-p )    (10) 

Independently  suspended  unsprung  masses  are  assumed  to  act  as 
point  masses,  and  for  purposes  of  dynamic  interactions  are  also  assumed 
to  be  constrained  to  motion  parallel  to  the  £  axis.   For  this  case,  the 
unsprung  mass  inertial  forces  acting  on  the  sprung  mass  are: 

%i  =  Miayi  ^M^ir+TxP-i^P^P^^x^P^PJ^^P^^^^Ca^P)^1^ 

» 

where 

H=-b  yi~'rz/2-  23  =  ^/4 

and   M^   =  /Vlt/f/2   for  i    =  1 ,  2  or  d/fc"  =  Mug/2    for   ^' =  3,  4. 

Inertia  forces  for  a  beam  axle  suspension  are  approximated  by 
means  of  a  thin  rod  representation  as  illustrated  in  Figure  3.2.   For  this 
case,  the  inertial  force  is  given  by: 


fxijL         *j 


-Hi  -!■ 


2) 
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Let    dm    =   X  d  £ 


Figure  3.2     THIN  ROD  APPROXIMATION  OF  REAR  AXLE  INERTIAL  PROPERTIES 
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and 


X1    =    a  y,  =  -pFsin.<^  -  Hco$(pF  £,    =   iF^£,^pFc.oS<pp-is^^F 

ty  -  -pF<pF  cosfr  +l<pF  5LncpF  lf  -    Jj  -pF<pF  si'nfc  -Z<pFC0S<t>F 

$  =   /V  &2  siifc  +1  <p* cos <pF 

~PF  Cpp  COS  <pF    +Q-4>F  SLYV<f>F 

ifs  —  -p$gCoS<frfz+JLcj>£'5>u<i<j>£  ^  -    &3-pcJ>eSM>(f>g-JLfac0S<l>£ 

-p$z  cos  fa  +£$*  sin.  fa 


Substituting  into  (12)  and  assuming  $>p    and   <$e  to  be  small 
angles  yields: 

Fkj,   =  MUP[v+uR-KP+j>Fch:$  -fcfc-lPf,  +2PpF4>f4>F   +ol(PQ+&) 

-t-f>F<t>F(P2+&2)-h(2F.i-Sl  +pp)(ClZr  P)~]  (13) 

hp  *  Mvi&v+iLRr1i>P+p42  4£-pfa-2PS3  +  2Pjo4t4k   -b(pa+iL) 

+p<tk(Pl+IP)  +   tez+4*-/>XQ&-P)  ]  (14) 

Substitution  of  the  appropriate  inertial  force  equations  into  (9) 
for  the  different  suspension  options  yields  the  equation  for  vehicle  accel- 
eration in  the  y     direction  as  defined  by  the  second  row  of  matrices  //£>// 
and  //w/   in  Section  3.1.5. 

For  the  case  of  sprung  mass  accelerations  in  the  £,  direction, 
interaction  with  the  unsprung  masses  takes  place  through  the  mechanism  of 
the  suspension  forces.   Therefore,  the  corresponding  equation  of  motion 
can  be  written  directly  as 
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Ms(ur+  vP-uQ)  =  -ZFv  +  L^S  +  M*  9  t°s  &  ccs  $  (15) 

This  equation  appears  in  the  third  row  of  matrices  \\J)\\     and  //£||  in 
Section  3.1.5. 

The  equations  defining  angular  motions  are  derived  in  a  similar 
manner.  That  is,  the  preceding  definitions  of  the  unsprung  mass  acceleration 
components  are  substituted  in  the  following  equations. 

-  Z^£Lq  CD50SM  <j>  (16) 

where  riyi     are  defined  by  (11)  or  (12)  for  either  independent  suspension  or 
beam  axle  suspension  unsprung  masses.   For  independent  suspension  £i   are 
given  by: 

i,  =   2f  tS[ 
2-2.  ~    ^/=  ~t~&z 


(17) 


and  for  beam  axles 


2/  =  zf  1-4 

2a   =  Zfi+<f3  (18) 


where  f£x£  are  given  by  (7)  and  (8),  ^-i     by  (17)  for  the  independent 


(19) 


suspension  and  by 
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for  the  beam  axle  suspension  option. 

-  ^^  *<:  ■*  2?M<:  Xigcosesi'/H}  +  J?/V1;  yL#*mO  (20) 

where  fr^i  ,   F^w   are  given  by  (7),  (8),  (11)  and  (12)  for  either 
suspension  type,  ^  =  <x    for   6  =  1,  2  and  X-i  =  ~t>   for   6  =  3,  4. 
For  independent  suspension, 

y,  =  tf/% 

yz  =  -  t>a 

Jk  =     ^ 

ft  --  -  W2 

and  for  beam  axle  suspension 

y,  -  -A  fr 
^^--  ~p4z 

Appropriate  substitutions  lead  to  the  fourth,  fifth  and  sixth 
rows  of  matricies  |PJ|  and   //£//  in  Section  3.1.5. 

3.1.2      Unsprung  Masses 

In  the  following  derivations,  the  general  procedure  that  is  applied 
to  aircraft  control  systems  in  Reference  1  is  followed.   That  is,  the  equations 
of  motion  are  obtained  by  application  of  Lagrange's  equation  of  motion  in  a 
moving  frame  of  reference: 
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=  ?  (21) 


*t  \*tJ        d$K 


where  T     =   kinetic  energy  relative  to  the  vehicle-fixed  axis  system  of  that 
portion  of  the  total  unsprung  mass  that  has    (L^      for  a  degree  of  freedom, 


ZF"   =■ =  generalized  force,  or  moment, 

W     =   work  done  on  the  unsprung  mass  system,  in  degree  of  freedom  fyk      , 
by  the  external  forces  which  act  upon  it,  including  the  effects  of  the 
inertia  force  field  produced  by  acceleration  and  rotation  of  the  vehicle- 
fixed  axis  system,  and  fa    =  generalized  coordinate  (i.e.,  degree  of  freedom) 
relative  to  the  vehicle-fixed  axes. 

Independent  Suspension  Unsprung  Masses 

The  kinetic  energy  of  an  independently  suspended  unsprung  mass 
relative  to  the  vehicle-fixed  system  (in  its  assumed  form  of  a  point  mass 
and  neglecting  wheel  rotation)  can  be  expressed  as 


(22) 


where  the  subscript  I     refers  to  the  mass  number  (1  =  right  front,  2  =  left 
front,  3  =  right  rear,  4  =  left  rear)  and  the  subscript  J   identifies  the 
front  or  rear. 

The  generalized  force,  <r   ,  is  obtained  by  means  of  determination 
of  the  virtual  work  done  by  external  forces,  including  inertia  forces,  during 
a  virtual  displacement  of  the  right  front  unsprung  mass: 


StV  =  (f>uL  +  SL  +^-2  cos  0  cos  0jS(Si)  +  <5(Wi) 


(23) 
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where   ™i      =  work  done  by  the  inertia  forces, 

Let  the  acceleration  of  the  right  front  unsprung  mass  relative  to 
the  space-fixed  axes  be  a.    ,  and  relative  to  the  vehicle-fixed  axes  be  £■ 
Then  the  acceleration  due  to  motion  of  the  vehicle  is  (  &  -  cL' )   and  the 
corresponding  inertia  force  is 


dFt'  =  -(&-&')  dm 


(24) 


The  components  of  df[        in  the  X     and  y      directions  are 
perpendicular  to  the  virtual  displacement,  S(Si)    ,   and  can,  therefore,  do 
no  work  during  the  displacement.  The  scalar  component  of  dFp       in  the  J 
direction  can  be  expressed  as  follows: 


dF>.  =  - 

9c 


w-+  y*P-  uQ  +  Xi(PK-Q)  +    yi(QX  +  p) 


-(fj  +Si)(r*+Qz) 


dm 


(25) 


Since  the  front  independently  suspended  unsprung  masses  are  treated 
as  point  masses,  the  virtual  work  due  to  inertial  forces  is  given  by 


SWi   = 


•  M, 


u-1 


w  +  vP  ~-uQ  +  H(PX-Q)  +   p(QK  +  f) 


(t*  +  *i)(r* +-<**)]*  (*i) 


(26) 


Application  of  Equations  (22),  (23)  and  (26)  in  Equation  (21)  yields 
the  equations  for  independently  unsprung  masses  in  the  matrices  Hull    and  jlE.ll 
in  Section  3.1.5. 
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Beam  Axle  Unsprung  Mass 

From  Figure  3.3,  the  motions  of  the  center  of  gravity  of  the  axle 
are  defined  by  the  following  relationships  (note  that  the  subscript  j 
is  used  here  to  indicate  either  the  front  or  rear  axle) : 

y  =  -/>y  siajj  i  =  2/  +Sf  +j>j  ct>s<pj 

y =  -  Pi  h  ^  to  i=  %  -/>/  &  sth  $/ 

y  -  />j(4j  5'"^/  - <j>j a%<j>j)       £  -  Sj -J>j(4>j  cosfy  +4>j  sirup/) 

The  kinetic  energy  of  the  axle,  relative  to  the  axis  system  fixed 
in  the  sprung  mass  can  be  expressed: 

r-   |  Iji)  +   |  Afc//>/#  +  Jf-.2/>j  ¥/*/'*>&)  (27) 

Evaluation  of  the  left  side  of  Equation  (21)  for  each  of  the  axle 
degrees  of  freedom,  and  a  subsequent  assumption  of  a  small  angle  restriction 
on  <j>j      ,   yields: 


J/2T\  bT 


Jt[    d4>j/  dcpj 


-  (Ij  +MUj'ftj)  ff  -  (Mujftjj)  &}       (28) 


The  work  done  by  the  inertia  forces  produced  by  vehicle  motions 
during  a  virtual  displacement,    <f  (  _.<£/  ) ,  is  determined  as  follows  (see 
Figure  3.3  and  Equation  (22). 

*W;        r  ^  I     *i  ,-      »\  d$  (30) 


i)    J-T,/z  J-7J/z 
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Figure  3.3  AXLE  REPRESENTATION 
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The  scalar  component  of  acceleration  of  element  aL      in  the  2 
direction,  due  to  motion  of  the  vehicle,  is  defined  by  the  following 
expression: 

a«-  'j  =  -Cs  +  Pv-Qa-2P/oj$}.  cos  $j  +  x(PK-q) 

-f/Oj  sen  0j  +  J  cos  fy)  (Q#+  P) 

( 311 

Substitution  of  (31)  into  (30),  integration  of  (30)  where 
%   ■=■    (X       for  the  front  and  X-      =  -b     for  the  rear,  and  subsequent 
assumption  of  a  small  angle  restriction  on  (Pj      ,  yields  the  portion  of  the 
generalized  force  that  corresponds  to  inertial  forces  produced  by  vehicle 
motions.  When  this  results  and  Equation  (29)  are  substituted  in  the  following 
relationship,  the  equation  of  motion  for  the  vertical  acceleration  of  the 
beam  axle  is  obtained: 


+  K~.  +'  f* 


d  far]   bt     s(Wj) 

+  Si  ■  +  Sitl  +  MUJ  g  cos  e  cos  <?>  (32). 


The  work  done  by  the  inertia  forces  produced  by  vehicle  motions 
during  a  virtual  displacement,  S(<Pj)    >   ^s  determined  as  follows  (see 
Figure  3.3  and  Equation  (24)  ). 

Y($~) =  ~J     cF* (pj€0s % ~ £ sin **)  * dFffas<n %  ***** *j )] 

T//z 

+  (**  -  *)(/**'*  *J+£  cos  *j  \d£  (33) 
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The  scalar  component  of  acceleration  of  element  cf£       in  the  Z. 
direction,  due  to  motion  of  the  vehicle,  is  defined  by  Equation  (31).   The 
corresponding  component  in  the  it       direction  is  defined  by  the  following 
expression: 

a^ '  -  y  =  v-  i-  vf u  -Pw  -2P  (Sj  -  /Ojjtj  sin  fy  )  -     (PQ  +6) 
+  (/OjSJn  gSj    +4  cos  fy)  (P*+  /? *) 
+  (%j  +*i  +/*J  ™s  tfj  -4 Stn  0j)(QK-f)  (34) 

Substitution  of  (31)  and  (34)  into  (33),  integration  of  (33)  and 
subsequent  assumption  of  a  small  angle  restriction  on  <^/   ,  yields  the 
portion  of  the  generalized  moment  that  corresponds  to  inertial  force 
produced  by  vehicle  motions.   When  this  result  and  Equation  (28)  are 
substituted  in  the  following  relationship,  the  equation  of  motion  for  the 
angular  acceleration  of  the  beam  axle  is  obtained. 


'3T  )        dT  S(Wj) 


+  /V  (35) 


3.1.3     Steering  System 

The  steer  mode  degree  of  freedom  is  defined  by  the  following 
relationships  with  the  assumption  that  inertial  coupling  effects  are 
negligible: 

lfrYF  +  TiY  +  T2y  =    Fyu.1  ( kj  cqsolKa  -FTcosti  cas&w) 

+      fytf,  (kiCOSotu  -  PTc06-fz   COS  QX3t ) 

-    FXU1  krfcosp^+fa cos \,)  -  fxai kz (ccxfth2 -t  4>i  cos rkz ) 
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where 


Tiy   = 


0     ,  faiW&er 

ClySQlff    ,  for  lfFl>Gf 


T     _  J   0      ,    fa  WfI  4JLf  ,ors0nt?s&n,fF 


The  assumed  form  of  the  steer  angle  stops  is  depicted  in  Figure  3.4 


STEERING  SYSTEM 
RESISTANCE 


k- n^-* 


-n 


— jijr— j 


FRONT  WHEEL 
STEER  ANGLE,  JTf 


Figure  3.4    STEERING  SYSTEM  RESISTANCE  VS  STEER  ANGLE 
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In  the  HV0SM-VD2  program  version,  torques  generated  within  the 
steering  system  by  gyroscopic  precession  of  the  spinning  front  wheels  are 
approximated  by  means  of  the  following  added  term  in  Equation  (36): 


4    • 


+ 


*   *    * 


(37) 


3.1.4 


Wheel  Spin  (HV0SM-VD2  Program  Version) 


The  inertial  coupling  of  the  drive  wheels,  through  the  differential 
gears,  and  the  associated  effects  of  drive-line  inertia  can  exert  a  significant 
influence  on  the  occurrence  of  wheel  lock.   Therefore,  it  was  necessary  to 
derive  differential  equations  for  the  rotational  motions  of  the  drive  wheels, 
that  would  include  these  effects.   Since  the  developed  program,  for  generality, 
permits  the  simulation  of  either  rear  or  front  wheel  drive,  the  derived 
equations  are  applied  to  both  ends  of  the  vehicle.   The  coupling  effects 
are  suppressed  at  the  nondriving  end  by  means  of  the  input  of  a  zero  value 
for  the  associated  drive  line  inertia. 


Differential  6ea^5 

OF  RATIO  64R)/ 


Figure  3.5        INERTIAL  COUPLING  OF  DRIVE  WHEELS 
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The  kinetic  energy  of  the  system  depicted  in  Figure  3.5  can  be 


expressed 


r  =  { inOMfy^ui  (o*&  i  Jdj  mfoj 


(38) 


Let  (RPSL    =     (AfOj 


~(m)j  t »; 


then   p*  \  l^(m)[Al^mU  ^  V^/ 


r  (&&+  2m\  M*  t  (&$?« 

4 


(39) 


Collecting  terms,  Equation  (39)  becomes 


r= 


Lz         a     . 


(ffi)f  + 


(zps)^  +  1-^^i.m\(m)^  (40) 


The  potential  energy  of  the  system,  V,  in  the  rotational  mode  (i.e., 
torsional  wind-up  of  axle  and  propeller  shafts,  elastic  distortion  of  tires) 
is  assumed  to  be  negligible.  Therefore, 

V  =0  (41) 

Energy  dissipation  may  be  expressed  in  the  following  form: 

D=  -12(TQ)DJ(R?S)D)+  [FcckL^rl2(T0)Bl  W\  4-  %LHK^in  -tfft^fffSk,     (42) 


Substitution  for 


and  collection  of  terms  yields 


(43) 
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Application  of  Langrange's  equation  of  motion  to  Equation  (40), 
(41)  and  (43)  yields: 


Iwi    +■ 


1dj(A$j 


TyCAtij 


4 


l 
Jit 


Cm\ 


1„-(A&; 


4 


\(m\„ 


=  -Fciki  +  nCTCdsi  +t?(Af0j(TQ)Dj  4-Kmi 


(44) 


J 


4    \yi 


cm, 


4- 


\S  cmh\ 


(45) 


--&U1  ^  +\2(TQ)3U1  +  6>(A&j(TQ)Pj-  *RMCt1 


The  rates  of  change  of  the  wheel  angular  velocities  are  integrated 


directly. 


(RPS)L  =  (RPS)Lo 


ft(RPS)L    Jt 


(46) 


where 


(m)u  = 


ae,i<>ux>fi0  +r&iDSLrifi0 


k 


ID 
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3.1.5     Summary  of  Equations  of  Motion  for 
Sprung  and  Unsprung  Masses 


3.1.5.1   Independent  Front  Suspension  -  Berm  Axle 
Rear  Suspension  Equations  of  Motion 


The  equations  of  motion  are  simplified  by  the  use  of  the  following 
combinations  of  inertial  terms: 


ZM 

(OX  - 
(%\  - 
Wt  - 

#>)«  - 

>t    - 

Mt  ■ 
(h)t  - 

(>A  - 

(tit  " 

(h\  • 
(*f)t ' 
(*$)t  - 
(h)t  - 


=   Ms  +  M<je  +Ml 


UF 


'ax 


=     Al. 


UF 


ft  +  ^fa+Sz)*^-^ 


+  Mu*(?«  *<6)**/*^*x  (}*  fS3) 


(47) 


(**)<  +  />  Mu/t  fa  +  S3  +/>) 

Mur*  (if  +  ^J")-***  h  fy  *Ss) 
^j-t  (SrSz)  -M„  />  <t«  (*«  +  />+*) 

M«*(*'+*L^sy+  *«*(**  +<°  +**) 

(s,  +4)  +  ***  (*3  -/>  **  ?«)] 
z  \^l!  (if  -i2)  - /fM/?  ^  (i3  ~/o  fa  to) 

z  \^a.(s,+Sz)-M««Jb(S3-fitoi>«) 


z 
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3.1.5.2   Independent  Front  and  Rear  Suspension 
Equations  of  Motion 


Xy     ~    X/ 


li   =■  Mup 


W+(W   +  Mull* +(%)*] 


y,    =     MuFO.-Muzb  C49) 

&  -    MvF(4+h)    +  Muz  (6^64) 


Inertial  Term  Definitions  for  Independent 
Rear  Suspension  Option 
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Matrix  Equations  of  Motion  for  Independent 
Rear  Suspension  Option 

39 


+ 


<3 

>^ 

i 

ex. 

>^ 

i 


+ 

I 

c*: 


V<5 


«£■ 

+- 

u* 

+■ 

a. 

+ 

OS" 

>^ 
/ 

s 

I 

.§ 

+ 

pi 

1 


^ 


4- 
Q, 


4- 

+ 
1 

5 

-t- 

H 

1 


c* 
*■ 


H 

k 

Hi 


Ob 

9 


-v- 


H 


§ 


H 


■+- 
Of 

Ho 

H 


H 

+ 

H 

H 

1 

I 

1 


•+- 

-+- 
a. 


S 


+ 
4* 

V? 

+ 

til 


4- 
I 


4 


U? 


-V- 


cc 


3m 


•si 
V/Y 

■V- 


-t- 


1  N 

•+ 


VP 

^ 

v^ 

tf 

s 

$ 

VA 

0 

8 

+■ 

+ 

cx 

Cc 

^ 

is 

<?>r 

^ 

^ 

?» 

|n 


j3H 


-t- 

+ 

8 

Q 
£* 

>= 

I 


Mi 


4- 


VSfIN 


+ 

%/> 

3 

VA 

CL 

b 


W 


uu 


40 


3.1.5.3   Beam  Axle  Front  and  Rear  Suspension 
Equations  of  Motion 


X'i  =     huF(<?+0fP )     +Muz( k?+pZfl ) 
r'xi  =     Mura(zFi-Sf)    -  Mwb (%&+<£&) 

l'yi  =  -MuFftfi  (zf+Pf  +*f)  ~Muefe fa (?£+/>£  +<&) 
y,   =       MufCL  -Muzb 

Y3  =      2T2  -  Mu/?/)/2   -Mvf/f 

|f6  =  2  M*  ^  iF  -f>r$FJF)  +  2-  MuJZ  (Jj*  -pzfo  fa) 

fe=  -^^PFJF^F'M^^2^^h^-fi^h) 


Inertial  Term  Definitions  for  Solid 
Front  Axle  Option 
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Matrix  Equttions  of  Motion  for  Solid  Front  Axle  Option 
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3.2 


Vehicle  Position  and  Orientation 


3.2.1 


Sprung  Mass 


The  transformation  matrix  from  the  coordinate  system  fixed  in  the 
vehicle  sprung  mass  to  the  inertial  coordinate  system  according  to  the  rotational 
sequence  TfT      ,       0     and  <P     is: 


COS  9  cos  1/r   i  -  COS  $  sin  Tfr  +  sin  $  sin  9  cos  f  i  Sen  j>  sin  f  +  cos  fl  sin  9  cos 


i! 


COS  9  sin  yr  i     cos  $  Cos  p  ■*  sin  0  Sin  0  sin  f  \-cosyr  sin  4>  +  cos  $  sin  9  sin 


-sen 


cos  0  sin  4> 


cos  9  cos  & 


where        0'    =  £?t  ,        q\      =   fa        ,   and        njf     =  lj/t 

In  order  to  obtain  the  Euler  angles  required  it  is  necessary  to 
integrate  the  components  of  the  rotational  velocities  (   P  ,  Q.     ,    £.      ) 
along  the  axes  of  rotation  of  the  Euler  angles.   These  components  are  given 
by: 


fl  =(a<SLn4>t  +  £cot4>l)sec&i 

Note  that  the  presence  of  the  tan  Q^       and  sec  0^      in  these 
expressions  result  in  a  singularity  at  Q^.     =  7L/2. 


(54) 
(55) 
(56) 


(53) 


An  indexing  scheme  that  allows  unlimited  angular  excursions  of  the 
vehicle  involves  the  use  of  three  sets  of  Euler  angles: 

(1)       ljrt    t       &±     >    4-t        are  tne  Euler  angles  that  define 

the  orientation  of  the  vehicle  with  respect  to  the  inertial 
axis  system  in  which  the  X/  and  u'     axes  are  horizontal 
and  £'     is  vertically  downward. 
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(2)  \ft       ,  Ot       >    fit         are  the  Euler  angles  that  define 

the  orientation  of  the  vehicle  with  respect  to  an  intermediate 
system  (the  indexed  axis  system) .   Note  that  this  axis 
system  is  fixed  in  space  either  coincident  with  the  initial 
orientation  of  the  vehicle  or  is  oriented  colinear  with  the 
vehicle  axes  at  the  time  of  most  recent  indexing. 

(3)  i^n     ,      Q'n      >    (p'/L     are  tne  Euler  angles  that  define  the 
orientation  of  the  indexed  axis  system  with  respect  to  the 
inertial  system. 

Thus,  it  can  be  seen  from  Equations  (54),  (55)  and  (56)  that  the 
rates  of  change  of  the  Euler  angles  that  are  calculated  are  written  with 
respect  to  the  indexed  axes  and  integrated: 


e±  -  eit+  f&Jt  (57) 

<K  =    <fi'n+f<ptdt  (58) 

Vt=r*+J%cli  (59) 


These  angles  are  then  related  to  the  vehicle  Euler  angles  with 
respect  to  the  space  axes  by  the  following  logic. 

If  O0    -    (j>0       =0  and  coordinate  system  indexing  has  not 
occurred, 

Qt--ei 
%  -  tl 
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If  &o      f   0,  <po       f   0,  or  coordinate  system  indexing  has 
occurred, 

(1)  For  first  time  increment  of  a  given  run  only,  set 

(2)  Let    llBn.ll        =    II  All  with  f    =t'n     ,      9    =  $U      , 

<fi      =     (}>  n,       (note  that  ll£>n,ll      remains  constant 
between  times  of  indexing) . 

Let  HCn.ll    =    M        witn  f   =  fl      ,  9   =  O't.     > 

<p     =  (j>'t 

Therefore,  _jl3n.ll    transforms  a  vector  in  the  indexed 
coordinate  system  (  V        )  to  the  space  fixed  coordinate 
system 

V'=llBn.Uf" 

and  llCn/l    transforms  a  vector  in  the  vehicle  coordinate 
system  (  y-     )  to  the  indexed  coordinate  system 

vr"-  IIUilv- 
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Therefore 


_  / 


V   =   IIBnlHICnllf 


and,  by  definition, 


II A 11=  II  5»ll -HUH 


(3) 


(4) 


By  inspection  of  the  elements  of     ilAlj       ,    (Equation   (53)    ), 
it   is   seen  that: 


tan  tt      =    Ait/A,, 


Therefore, 


ift  =  ardan 


&2lCn   +  322^21  +■  ^Z3^31 


Number  of  revolutions   is  integer  value 

1V     2^ 


(60) 


(5)    Define  four  possible  values  for  testing 


TRYPSI. 


TRYPSI. 


TRYPSI, 


TRYPSI 


%       +2%(IV) 

ft     +1b+  Z1t  CIV) 
%    -It  +  Z1t(iv) 
tt    -ZK+  Zfl  (TV) 
ft    +Ilt+  ZTt(ZV) 


for  positive  v4  , 
for  negative  rt     . 
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(6)  Test  each  of  the  four  possible  angles  to  seek  correct 

one.      For         I    =   1,    2,    3,    4 

DlFFt   =  TRXPSIi  -  ft-, 

(a)  If         IP£FF[I*£     ,   TRYPSI.    is  the  correct 
one.      SET      ^     =  TRYPSI.  . 

(b)  If  none  of  the  four  values  of     jV^FFcl^   %       t 
let    rV    =ft,,/'2lt    +      i    sgn   (^/flt  ). 
Return  to    (5) . 


(7)  (a)      For       lifotfi  I       0. 


07, 


Qi  =  arcton 

(b) 

Qt  -  arctan 


-cosWtfdaiCn+faCii+BmCt,) 


For0.7<      ISLnlfel       . 

-6ln%(B31C11+BnCl1+333C3ij 


(8) 


A  =  arctan 


&31  ClZ  +  &32  ttz  +  83$  C51 


531  0^  +  532  Cz3  *-  3^3  Ci3 


(61) 


(62) 
(63) 


(9)    Repeat  steps  (4)  through  (6)  for    €>t  ,     (fit 

The  values  of  the  angles   1ft     >        @-b    »  <fit       at  the  time  of 
indexing .  are  saved  as  f^n,    ,       &/i    ,       <r*t-  ,  and  the  angles   ft   ,  Ot      , 
and  Yt.       are  set  to  zero. 

The  position  of  the  sprung  mass  center  of  gravity  is  determined  by 
integration  of  the  velocity  components  in  the  space-fixed  axis  system, 
V-      ,    V      ,   and  W      ,  where 
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uS 

« 

V' 

= 

A 

• 

V 

w' 

•u/ 

(64) 


1L     f     V        ,  and  1^~   are  obtained  from  the  differential 
equations  or  initial  conditions,  and 


y'c- 


9  C 


ICO   +    /*"' 


dt 


'dt 


<tt 


(65) 
(66) 

(67) 


3.2.2     Unsprung  Masses 


3.2.2.1   Independent  Front/Solid  Axle  Rear  Suspension 


Wheel  Center  Location 


1 
X; 

t'c 

4 

— 

HL 

4 

zfc 

A 

A 

<i'z 

= 

% 

i'i 

llc 

X's 

^ 

A 

— 

V'c 

2'i 

*L 

+■ 


II   A|| 


HI  A 


+  \\A 


CL 
CL 

~(TF/2  +ATHFz) 
-b 

yi  -/fa 


RIGHT 
FRONT 


LEFT 
FRONT 


RIGHT 
REAR 


(68) 


(69) 


(70) 
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*'* 

i'c 

3i 

— 

a 

HMII- 

*4 

A 

2r+/>-(T(>/z)<I>2+J3 


LEFT 
REAR 


(71) 


where  AThFi   ,  and  ATt+Fz      are  the  front  half-track  changes  from  the  static 
values  of  Tp/±    interpolated  from  the  AThf      table  as  a  function  of  suspension 
deflection   o,   and   <^i    . 


Wheel  Orientation 


4>3  ~  (pf  =  ^ 


^/    ,    & 


t  ,    fz 


fi--f+  =  £fzs<Pe. 


,  from  differential  equations  of 
motion  or  initial  conditions; 

by  interpolation  of  tabular  inputs  of 
<pc     as  a  function  of  ^4    ; 

by  interpolation  of  tabular  inputs  of 

fF        as  a  function  of  time, 
or  from  differential  equation  of 
fF       if  curb  routine  is  active; 

(Rear  axle  roll  steer,  positive 
values  of   ^^   produce  roll 
understeer) 


(72) 
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3.2.2.2   Independent  Front  and  Rear  Suspension 


Wheel  Center  Location 


*; 

A 

ai 

— 

Jl 

2; 

i'c 

i-WA 


a 

TF/2+ATHF/ 


(73) 


4 

*'c 

9k 

- 

i)L 

Z'l 

i'c 

+  UA 


A 

K 

% 

«=* 

'i'c 

2'i 

*'c 

f-WA 


4 

i'c 

A 

— 

I'c 

% 

*c 

+  \\A 


a 


~(TF/U/STHFL) 


Zp+Jl 


-b 

Te/i  +  ATwz 

-b 

•fcfc+ATm) 


(74) 


(75) 


(76) 


where  A7#£3  and  -AT^z^  are  the  rear  half-track  changes  from  the  static 
value  of  Tfz./z.      interpolated  from  the  ffirifz      table  as  a  function  of 
suspension  deflection  £?>      and  <£f 
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Wheel   Orientation 

<j>1  ,   <?!.        -   interpolate  <pc    VS.     Si  }  <$z 

^3 ,    (j>q       -   interpolate  <Pcr   ITS.     J$  >  S+ 

fi  ,    yz        -   interpolate  ff    W>.  «■     or  from     ^> 


(77) 
(78) 


where  Kj$     ,      Kj$f    f    *~<£sz    ,     ICf$3    are   input  coefficients, 


3.2.3 


Beam  Axle  Front  and  Rear  Suspension 


.Wheel  Center  Location 


*', 

*c 

% 

— 

% 
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II  a| 

• 

z; 

ic 

4 

< 

ti 

- 

a 

Z'z 

i'c 

+■ 


All 


a. 

TF/2-/>F?F 

£F+f>f:t(TF/l)(pf  +<$, 
CL 

if  +/>F-(TF/l)<j>F+Jl 


(79) 


(80) 


*3 

Xc 

^ 

— 

i 

*3 

2L 

*  \\A 


Tz/2-pfa 
HR+p+(Td2)4e.  +4 


(81) 


A 

< 

ti 

~ 

y'c 

i'4 

Z'c 

+  IIA 


-b 


(82) 
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Wheel  Orientation: 


4i  >  4z- 

-      <Pf 

4z  >  4i 

=       <t>R. 

f,  >  T*- 

-        % 

fl   >V4 

Kg 

Hf       from  interpolation  or   T/^   equation 

Z  (83) 


3.3 


Vehicle  to  Ground  Orientation 


The  directional  cosines  of  a  line  perpendicular  to  the  ground  in 
the  local  terrain  region  under  wheel    u      are  calculated  by: 


cos  &62>. 

cos  A*; 

cos  ty 


=  Cos  <pc.  Sin  &a. 
-  cos  Oc.  cos  4>6. 


(84) 


where  &&l    , 
in  Section  3.4, 


'6-^   are  the  ground  slopes'  under  wheel   L       as  determined 


A  normal' to  the  wheel  plane  of  wheel    i,      is  determined  by  the 
directional  cosines: 


cos  ccyw. 

cos  /Syw. 

= 

A 

• 

cos  Xyur. 

•  • 

-  s  in  Yi 

COS  0C-     COS  f; 

5m0.  cos  ]//• 


(85) 
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By  subtracting  the  angle  between  the  two  above  normal  lines  from 
%/Z,      >    tne  camber  angle  relative  to  the  ground  is  obtained: 

^   =  %lz  ~arc  cos  \cos«y„c  cos*S3,  +  cos/3y„.  cos/fGZ,. 

-h  cosiryur.  co$*et^  C86) 

(Note  that  a  +  sign  for  (pc&l     corresponds  to  camber  thrust  in  *-i/  direction.) 

The  directional  components  of  a  line  perpendicular  to  both  the 
normals  to  the  wheel  plane  and  ground  plane  are  calculated  by  taking  the 
vector  cross-product  of  the  two  normals: 


D, 


<3  • 


cos  /?„/. 

cos  jfa. 
cos  irati. 

%     - 

cos  Jfa^ 

cos  r6x:_ 

cos*.yuk 

cos  fiyor. 

cosocGi,. 

C0S  fi&t 

Cos  or 


yurc 


COS  ot. 


GZ'C 


(87) 


Front  and  rear  wheel  steer  angles  in  the  plane  of  the  ground  are 
obtained  by  modifying  the  respective  steer  angles  relative  to  the  vehicle 
by  the  cosine  of  the  angle  between  the  assumed  axis  of  rotation  of  the  wheel 
(in  steer)  and  the  normal  to  the  ground  plane  at  the  individual  wheel,  where 


cos  odiUn: 

0 

COS/?j*rc 

= 

A\ 

• 

-  s'm  ft 

cos^. 

cos  fa 

(88) 


are  the  directional  cosines  of  the  assumed  axis  of  rotation  of  the  wheel 
(kingpin  axis  assumed  to  lie  in  wheel  plane) . 
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Thus: 
r/    =  fi(C0S  «^  COS  *m  +  COS  f^   Wfyuri    +***'**&,  COS  *}ur.)  (89) 

In  order  to  determine  appropriate  slip  angles  to  be  used  in  the 
calculation  of  tire  side  forces,  the  forward  and  lateral  velocities  of  the 
tire  contact  point  in  the  plane  of  the  ground  under  wheel  L      must  first  be 
calculated.  This  is  done  by  determining  velocity  components  along  the 
vehicle  -yL       and  u       axes  projected  onto  the  ground  plane  under  wheel 


The  velocities  of  the  tire  contact  point  along  the  vehicle  axes  are 
first  calculated  by  the  expressions  given  for  each  suspension  option. 

Independent  Front/Solid  Axle  Rear  Suspension 

Forward  velocities  of  hubs  (along  X  -axis) 


U,   =  U.  -  (TF/2.  +ATm)Z   +  dF+-S,)a 


Lateral  velocities  of  contact  points  (along  it     -axis) 


Tt-     V  +  al~(2F^P-(klcost^(P+^Jd:+     ^t 


(90) 


74=  V-bR.-  (Ze+f3)P-(p+TRfa/2i-h3cosVte)(4a+p)  C91) 

Up  v-bZ-  teetS^P-Cp-Tth/l+h+c&YwXbi+P) 
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Vertical  velocities  of  contact  points  (along  £  -axis) 


W^K-aA  +(%/l+ATHF,)P  +i,  i-k/axfylP+diji) 


dS, 


w2^  im~<xQ  -  (tf/2+-ATHf2)p  +j2  +  kL  cceVte(P+  ~-  sz)  (92) 


w- 


=  i^-tbQ  f  4  -  tp4>z-\/2  -hzax^dz+p) 


Independent  Front  and  Rear  Suspension 
•Forward  Velocity  of  Hubs: 

u.,  -   u    -  (TF/2+ATHFl)Z*-  (zF+Sl)4 

dz   =  a  +  lTF/1  +  t±THfl)tZ  ±UF+sz)Q 

Lateral  Velocity  of  Contact  Points: 

v,  -  v- +all-(zf+sl)P-(klcx>t>u)(P+&2,)  +  dAJmt 

ddi  J  S,       ' 

V2  =  if  +  *£-(zFt£)p  -(h1co5Xn2XP-£rF  *)  "    #**A 
2T3  =  IT  -  lo£~  U^65W-(k3CtZth3)(P+^  S3)  +    ^fe  i3  (94) 


1M  --  nr  -  bll-te^P-Ch+atrtoXPt  #£)  -  ^*  4 


fj<^>>       J/47 
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Vertical  Velocity  of  Occupant  Points 


vrt  ^nAr-aa  +  (TF/Z^Tm)P  +$,  *  h,  cosfa,  (Pi-  ^  S,) 
>r2=  Iat-ciQ  -  (TF/2tATHFZ)P  +S2  i- h.zcos^(P*   ^  Si) 
2ua  =  70-^4*  (Tt>l2+hTm)P+S3  +h}CDSfa(f>+   ^f  S9) 
fy  ^  iff +b<l-(re/2+ATm)Pt<b+ l^te fa(P  +  ^   £f) 

Solid  Front  and  Rear  Axle 


Forward  Hub  Velocity: 

U,  -  U  -  (TtJ2-f>F4>F)  +  (zFtJl  +  TF4*F/2)a 
V^  -  %L  +  (TF/Z+/>F<k)+feF+S,  -TFfa/2)  a 

u^u-  tyi-p  $>  +  &t  +s3 + Tt$/2 )  a 

Lateral  Contact   Point  Velocity: 
Vf  =  7T  +  a£-  (ZF+J,)  P-(pF+  Tpfc/Z  +k,  CDSi^(P^4f) 

v1=v+ol(Z~  (zf+5,)?  -  (pF-TF4>F/<2+k1  a>$^)(P+4>F) 
Vb  ^ir-bH-  fei+sjp  -  (/>+%cl$L  +  h3&5  XwVP+i) 
Vj  =1r  -bR.  -  (zR+<%)P  -  (p-T^/2  +hfa&jftof)(Pr&i) 


(95) 


(96) 


(97) 
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Vertical  Contact  Point  Velocity: 

%)z-  W  +J,  -aGt-  (pF<pF  *-TF/Z-h%CDSfih2)(pi-$F) 
m,  =  wt/3  +bQ.  ~  (p4n  -T&/1-  k^cosfaJCPtfa) 
Wa  -   W  +  S3  +  bQ.   -  (p<j>z  +\I2 -  h4-0SPM)(Ptfa) 


(98) 


The  angle  between  the  vehicle   X  -axis  and  its  projection  onto 
the  ground  plane  is  then  determined.   By  taking  the  vector  cross-product  of 
the  U       -axis  and  the  normal  to  the  ground,  the  directional  components  of 
the  X-      -axis  projection  in  the  plane  of  the  ground  are  obtained. 


**i 


**i 


C*i      = 


CaS  *4Mi 

COS   ty 

cos  fa*. 

COS    CCy 

C0SctGg, 


COS 

* 

COS 

***} 

COS 

"y 

cos 

<*£*;• 

COS 

4 

COS 

^te'i 

(99) 


where 


COS         Oty 

Cos  jSy 
COS     *y 


0 

1 

0 


(100) 


are  the  directional  cosines  of  the 


V 


-axis. 
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The  angle  between  the  ~£     -axis  and  its  projection  onto  the 
ground  plane  is  then  determined  by  the  vector  dot  product: 

~  (Cosoex)a.z     +  (cos  #*)  bz.   y  (cos  *x)  c*. 

COS   &x&i     =   ■ j     z (101) 


and, 


*f*  ft*;  =    s9* 


COS  £Y  - 


'Xi 


y*i£.  /  />%.  7 c*. 


where 


cos  ocx 
Cos  /S^ 
Cos  Xx 


/ 

A 

• 

0 
0 

(102) 


are  the  directional  cosines  of  the  %     -axis. 

The  forward  velocity  of  the  contact  point  in  the  ground  plane  is 
calculated  as: 


H;     =  *«  C0S   eXG£  -   **  &*  Ox&£ 


(103) 


Similarly,  the  directional  components  of  the  V      -axis  projection 
onto  the  ground  plane  are: 


*Vi     = 

cos  2rx 
cos  >&i,. 

cos  /SK 

COS  /fai,. 

hi  - 

COS     Xy, 

COS  oCfia.>. 

COS  fy 

cos  *6g,£ 

<wt    ■ 

cos  /8-t 

COS  jSG# 

COS     Otyt 

COS  oL^jit. 

c 

(104) 
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And  the  angle  between  the  -  M      -axis  and  its  projection  onto  the 
ground  plane  is: 


and 


cos  <tvs.  =  ("*  1*1  aJi  *  (C0SSy)  gg  *  fafAlfZL 


(105) 


plane  is 


s$n  0vct.     =  •*?<* 


cos  }y  - 


e*t 


Finally,  the  lateral  velocity  of  the  contact  point  in  the  ground 


3.4 


Tire  Contact  with  the  Ground 


(106) 


3.4.1 


Ground  Contact  Point 


Determination  of  the  tire  "ground  contact  point"  is  accomplished 
by  passing  planes  through  the  wheel  centers  (  p'^     >     li'i      >     t1 '       ) 
perpendicular  to  both  the  wheel  and  ground  planes  at  the  individual  wheels. 
At  each  wheel,  the  point  that  lies  in  all  three  planes  is  designated  the 
"ground  contact  point".  The  characteristics  of  the  ground  under  wheel 
are  determined  in  Sections  3.4.2  through  3.4.5. 


Thus, 


KH    =    %'i  cos  oLy^.    ¥■   y'.  cos  fy^   +   Vi  cas  ^yur,:     (1 


07) 


% 


%i  cos  *0i<.     /   ViCosfi&z<.    ■+■  3%  cos  *9Z  ,.      (108) 

(109) 
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where  Equation  (107)  is  the  equation  for  the  plane  of  wheel  i      ,    (108)  is 
the  equation  of  the  ground  plane  through  the  terrain  point  directly  under  the 
wheel  center,  (  %'L    ,   v£  ,  2;'L      ),  and  (109)  is  the  equation  of  the  plane 
through  the  wheel  center  perpendicular  to  both  the  wheel  plane  and  ground 
plane. 


Let 


4 


COS  <£ 
COS   Of, 


w 


cos  a 


yt4r- 


COS  * 


yt*c 


<?2V 


cos  a 


oz'c 


cos  }, 


tf?V 


D>< 


DH 


(110) 


Then  a  point  that  lies  in  all  three  planes  is  determined  from 


Ik  II 


1 

x&pt 

V 

y'&Pi 

= 

** 

$GPi 

*si 

(111) 


The  radial  load  on  the  tire  is  subsequently  found  by  determining 
the  deflection  of  the  tire  at  the  ground  contact  point. 


A<-  -  *  ph^y  *  (yi-*„£y  +  (*i  -  *ij 


(112) 


where  A  i        is  the  distance  between  the  wheel  center  and  ground  contact 
point;  and 


hi 


A  - ,  for  A;  <  R„ 
RW)     for  f?„±Li 

where  h^i      is  the  rolling  radius  of  the  tire. 


(113) 
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Finally,   the  radial  tire  force  is 


O,    /br  (/?„-/z<)=0 

*r  (**>  -A J  ,  for    0<  (z„  -fa)  <  <rr 


(114) 


The  directional  cosines  of  the  line  of  action  of  the  tire  radial 
force  for  wheel  I         are: 


COS  dR.   - 


A; 


cos  A.   -  *»<- 


cos  X*.   s  igSlllL 

A, 

The  corresponding  direction  cosines  in  the  vehicle  coordinate 
system  are  obtained  from 


(115) 


COS  oLl.   I 
1 

c°s  t^ 

= 

AT 

• 

cosfAi   1 

COS  ote£ 
COS  fa. 

Cos  *Ki 


(116) 


3.4.2 


Flat  Terrain  Ground  Characteristics 


The  default  terrain  is  defined  as  being  flat  and  level  at  an 
elevation  of  0,0.   Therefore,  the  ground  characteristics  under  wheel 
are  defined  as: 


2-&L 

= 

0 

0 

<P&i, 

= 

0 

0 

&6L 

= 

0 

0 
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3.4.3     Variable  Terrain  Ground  Characteristics 

Varying  terrain  is  defined  by  the  use  of  input  terrain  tables. 
Therefore,  it  is  necessary  to  determine  the  characteristics  of  the  terrain 
under  wheel  _  l-        as  defined  by  Zs-i  ,     4>sri     ,     &&L       by  interpolation 
of  terrain  table  data*. 

The  calculational  procedure  for  determining  the  elevation  and 
slopes  of  the  terrain  under  each  wheel  of  the  vehicle  from  tabular  input 
data  is  described  in  detail  below. 

Step  1 

The  highest  numbered  terrain  table  (1  through  5)  applicable  to  the 
wheel  is  determined  by  sequentially  testing  if  the  wheel  is  located  within 
the  X'  and  Y'  bounds  of  each  table.  Note  that  in  regions  of  overlapping 
tables,  the  program  will  use  the  terrain  as  defined  in  the  highest  numbered 
table.  Therefore,  if  it  is  desired  to  use  a  small  or  "fine  grid"  table  to 
describe  in  more  detail  terrain  that  is  also  within  the  bounds  of  another 
"coarse  grid"  table,  the  number  of  the  "fine  grid"  table  must  be  larger  than 
that  of  the  other  table. 


Step  2 


The  particular  grid  segment  within  which  the  wheel  is  located  is 
determined  and  the  corner  points  labeled  as  shown  in  Figure  3.6. 

Step  3 

Determine  if  an  interpolation  boundary  cuts  through  the  segment. 
For  Y  boundaries,  a  boundary  cuts  through  the  segment  if 
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NOTE:    1,2,  3,  4  ARE 

CORNER  POINTS  OF 
TERRAIN  SEGMENT 
IN  WHICH  WHEEL  i 
S  LOCATED 


POSITION  OF  WHEEL  i 


Figure    3.6     TERRAIN  TABLE  GRID 
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where: 


Yi  <  YdD&y  <  Y3 

(Note:   Subscripts  1,  2,  3,  4  refer  to 
corner  points  of  grid  segment.) 

An  angled  boundary  cuts  through  the  segment  if 

Mi  <  Xlc£PT  -  X<U     ,   or 
^lcbpt-^61   and  XgcEFT  ^   *64      ,  or 

/tfspT  -  XaKY  +  (  Y&i  -  Yd)  cln.  'faoRr  (117) 

Yzcepr  =  Xbpxy  +  (YU  '  Ye)  ctn  Yep/ty 


Step  4 


,/ 


If  no  boundaries  cut  through  the  segment,  caluclate  Z&L       by  two- 
way  interpolation  using  the  elevations  of  the  corner  points  of  the  segment 

H&f        ,    Z&z.      ,  £&*,       and  £  &4    . 


Then 


*.-«^"-'f^^j  (120) 
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<?w= 


ardaJ^i^hh    sec  9^        .fcrrf-ySf  (122) 


where  £<w,2   and  z-^^-f   are  "intermediate"  terrain  elevations  at  points 
on  the  Y*  sides  of  the  terrain  segment  at  X-  I      as  determined  in  the  first 
of  the  two-way  interpolations. 

Step  5 

If  there  is  only  an  angled  boundary  cutting  through  the  segment,  a 
determination  is  made  as  to  whether  or  not  the  wheel  has  crossed  the 
boundary.  The  wheel  has  crossed  the  boundary  if 

i&W+  (  y'L  ~  Y3)  dnfaPZY    <    x'L  (123) 

(a)    If  the  wheel  has  not  crossed  the  boundary, 

(1)  starting  with  corner  point  2,   the  first  grid 
point  along  the  Y'   »  Y'p2  ec*Se  °f  tne  segment 
for  which  X'<XLCEPT  is  found  and  designated 

as  point  2.  Corner  point  1  is  then  relabeled 
accordingly  (one  AX  increment  toward  X  ) . 

(2)  Starting  with  corner  point  4,  the  first  grid 
point  along  the  Y*  _,  Y'   edge  of  the  segment 
for  which  X'<XRCEPT  is  found  and  that  point 

is  designated  point  4.  Corner  point  3  is 
relabeled  accordingly  (one  AX  increment 
toward  X  ) . 

D 
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(3)  If  X'        =   X'      ,      i.&i     ,        Q&L_    and       <P&u     are 

computed  as   in  Step  4. 


If  X?G2   =   X'G3 


,   =  arttaj*'6*'  *'*+  ) 


96l  =  arctcm.(  -^~~  i 


Zii  =   *«  +  fy-  -Y^^&m^  -(?i-^2)tan.9&L       (126) 


I£  X'ci   =   X'G4' 


9a  ^arcUn^61^62'  ) 


(i'&4  -  ^4? ) 

AY 


^  =  ardan.  [       *w  -  Sec0s-j  (128) 


£<srt   =  Equation  (126). 

If  X'    <  X'   ,  designate  the  grid  point  one 
A  X  increment  toward  X  from  point  3  as 

D 

point  5.   Then 

&6-L  =  Equation  (124) 

<f>6i  -  archxn  [^L^sec  ft]  (129) 

z&l  =  Equation  (126) 


If  X'   >  X*   ,  designate  the  grid  point  one 
GI      G4 

AX  increment  toward  X  from  point  1  as 

D 

point  5.  Then 
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9&-L  =   Equation    (127) 

fe  =■  ordan[^k^}  sec 9^  (130) 

2^.  •  =  Equation  (126) 

(b)    If  the  wheel  has  crossed  the  boundary, 

(1)  Starting  with  corner  point  1,  the  first  grid 
point  along  the  Y'  , ,  Y'   edge  of  the  segment 
for  which  X'>XLCEPT  is  found  and  designated 

as  point  1.  Corner  point  2  is  labeled  accordingly 

(one  AX  increment  toward  X„) . 

b 

(2)  Starting  with  corner  point  3,  the  first  grid 
point  along  the  Y'  _,  Y'-,.  edge  of  the  segment 
for  which  X'>XRCEPT  is  found  and  designated 

as  point  3.  Corner  point  4  is  labeled  accordingly 

(one  AX  increment  toward  X„) . 

b 

(3)  If  X'G1  =  X'G3,   H^6  ,   Q&L      and  (p&i   are 
computed  as  in  step  4. 

If  X'g2  =  X'g3' 

&Cr0  =  Equation  (127) 

<P&L  =  Equation  (125) 
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If  X'G1  A  G4» 


^fi-i.  =  Equation  (124) 

A)    •  =  Equation  (128) 

•z'  .  =  Equation  (131) 


If  X'  «  <  X'r  *,  designate  the  grid  point  one 
AX  increment  toward  X  from  point  2  as 
point  5.  Then 

q    ,  =   Equation  (127) 

4>ti  =  arc-tan  (  (?kl*f^  sec 9^  )  C132) 

-I'    .  =  Equation  (131) 

If  X'  .  >  X'   ,  designate  the  grid  point  one 
AX  increment  toward  X  from  point  4  as 
point  5.   Then 

q    .  =  Equation  (124) 

4>&L  =  arctan[  ^'J^  sec 9&t)  '133' 

£'f  =     Equation    (131) 
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Step  6 

If  a  Y'  boundary  cuts  through  the  segment  and 


(a)  if  Y1.  1  YDDRY'  tne  §rid  segment  one  AY  increment 
toward  Y  is  used  with  the  corner  points  designated  as 
shown  in  Figure  3.6.   Proceed  as  in  Steps  3  through  5. 

(b)  If  Y'.  >  Y    ,  the  grid  segment  one  AY  increment 

toward  Y  is  used  with  the  corner  points  designated  as 
b 

shown  in  Figure  3.6.   Proceed  as  in  Steps  3  through  5 
using  the  following  expression  for  the  roll  slope,  S^j^    , 
in  Step  4. 

**  ~~ arcian  ( %^fr sec  &^)    ' w  *  <  Y« 

When  boundaries  cut  through  the  gird  segment  in  which  a  wheel  is 
located,  the  elevation  and  slopes  of  the  terrain  under  the  wheel  are  computed 
from  data  supplied  for  nearby  terrain  segments.   To  minimize  errors  that 
might  result  from  use  of  grid  points  too  far  from  the  actual  terrain  segment 
in  which  the  wheel  is  located,  it  is  recommended  that  the  angle,  ifopfZY    *. 
be  restricted  to  the  following  range: 

(ft-  arctan  l£  )    2>  fMy  >    a titan  £L 

A  minimum  of  two  tabular  values  between  like  boundaries  (i.e.,  two 
angled  boundaries  or  two  Y"  boundaries)  or  between  a  boundary  and  the 
beginning  or  end  of  a  terrain  table  must  also  be  provided  for  proper 
operation  of  the  subroutine. 
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3.4.4     Curb  Impacts 

The  use  of  the  curb  impact  subroutine  is  determined  by  a  series 
of  logical  tests.   If  no  curb  impact  input  data  are  supplied,  then  this 
subroutine  is  not  activated  during  the  course  of  the  computations.  However, 
if  input  data  are  supplied,  then  each  wheel  center  is  tested  for  proximity 
to  the  curb  immediately  after  the  vehicle  position  and  orientation  subroutine 
has  been  completed. 

if  Rw  -  (yit  -y'L)  »  then  for  wheel  l 

1.  use  tabular  input  of  ijfp       until  set  to  differential 
equation  of  Iff       by  any  one  wheel, 

2.  bypass  terrain  input  and  set  «£   =  &&      =  y&     =  0, 

3.  retain  coefficient  of  friction  =  .AMU^ 

If  /<W  -  (v'i  -  Vcl)     y   where  Vcu         is  the  location  of  the  start 
of  the  last  curb  slope  then  for  wheel  ^ 

1.  continue  to  use  differential  equation  for  determination 
of  fF       , 

2.  bypass  terrain  input  and  set 

3.  reset  coefficient  of  friction  to  the  nominal  ground 
value 
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4.  if  all  wheels  pass  this  test,  reset  calculation  time 
increment  to  At, 

5.  return  to  vehicle  to  ground  orientation  subroutine. 

If  ty  1 1-  y'ci.  )-  R.w-  (yif  -y't)  then  for  wheel  'l       >   enter  into 
curb  impact  subroutine  calculations, 

1.     use  differential  equation  for  determination  of  fF 
with  most  recent  tabular  value  as  initial  condition, 


set  calculation  time  increment  to  At  , 

c 


3.     modify  the  tire-ground  friction  coefficient  by  the 

curb  friction  multiplier,  i.e.,  -^i  =  (AMUi)  (AMUC) 

Let  the  matrix  \\Aj\\  be  defined  as  follows: 


A- 


COS^u  COS&j 

-$LafL 

GOSfi  SirtQj 

cos4>l  Sirifi  cos  Gj+sLntyi  sinGj 

C05<pLcosfL 

C05<j>L  SLniisLnty  -  Sinji  C059j 

Sin 4>l  5Ln>fi cosfy  -  cos <f>i  si-n,  Qj 

sinfacosfc 

e'u^(j>i  stiifi  5ut  Qj  *■  co5(pL  cosOj 

(134) 


where 

9Js4(j)  >    j='Ur15,->0;1,-     •25,26 

(Note  that  the  above  matrix  corresponds  to  the  sequence  <pj     ,  fj       ,  % 
which  will  apply  the  camber  angle  <j>L        as  a  rotation  about  an  axis  parallel 
to  the  sprung  mass   /'    axis.) 

Let  the  matrix   HBI|   be  defined  as  IIAil'llAjll    ,  and  #23  ,  ^33 
represent  elements  of  l|BI|    with  the  usual  notation. 
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Solve  each  of  the  following  expressions   for    h-j 


1. 


2. 


5. 


7, 


«i     -     ~ 


i'c 


h!<    - 


hi:    ■=- 


h!i  = 


kV  = 


fcV  = 


k:    = 


5a3-  023  -fan  fez 

&2$  -  #23  to  ^c3 


M  y'a^y'^yc 


C3 


,for  y'c^y'itfjU 


&r  $<+±  alty'c 


^b 


jfor   $'o5  ^  &  <  y'c* 


•for  y^^y'c 


(135) 
(136) 
(137) 
(138) 
(139) 
(140) 
(141) 


Using  each  of  the  values  of      h_$       solve  the  following  expression 
for  three   sets  of   (     %'■     ,      y'j      ,    %>:      ) , 


*;• 

0 

y'j 

= 

y'i 

•h 

3 

e 

0 

tj 

h 

j 

(142) 


Using  the  solutions  of  (  /■'/    ,  y'j      ,    Z.'j  )  for  each      , 

test  the  value  of  h.^     for  the  range  corresponding  to  its  y j       solution 

shown  above.   Retain  only  k./     ,     Y-'j       ,  y'j       ,    i.  /  for  which       is 
in  the  correct  range  for  its   It/    soluti( 


.on, 
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Using  the  retained  values  above: 

If  k.j'^£iA/c     ,  then  this  individual  radial  spring  is  not  in 
contact,  index  J       by  +1  if  j       >  26  proceed  to  calculation  of  radial 
tire  forces. 

If   hJj    <   £\a/l       ,  calculate  and  store  the  following: 

x'-  -  %'. 

COS  #;   =       c    -    J 

y'  -  y' 
*j 

COS  *;     =  iLl£/  (143) 


Index  j       by  +1  or,  if  .  j      >   26,  proceed  to  calculation  of  radial  tire 
forces. 

Calculation  of  radial  tire  forces: 

For  each  rtj    ,  look  up   ry  in  input-generated  table  of  Ff 
versus  (£wL-hj),   calculate: 

1        £*S    J  '  (144) 

J— 24 


*U 


J '-26 


(146) 


n 


and 


%  -  H(L'**iY  +  (LF*9tf  +  (Lr*fiy 


(147) 


Therefore,  the  directional  cosines  of  the  line  of  action  of  the  tire 
radial  force  are: 


COS  off. 


Cos  fa.     -  - 


Cos  **,  = 


Z-r  F**', 


Fm 


Lf* 


*y'i 


/> 


LF*£ 


(148) 


An  "equivalent"  ground  contact  point  and  "equivalent"  ground  slopes 
are  now  determined  for  calculation  and  application  of  side  and  circumferential 
tire  forces. 


The  equivalent  rolling  radius  of  the  tire  is 


*t~ 


w 


n  Get    .,   (*r  ~  0  °t 


,  for  Ai  <  (R„  -  a-7 


(149) 


and' the  corresponding   V   and   2    coordinates  of  the  equivalent  ground 
contact  point  are: 


y&r£  *  y-  r  ^  **-fac- 

I'eri   B     $'c    +  **;    cos  £,. 


(150) 
(151) 
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"Equivalent"  terrain  slopes  for  this  case  of  nonplanar  terrain 
contact  are  established  by  the  use,  as  a  normal  to  the  equivalent  terrain, 
of  the  projection  of  the  normal  to  the  actual  terrain,  at  the  equivalent 
ground  contact  point,  in  the  plane  determined  by  the  radial  force  vector 
and  a  normal  to  the  wheel  plane.   In  this  manner,  the  line  of  intersection 
between  the  equivalent  ground  plane  and  the  wheel  plane  is  made  to  remain 
perpendicular  to  the  radial  tire  force  vector. 

Determination  of  the  slopes  of  the  equivalent  ground  plane  is 
accomplished  by  means  of  application  of  the  following  relationships: 


q>&i    =  arcsun 


-  Dm 


^DNfL^W2zL^DN32i  )  (152) 


0*1  =  *****    ( ^[|   )  (153) 


'L 


where 


DHk  -Oik  $Ui<PU  -*&  cos  fa 
DH2l  =  -bid  CDs4>6i'(ai-hcl  )SLn4>^L 
VN\  -  (<£+%) cos 4'6i  ±bLcL  sin<p'6i 


(154) 


<f>6-L       and  @GrL     are  the  equivalent  ground  slopes,  <J &l    is  the 
actual  ground  slope  at  the  equivalent  tire  contact  point  (note  that  Qg-i   =  0 
for  the  case  of  curb  simulation),  and 
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h  ■ 


cos  /svw,. 

cos  **. 
cos  J>yur. 

COS  Oijp. 

COS  qiyur£ 


COS  ?£. 

cos  *yu,. 


'COS  c6~. 

cos  ocyur. 

COS  /3yurc 


(155) 


are  the  directional  components  of  a  line  perpendicular  to  both  the  normal 
to  the  wheel  plane  and   Pg^ 


3.4.5 


Road  Roughness 


The  road  roughness  algorithm  is  based  on  the  assumption  that 
roughness  data  is  in  the  form  of  small  elevation  changes  from  the  ground 
plane  and  is  supplied  with  a  constant  increment  between  elevation  points. 
The  roughness  data  is  assumed  to  vary  in  the  X'  direction  and  is  constant 
in  the  Y'  direction.   Tire  enveloping  .effects  are  approximated  through  the 
use  of  the  distributed  radial  spring  tire  model  in  a  manner  similar  to  that 
used  in  the  curb  impact  model. 


Let 


A„   be  the  distance  increment  between  roughness 

b 

elevation  points 

Z'   be  the  mth  elevation 
Gm 

The  overall  solution  procedure  involves  determining  the  intersection 
between  radial  spring  j     and  the  straight  line  connecting  elevation  points 

yYi.    and   fit  +1.   If  the  point  of  intersection  lies  between  the  limits 
defined  by  the  two  elevation  points,  the  solution  is  legitimate  and  the 
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radial  spring  deflection  and  force  are  computed.   If  the  intersection  is 
outside  those  limits  r/L     is  incremented  and  the  procedure  repeated  until  an 
intersection  is  found.   Note  that  the  terrain  outside  of  the  limits  of  the 
data  is  assumed  to  be  flat  and  level  at  an  elevation  of  i,1^    =   0.0. 

The  transformation  matrix  between  the  radial  spring  J      and 
space,  II 5 II        ,  is  computed  as      ll 311-    HAII'llAjll         where  //A/// 
is  given  by  equation  (134). 

The  intersection  between  radial  spring  j      and  the  straight  line 
between  elevation  points  Wl      and  Ht       +1  is: 


/ 


*/    = 


4        3>t>  f  £'&fi/y-i  ~~^6vn 

1       B33{  && 


XfLt 


$3$ 


t  I  _    -£/ 


■yvr\  -    A6M 


A& 


(156) 


If  tfrft-  Y-J^KfrMfl     >  tne  intersection  occurs  within  the  limits 
for  this  terrain  segment  and  the  distance  between  the  wheel  center  and 
intersection  point  is  computed: 


HI,-- 


X}-    Xi 


(157) 


6n 


If   ^j  >  Kw l        >    the  spring  is  not  deflected,  and  calculations 
proceed  for  the  next  radial  spring. 


If  If  a '  <  Rtfl  ,  the  direction  cosines  for  this  spring  are 


computed: 


COS  o(j    - 


CDSfij    = 


COS  %j   - 


X  i      X-j' 


*;■ 

y; 

•-# 

/t'y 

zi 

-*} 

(158) 


A.V 
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and  the  radial  spring  force,  f  a      ,  is  obtained  through  interpolation. 

When  calculations  have  been  completed  for  all  radial  springs,  the 
resultant  force  components  are  computed: 


XFzz!C    -     ZFj&sXj  (161) 


(159) 
(160) 


and 


hi     -    J&fkx'i)*  *■  C2t*w)*  +  (SFeti)2  (162) 

Therefore,  the  directional  cosines  of  the  line  of  action  of  the 
tire  radial  force  are: 


CDi'Wjii    *  - 

ZFfr'L 

Fki 

cpsfa  =  - 

-2  nzy'L 
Fee 

CDSXfr   =   - 

Z'fke't 

(163) 

An  "equivalent"  ground  contact  point  and  "equivalent"  ground  slopes 
are  ,now  determined  for  calculation  and  application  of  side  and  circumferential 
tire  forces. 
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The  equivalent  rolling  radius  of  the  tire  is 


£ 


wL 


flL=      * 


£ 


VJl 


-~     ,    for  (eWL-  afi )  <  kL 


Fee  (hrL-l)ofL 


V^-Tl 


*-TL 


,-for  hi^C^Wi-cm-i) 


(164) 


and  the  corresponding  coordinates  of  the  equivalent  ground  contact  point 
are: 


£.'c-  +  hi  cos  X& 


(165) 
(166) 


"Equivalent"  terrain  slopes  for  this  case  of  nonplanar  terrain 
contact  are  established  by  the  use,  as  a  normal  to  the  equivalent  terrain,  of 
the  projection  of  the  normal  to  the  actual  terrain,  at  the  equivalent  ground 
contact  point,  in  the  plane  determined  by  the  radial  force  vector  and  a 
normal  to  the  wheel  plane.   In  this  manner,  the  line  of  intersection  between 
the  equivalent  ground  plane  and  the  wheel  plane  is  made  to  remain  perpen- 
dicular to  the  radial  tire  force  vector. 


Determination  of  the  slopes  of  the  equivalent  ground  plane  is 
accomplished  by  means  of  application  of  the  following  relationships: 


<4r/  =  orcein. 


-  VI&L 


6&i  =  arciam 


\dU1l+DN2l+DM}i 


I  DM3t 


(167) 
(168) 
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where 


DM1 ■    =     (Cl±bl)i>in9k   '  <XiCLC0S0'&L 


(169) 


ffij.^;   and  Q&i         are  the  equivalent  ground  slopes,  0^    is 
the  actual  ground  slope  at.  the  equivalent  tire  contact  point  (note  <p&i    =   0 
for  the  case  of  road  roughness),  and 


a;   = 


bC 


CL    = 


COS  t&i  CjDZ  oift. 

cos  'iywi       co<>  otyn  j 


(170) 


are  the  directional  components  of  a  line  perpendicular  to  both  the  normal  to 
the  wheel  plane  and  F&l 


3.5 


Suspension  Forces 


Suspension  forces  include  effects  of  viscous  and  coulomb  damping, 
suspension  stops,  auxiliary  roll  stiffness,  anti-pitch  linkages  and  jacking 
forces,  and  are  calculated  based  on  suspension  displacement  and  velocity. 
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3.5.1     Front  Suspension 

The  displacement  and  velocity  for  front  suspension  force 
calculation  are  dependent  on  the  suspension  type.   For  independent  front 
suspension  they  are: 

Si  =  sL 

5L    =    Si 

For  a  beam  axle  suspension,  they  are: 


Si.   =    s,  -c-tf  5^/2 
k     =  4  '(-V1  T*F4>F/2, 

The  front  suspension  forces  are  given  by: 


where 


=-  M5  Q     is  the  static  sprung  mass  weight  at  wheel  ^ 

9 

Cp  S[  is  viscous  damping  force  at  wheel  I 

hfi  is  the  coulomb  damping  force  at  wheel  t- 

given  by 


fc;   =    ( 


HF 


CfS^(Si)       ,  for  \SL\>GF 


(171) 


(172) 


%  -  ^  "s $  -  *4  -  t<n  ~ ?zn  -Gk,-F*n  +  fa  (173) 


(174) 
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is  the  combination  of  suspension  force  due  to 
displacement  from  equilibrium  and  suspension 
stop  force  as  illustrated  in  Figure  3.7  and 
given  by: 


F2FL    =    1 


C$i 


F  ~Jl 


,    &r  HF£>SL>SlFC 


Kf  \  +  KF£  (Sl~SIF£)  f  ^  ££  -iVf ) 3       ,-Pcr    SL  >SlFB  and  sgnfi  -  5gn.S^ 


(175) 


f>£t      is  the  suspension  force  at  wheel  I      due  to 
auxiliary  roll  stiffness  given  by 

f    .  g  -(-tf\  z     '  1  for  ^dependent  front 

'       suspension 


^     =      -M/ 


gp4 


for  beam  axle  front 
suspension 


.176) 


(177) 


^APl 


is  the  anti-pitch  force  at  wheel  l      given  by: 


f*L^  -m\ 


FcLht 
1t 


CoS(j>i  cos% 


(178) 


where  (APpJ^r       is  the  anti-pitch  coefficient  inter- 
polated from  the  front  anti-pitch  table  as  a  function 
of  suspension  deflection 


83 


f, 


JF  I 


is  an  approximation  of  the  jacking  force  induced 
on  an  independent  suspension  at  wheel  ^  given 
by: 


J4i 


SATafc 


fm  .-OtoOtt-fcaQjhi  g  -C-lff^  JJ~     ("9) 


Note  that  -Sfp-i     -   0  for  beam  axle  suspensi 


on, 


^7  i   <£  1  Sj  > 


^ 


2  »  u2  »  J3?  ^ 


Figure  3  7    GENERAL  FORM  OF  SIMULATED  SUSPENSION 
LOAD-DEFLECTION  CHARACTERISTICS 
(REFERENCE  2) 


3.5.2     Rear  Suspension 

The  displacement  and  velocity  for  the  rear  suspension  are  similar 
to  those  of  the  front.  For  an  independent  rear  suspension  they  are: 


4-  =  k 
$  -  k 


The  rear  suspension  forces  are  given  by: 
(L 

via 


where 


C 


*  9 


n(iL)     >  for  l^  ~G* 


(180) 


And  for  a  beam  axle  suspension: 

*i  =  is  -(-if  %  fa /I 

5  =   £-H/S&/2  (182^ 


—  /4sgr  -  ^ J-  -/^  -^  -^.  -^  f  /^.  (183) 


r-; — 7v  /VI5  J  is  the  static  sprung  mass  weight  at  rear 
wheel,  L 

m 

Cp  Si  is  the  viscous  damping  force  at  rear  wheel  l 

t'mi  is  the  coulomb  friction  force  at  wheel  t- 

given  by: 


8S 


hz-L 


is  the  combination  of  suspension  force  due  to 
displacement  from  equilibrium  and  suspension 
stop  force  as  illustrated  in  Figure  3.7  and 
given  by: 


kti     =  ■< 


^Rl  $L 


/■for   S1RE>  SL>SZ1ZC 


&  S,  +  <ac"ft  -Size)  +  Kk  CSl  -Slzcf     ,  for    £ < &&  ***  $3*  4  =  *3«  & 
hSi+  fa(Sc~SL&)+Xie(fi-Jlfi£p    ,for   SL  >Sl^  and  #*.-£  =  sgn  % 


(185) 


{-■AR.L  is  the  suspension  force  at  wheel  i       due  to 

auxiliary  roll  stiffness  given  by 


FARi  =  (-»L 
Fag  =  -M 


for  independent  front 
suspension  (186) 

for  beam  axle  front 
suspension  (187) 


'APl      is  the  anti-pitch  force  at  wheel  <-       given  by: 

Fa  kL 


FAo;    =    -fAfc) 


-L±k. 


"Af>6 


'4         /2 


coS<fo  cos l^i 


(188) 


where   {APjzJfy       is   the  anti-pitch  coefficient   inter- 
polated from  the  front   anti-pitch  table  as  a  function 
of  suspension  deflection 
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f- 


7f  L 


hH 


is  an  approximation  of  the  jacking  force  induced 
on  an  independent  suspension  at  wheel   ^  -  given 
by: 

Note  the   fjpi      ~   0  f°r  beam  axle  suspension. 


Note  that,  the  Anti-Pitch  coefficient,  APf       or   A/^   ,  is  the 
ratio  of  the  positive  or  negative  jacking  force  (assumed  positive  for  anti- 
pitch  effects  at  the  front  and  rear  suspensions  for  forward  braking)  acting 
on  the  wheel  to  the  corresponding  moment  of  the  applied  circumferential 
force  about  a  line  parallel  to  the  vehicle  M        axis.  The  units  of  the  anti- 
pitch  coefficient  are  lbs/lb-ft.  Tabular  entries  of  the  coefficients 
correspond  to  given  displacements  of  the  wheel  centers  from  their  design 
positions. 


Front  and  rear  suspension  jacking  force, 
by  the  following  derivation. 


Fjp  i       ,   approximated 


The  -u       and   g.   coordinates  of  the  ground  contact  point  of 
front  wheel  L       may  be  expressed  as 


(190) 


where 


COSfc  6LYiQt 


COSfc  SLn&i 

cos  4>i  S'vfi^fi  slnOi  -  Sirups  co$9L 

sin<f>i  siriljfi  sinOi  *■  cos^  cqsQl 


(191) 


37 


where 


Equating  the  work  done  by  the  tire  forces,  Fyu.L  and  fzu  i  , 
to  that  done  by  the  suspension  force,  Si  ,  in  a  virtual  displacement  of 
the  suspension  A  S'L        yields: 

c     _  r    9y  ,    r  3i 

- 5i  ■ f^  n   +  ^  aL  (192) 

The  suspension  force  is: 

and  the  jacking  component  is  given  by: 

PjK  =  -fe+^)  W^S^^te^ijf  "f-fl^  ^f  (195) 
Similarly  for  the  rear:    <-    =  3,  4 

^  -  ^^^^.-^.^^J^^  -^-;>^  ijfe*     (196) 

where  A  7"w£    ,  A  77te<:     are  the  half-track  changes  input  as  tabular 
functions  of  suspension  deflections. 

3.6       Braking  and  Tractive  Torques 

3.6.1     HV0SM-RD2  Version 

Braking  and  tractive  torques  in  this  program  versions  are  inter- 
polated from  input  front  and  rear  torque  tables  (  TQp      ,    TQ&      )  as 
functions  of  time. 
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3.6.2 


HV0SM-VD2  Version 


A  general  outline  of  the  procedure  used  for  generation  of  wheel 
torques  (both  tractive  and  braking)  is  depicted  in  Figure  3.8. 

The  form  of  the  open  loop  control  inputs  (tabulated  as  functions 
of  time)  for  speed  control  is  illustrated  in  Figure  3.9.   The  detailed 
algorithms  for  generating  wheel  tornues  is  duscussed  in  the  following 
sections. 

3.6.2.1   Braking  Torque 

The  brake  torques  at  the  individual  wheels,  (TQ)s'       >   are 
calculated  for  the  specified  brake  type  (i.e.,  1,  2,  3,  4)  based  on  analytical 
relationships  given  in  References  3  and  4  using  the  input  values  for  brake 
pushout  pressure,  Pj0  ,  and  the  instantaneous  values  of  brakeline 

pressure,   pj  ,  and  the  lining  fade  coefficient,  0-F)^      .   The 

following  relationships  are  applied: 

1.       ic      vs.  ■£    interpolated  from  direct  tabular  input. 


Pf  2  P 


C 


Pz  =      \ 


Fc        ,  -for     O^  Pc^=  ft 
Pi+Xlfe'ffi    ,  for    P^PC±P% 
3*Wft-P/>**ite-8J    ,   for  PC>PZ 


(197) 
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P-fW 


THROTTLE 
SETTING 


i 

i 

1  (TS) 

i 

• 

1 

1   1 

TORQUE  GENERATING 
ROUTINES 

PROPORTIONING 

VALVE 
CHARACTERISTICS 

1 
BRAKING      TRACTIVE 

1 

I 
1 

1 

BRAKE  TYPE(S) 
AND  PROPERTIES, 
TORQUE  VS.  PRESS., 
FADE 

:            : 

: 

■ 

TRANSMISSION 
RATIO 

(TR) 


ENGINE  TORQUE 

VS.  SPEED  FOR 

WIDE  OPEN  THROTTLE  (WOT) 

AND  CLOSED  THROTTLE  (CT) 


(TQ)B.(lb-ft) 
i  =  1,2,3,4 


(TQ)Dj(lb-ft) 
j  =  F,R 


Figure  3.8    GENERATION  OF  WHEEL  TORQUES 
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t  Pc (TS) (TR) 

Time  Master  Cylinder  Throttle  Setting  Transmission  Ratio 

Sec  Pressure,    PSIG  0. 00  <    TS     <1.00       (For  Disengaged 

0<    P       <  1000  .     __         wrYr  Clutch,  *  Enter 

—      c     —  1.00=    WOT  „n         n    nnn\ 

0.  00  =     CT 


TR     =  0.000) 


0.500  _    TR      _  5-000 


Exam  pi 

es :    • 

2 

.57 

1st 

1 

.83 

2nd 

1 

.00 

3rd 

0 

.70 

OD 

For  engine  RPM  <  500  during  a  brake  application  the  transmission 
ratio  is  set  to  zero.     The  effect  is  the  same  as  a  disengagement  of 
the  clutch.     Note  that  disengagement  by  the  driver  at  other  values 
of  engine  RPM  (e.g.,   between  gears  or  while  braking)  can  be 
simulated  by  directly  entering  a  value  of  0.  000  for     TR 


Figure  3.9        CONTROL  INPUTS  FOR  SPEED  CONTROL 
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3.     Brake  torques  are  calculated  from 


TYPE  1  BRAKE  -  Drum  Type  with  Leading  and  Trailing  Shoes 


[0 


(TQ)&r 


in 


for    (Pj'r  ?;,)■<  0 

-for    (Pj-Pjo)^0         (L6-FT) 


(Pj-fttbjGijtyafy 


(198) 


TYPE  2  BRAKE  -  Drum  Type  with  Two  Leading  Shoes 


(TQ)si  = 


O         ,     -for    lP,-Pi0)+0 


~Z     (Pj-  Pjo)  &4}   &1j 


Gi3  %  Mi 
U-6ZjG3j(LF)L 


,  tor  (Pj-Pjo)zO 
(L3-FT) 


(199) 


TYPE  3  BRAKE  -  Bendix  Duo  Servo 


H,  - 


0        ,     for    (Pj-Pj>)<0 


i  n 


TZW-frify 


U 


&Gj  +67j$i*$ 


&Si  6iij  &itj  Mi 


(ziaGmGihMl 


'tj  ^>8i  «7* j 


&9i&10<    1-&2j&l2i(LP)i 


'9^10)    '~^2j^i2j 


,    for(P/-Pj,)zO    (I*-*) 


(200) 
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fo&f  = 


TYPE  4  BRAKE  -  Caliper  Disk 


'    0  ,      for   (Pj-P^tO 

\   (P^Pio)  OajGyGtojUfy  ,  for  (P,-  -fo)  >0  (lb-pi)  (201) 

Where  &i j       ,    (?2j        ,    •-.  &3j         are  taken  from  the 
input  table  of  brake  data  (see  Figures  3.10  through  3.13, 
and  (LF)i        is  interpolated  from  the  brake  fade  table 
following  calculation  of  ti 


Because  the  brake  torques  must  always  oppose  wheel  rotations,  the 
algebraic'  signs  of  the  brake  torques  are  set  to  be  opposite  those  of  the 
wheel  rotations: 


Sgn(iGl6L    =  (r)Sqn.(£PS)L 


With  a  digital  integration  procedure,  the  simulation  of  coulomb 
friction  can  produce  instability  at  small  values  of  velocity.   Since  the 
friction  force  is  held  constant  during  a  calculation  increment,  it  can 
produce  an  overshoot  of  zero  velocity  and,  thereby,  add  energy  to  the 
simulated  system.  To  avoid  such  a  problem  in  the  present  simulation  of 
braking  torques  opposing  wheel  rotations,  a  sequence  of  calculation  logic 
has  been  developed  to  limit  the  applied  brake  torques  at  small  wheel 
rotational  velocities  (Figure  3.14).  The  analytical  basis  for  the  logic 
is  presented  in  the  following. 
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DRUM 
ROTATION 


6,     = 


<2         ~ 


(^/v 


*4 

Or 

P 


(**). 


8 


Lever  arm,   inches 

Actuation  constant,   assumed  to  be  equal  for  the  two 

shoes.      (Note  that     &z    m     1.42  in  Chrysler  Products,    Ref .  3  > 

Coefficient  to  permit  change  of  lining  friction  at  elevated 
temperatures 

Effective  lining -to -drum  friction  coefficient  at  design  temp. 

Cylinder  area  -  leading  shoe,   in. 

Cylinder  area  -  trailing  shoe,    in. 

Hydraulic  pressure,    psig 

Push-out  pressure,    psig 

0,  far  (P-Po)±0 

f)  X         Fw#  J       (  o) 


Tj(P-Po)  61^63  Up) 


tz 


03'FT) 


Figure  3.10   TYPE  1  BRAKE  -  DRUM  TYPE  WITH  LEADING  AND  TRAILING  SHOES, 
UNIFORM  OR  STEPPED  CYLINDER 
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DRUM 
ROTATION 


(P-po)  G+ 


(f>-f%)&4 


S7f 

(LF) 
©3 


Po 


(T*)B  - 


Lever  arm,   inches 

Actuation  constant 

Coefficient  to  permit  change  of  lining  friction 

at  elevated  temperatures 

Effective  lining -to -drum  friction  coefficient  at  design  temp. 

Cylinder  area,   in. 

Hydraulic  pressure,   psig 

Push-out  pressure,   psig 


'  0,  for  (P-f>)4  0 


l(P"P0)&1^4 


&z  &3(L/=) 
1-&z&3(LF) 


for0*(P-Po)    (L5-FT) 


Figure  3.11   TYPE  2  BRAKE  -  DRUM  TYPE  WITH  TWO  LEADING  SHOES,  TWO  CYLINDERS 
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DRUM 
ROTATION 


G7  ~  (/.  0+)(DR(JUA  radius) 


'  0,    -for  (P-P0)  40 


tok  = 


B 


Tz  (P~Po)  G* 


G6+G7  sing 

1  + 


Gf-G6-G7stng 

GtGaGr2(lF) 
09Gfo\l-GzGfz(LF)-\ 


G7  sin^  + 


,   for  0^fP-Po)  (L5-FT) 


Figure  3.12  TYPE  3  BRAKE  -  BENDIX  DUO  SERVO 
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sfcr  cc 


^/3 

P 


(TQl 


Total  cylinder  area  per  side  of  disc,    in. 

Mean  lining  radius,   inches 

Effective  lining-to-disc  friction  coefficient 

Coefficient  to  permit  change  of  lining  friction 
at  elevated  temperatures 

Hydraulic  pressure,    psig 

Push-out  pressure,    psig 

0,   for  (P-P0)  <0 


Figure  3. 13    TYPE  4  BRAKE  -  CALIPER  DISC 
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Solution  of  Equations    (44)   and    (45)    of  Section  3.1.4  for     (gP5)i 
and  (£P5)Ul     yields 

±  (ZPS\   =     X«J   [-fcik^(TQ)3^(ARl(TQ)Dj]t  z2j-2 [-£ ■„ hUl  -FclkL  +  11(TQBl - TZfaJ]  (202) 

di(m)li,=~ 1 ; 


•For  very  small  values  of  (lZP5)i     ,  (^.PS)lh       (i.e.,  smaller  than 
an  arbitrarily  specified  threshold  value  of   5g    radians/second)  it  is 
desired  to  prevent  the  applied  brake  torques  from  erroneously  accelerating 
the  wheels.  To  accomplish  this  objective,  the  values  of  (TQ)ql     ,    (TQ)&i+i 
are  limited  to  those  required  to  make  the  numerators  of  Equations  (202)  and 
(203)  equal  to  zero  (with  the  exception  of  the  tractive  term). 

From  Equation  (202),  the  limiting  value  of  brake  torque  is 


where  T   m} 


A„  =  +—± 


x   +  %W2 


IS 


Similarly,  from  Equation  (201),  the  limiting  value  of  brake  torque 
MuV,  -  -  fii4l  kLt1   -   -  7i3  (Fc-  kL  -  n  TOBl  )       .         (205) 
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The  value  of  T^-g       distinguishes  between  the  driving  and  non- 
driving  ends  of  the  vehicle  (i.e.,  for   XVi      =   0>   ^-3   =  0)  • 

The  need  for  the  inclusion  of  inertial  coupling  with  the  other 
wheel  of  a  driving  pair  is  illustrated  in  Figure  3.19.  The  traversal  of 
the  patch  of  ice  by  braked  wheel  I  can  reduce  the  circumferential  force, 

.  Pc  i  ,  permitting  the  brake  to  rapidly  decelerate  the  wheel  to  a  small 

rotational  velocity.   If  the  inertial  coupling  term  were  not  included  in 
Equation  (202) ,  the  value  of  (T^)g^        would  be  set  to  the  low  value  (  fc-L  ^  ) 
when  (RP5)i     became  smaller  than  5&      .   In  this  case,  however,  because  of 
the  inertial  coupling  terms  in  Equation  (202),  the  large  value  of  Fc  i+-i 
could  then  accelerate  wheel  i         in  reverse  rotation  unopposed  by  the  true 
brake  torque  applied  to  wheel  i        .   Therefore,  the  inclusion  of  inertial 
coupling  effects  in  the  determination  of  a  limiting  value  of  CfQ)&i     i-s  * 
seen  to  be  a  necessity. 


4r(VS\ 


&(**)i 


Figure  3.15     EFFECTS  OF  INERTIAL  COUPLING  OF  DRIVE  WHEELS 
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In  the  case  where  both  wheels  have  rotational  velocities  less  than 
the  specified  threshold  value  of  Se>        rad/sec,  the  simultaneous  satisfaction 
of  Equations  (204)  and  (205)  requires  that 


^a\i    ~    77ifcl  h'L  ,  and 


12 


(206) 


j 

Note  that  this  case  is  the  same  as  the  one  in  which    ±Dj      =   0. 

If  both  wheels  have  rotational  velocities  greater  than  J&      rad/sec, 
the  limitation  of  brake  torques  is  bypassed. 

The  effects  of  brake  fade  are  approximated  by  means  of  interpolation 
of  a  lining-to-drum  fade  coefficient,  (LF)i       ,  that  is  entered  as  a  tabular 
function  of  brake  temperature.  The  brake  temperatures  are  approximated 
by  the  following  relationships: 

£  =  ?U     *  titti)      ,  °F  C207) 

where 

A*  -  A-       ;  * 

&iSj  &I6J  (208) 


&t  =     ^~    -^J'*(tL-tA)     M       ,   stu     (209) 

1-^4         for  j~R 
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3.6.2.2   Tractive  Torque 

The  engine  torque  is  interpolated  from  tabular  entries  of  engine 
torque  vs.  engine  RPM,  for  the  given  throttle  setting,  as  shown  in  Figure  3.16, 
(The  throttle  setting  is  entered  as  a  tabular  function  of  time.)   It  is 
assumed  that  the  engine  torque  increases  linearly  with  throttle  setting 
between  the  values  entered  for  closed  throttle  and  for  wide  open  throttle. 
Thus,  the  engine  torque  is  obtained  from 


(TQ)E    =  (CT)+  (TS)    (WOT-CT) 


O^TS^IO 


(210) 


ENGINE 
TORQUE 
(TQ) 


WIDE  OPEN  THROTTLE  (WOT) 


4000 


CLOSED  THROTTLE    (CT) 
Figure  3.16     INPUTS  TO  DEFINE  ENGINE  PROPERTIES 
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The  engine  speed  in  revolutions  per  minute,  (fcPtf.E)    ,  is 
calculated  from 


f500     ,  or 


whichever  is  larger,  for  (fifty        ,  (RPSJtr/       at  driving  end  of  vehicle. 

To  permit  the  simulation  of  a  standing  start,  the  logic  depicted 
in  Figure  3.17  has  been  incorporated  in  the  program. 

3.7      Tire  Forces 

3.7.1     HV0SM-RD2  Version 

The  side,  braking  and  traction  forces  are,  of  course,  related  to 
the  tire  load  normal  to  the  plane  of  the  tire-terrain  contact  patch, 

Frl        ,  rather  than  the  radial  tire  load,  </z.L       .  Therefore  it  is 
necessary  to  find  the  value  of  f%.^      corresponding  to  the  radial  load, 

f~(Li        ,  and  the  side  force,  Fsi       .  The  components  of  the  external 
applied  forces,    t~/Zi       and  -  f~5L      ,  along  the  line  of  action  of  the 
radial  tire  force,  f~/zi         ,  are  depicted  in  Figure  3.18.  These  force 
components  must  be  in  equilibrium  with  F&L       »  such  that 

^I^-^cg-    *'Fs/M#i0c<k.    =*«;  (212) 
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NO 


YES  (EITHER  A  STANDING 
START  OR  A  BRAKE 
STOP) 


Is  10  psig  <  P 


NO 


YES 

(BRAKE  STOP) 


SET  (TR)  =  0 
(DECLUTCH) 


r 


INTERPOLATE  (TR) 
FROM  INPUT  TABLE 


(TQ)Dj  =  (TQ)E(TR) 


Figure  3.17    LOGIC  TO  PERMIT  A  STANDING  START 
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PLANE  OF 

TIRE-TERRAIN 

CONTACT 

/? 

PATCH 

1. 

A 


£6. 


4 


CAMBER  ANGLE 
RELATIVE  TO 
TERRAIN  (SHOWN 
NEGATIVE) 

SIDE  FORCE 
PRODUCED  BY 
COMBINATION  OF 
SLIP  ANGLE  AND 
CAMBER  ANGLE 

=  1,2,3.4 


Figure  3.18    VECTOR    SUMMATION    OF    FORCES    WITH    COMPONENTS    ALONG    THE    LINE 
OF  ACTION    OF    THE    RADIAL    TIRE    FORCE    (VIEWED    FROM    REAR) 


Solution  of  (209)  for   r>£   yields 


Fe:  "  F, 


«rFRL^0cc,.-Fs.i^0cGf_ 


(213) 


Since    f&L        is  required  for  the  determination  of    f-^i.      ,  an  initial 
approximation  of    f^£    i-s  obtained  by  extrapolation  from  the  previous 
time  increment.   Following  the  calculation  of  F$(_       in  the  current  time 
increment,  an  iterative  procedure  is  employed  to  correct  both  Ftz^       and 
PSL 
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Side  Forces 

The  side  force  calculations  are  based  on  the  small  angle  properties 
of  the  tires  which  are  "saturated"  at  large  angles.  Variations  in  the  small- 
angle  cornering  and  camber  stiffnesses  produced  by  changes  in  tire  loading 
are  approximated  by  parabolic  curves  fitted  to  experimental  data  (Figure  3.19). 
The  small-angle  cornering  stiffness  is  taken  to  be  the  partial  derivative  of 
lateral  force  with  respect  to  slip  angle  as  measured  at  zero  slip  angle  for 
various  tire  loads.  The  upper  plot  of  Figure  3,19  depicts  a  parabola  fitted 
to  the  small-angle  cornering  stiffness  as  a  function  of  tire  loading,  in  which 
the  cornering  stiffness  varies  as, 


*.0-Ao+Af^.-  -47-K) 


(214) 


The  lower  plot  of  Figure  3.19  depicts  a  similar  fit  to  small-angle  camber 
stiffness,  the  partial  derivative  of  lateral  force  with  respect  to  camber 
angle  as  measured  at  zero  camber  angle,  in  which  the  camber  stiffness 
varies  as, 

The  fitted  parabolic  curves  are  abandoned  for  tire  loading  in 
excess  of  J~Lj- A2       (see  Figure  3.19),  where  the  approximate  range  of  the 
coefficient  is  0.80  <  fZ-j-       <   1.15  (  JZT        is  an  adjustable  item  of 
input  data) .  The  side  force  properties  are  then  treated  as  being 
independent  of  tire  loading.  The  use  of  ftT        is  necessary  to  avoid 
artificial  reversal  of  the  slip  angle  forces  under  conditions  of  extreme 
loading  (i.e.,  where    Ax       «     rz-c      ).  Actual  properties  of  tires 
in  this  range  of  loading  (i.e.,  extreme  overload)-  are  not  known. 
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For  the  case  of  zero  traction  and  braking  the  side  forces,  for 
small  slip  and  camber  angles,  can  be  expressed  from  (215)  and  (214)  as 


(*J 


CAMBER. 


rCG>: 


(216) 


('«) 


SLIP 


AlFgt.(kF^-A2)-A0Ai 
A2 


Vs 


UA.        ri 


(217) 


where  the  sign  convention  corresponds  to  the  right-hand  rule  applied  to  the 
system  depicted  in  Figure  3.1. 

The  tire  model  must  handle  extremely  large  camber  angles  relative 
to  the  tire-terrain  contact  planes.  Applicable  tire  data  are  not  known  to  be 
available.   Therefore,  the  assumption  has  been  made  that  the  camber  force, 
for  a  given  normal  load,  will  reach  its  maximum  value  of  45  degrees  of  camber. 
In  accordance  with  this  assumption,  a  parabolic  variation  of  camber  force 
with  camber  angle  is  simulated  with  the  peak  occurring  at  45  degrees  (see 
Figure  3.20).  With  the  assumed  large-angle  camber  characteristic  depicted 
in  Figure  3.20,  Equation  (216),  for  the  complete  range  of  possible  camber 
angles,  becomes 


N 


CAMSfft. 


A. 


[*c«,-:f  4*J*caJ 


(218) 


To  permit  the  use  of  the  nondimensional  slip  angle  concept  which 
"saturates"  the  side  force  at  large  slip  angles,  an  "equivalent"  slip  angle 
is  defined  to  approximate  camber  effects 


fl!m\ 


A4[A|/5ii(Fi,-A,)-A.At] 


&•<- 4- £•,.!*., 


(219) 
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CAMBER 
ANGLE 
FORCE 


APPROX.  RANGE 
OF  AVAILABLE 
DATA 


f#.  =    CONST. 


45s 
CAMBER    ANGLE 


9(f 


Figure   3.20    ASSUMED  VARIATION  OF  CAMBER  FORCE  WITH  CAMBER  ANGLE 


Note  that  the  selected  analytical  treatment  of  camber  angles  subjects  the 
camber  force  to  the  saturation  effects  of  the  slip  angle,  superimposed  on 
the  assumed  behavior  shown  in  Figure  3.20.   While  the  assumption  depicted 
in  Figure  3.20  may  be  shown  to  be  in  error  when  appropriate  tire  data 
becomes  available,  it  was  found  to  be  necessary  to  reduce  the  "equivalent" 
slip  angle  of  large  camber  angles  to  avoid  an  unrealistic  predominance  of 
camber  effects  on  steeply  inclined  terrain  obstacles.   On  the  basis  of  the 
comparisons  of  predicted  and  experimental  responses  that  have  been  made 
to  date,  it  must  be  concluded  that  the  selected  analytical  representation  of 
large-angle  camber  effects  is  at  least  adequate. 
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Using  definition  (219),  the  resultant  side  force  for  small  slip 
angles  and  the  entire  range  of  camber  angles  can  be  expressed  as 


<- 


A  2. 


(220) 


Application  of  Equation  (220)  to  the  nondimensional  side  force  relationship 
(see  Figure  3.21)  yields 


m-)-7? 


tea 


-  A  -ifii  \fii 


+   27  A" 


(221) 


uw^ 


where    f^    =  resultant  side  force  for  entire  range  of  slip  and  camber" 
angles,  and 


A-- 


K 


(%)■ 


(222) 


■yVL4J& 


Large  values  of  the  slip  angles,  particularly  in  skidding,  make  it  necessary 
to  use  the  arctan  (  V&L /UCfL)   rather  than  (  Uei/xl&i) .      Also,  in  cases  where 
reversal  of  the  vehicle  velocity  has  occurred,  it  has  been  found  to  be 
necessary  to  control  the  algebraic  signs  of  the  slip  and  steer  angles. 
Equation  (222)  requires  that  Equation  (220)  be  modified  as  follows: 


A    = 


Aif*i(F*-At)-AA 


hz  (Hi)  mat       J  _ 


v& 


areian   -. —  -  f«n%)t'fj3i 


(223) 
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f(/S) 


?-Tfl/*hzT/>3,S"l/'(<3 

t  ,  for  3  <  J/3  J 


NONDIMENSIONAL 
SIDE  FORCE, 


BOUNDS  FOR  CUBIC  FORM  OF  fffi) 


1 


Figure  3-21     NONDIMENSIONAL  TIRE  SIDE-FORCE  CURVE 
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The  "friction  circle"  concept  is  based  on  the  assumption  that  the 
maximum  force  that  can  be  generated  by  the  tires  in  the  plane  of  the  tire- 
terrain  contact  patch  is  equal  in  all  directions.   With  the  use  of  the 
"friction  circle"  concept  (see  Figure  3.22),  the  maximum  side  force  can 
be  expressed  as 


Wu-/*w-tf 


(224) 


where   Fc-C        =  circumferential  tire  force  (i.e.,  traction  or  braking)  at 
wheel    i       ,  in  pounds  which  is  determined  from  interpolation  of  the  input 
tables  of  applied  torque. 


INCREASING 
TRACTION 


TRACTION    *  r 


INCREASING  . 
BRAKING      [ 


HEADING 


LOCKED  WHEEL  PRODUCES 
M  P«t  OPPOSING  VEHICLE 
MOTION 


^^utri'T^V'} 


SLIP  ANGLE  IN  TIRE  TERRAIN 
CONTACT  PLANE 


Figurel  3.  22      FRICTION  CIRCLE  CONCEPT 
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The  detailed  calculation  procedure  is  shown  below. 

The  determination  of  the  tire  side  force,  ii>L        ,  for  the  current 
time,  increment  requires  the  use  of   r^i        ,  and  a  value  of  f&C        i-s» 
therefore,  not  available  at  this  point  in  the  calculation.   For  this  reason 

F*C         from  the  previous  time  increment,  is  initially  used.  The 
subsequently  calculated  value  of   f<,±       is  then  used  to  recalculate  F^       , 
where 


FRl  -  Frl  Sec  cj>c&.L  -  HLt<*n(f>c&L 


The  following  tests  are  employed  to  determine  if  the  tire  has  lost 
contact  with  the  ground,  or  if  the  force  normal  to  the  ground,  F/z/         > 
is  calculated  to  be  a  negative  quantity  (this  condition  can  be  produced  by 
the  extrapolated  value  of  F$^         under  conditions  in  which  the  side  force 
reverses  signs) . 

If  FfLi         =  0,  or  (FfLi-F%$in-$a*i)   <_  0,  bypass  calculation  of 
the  tire  side  and  circumferential  forces,  set  f$L     =    '^i      =    '&£      =  0. 

The  circumferential  tire  forces  are  calculated  from  the  following 
relationships: 

1.  U        =   — /P~~  >  f°r  ..  P   =  1>  2  (front  wheels) 

r  lZ(r<8)e. 

>l       =   '  77 .        ,    for         L  =  3,   4    (rear  wheels) 

2.  A.      If        7~i      =  0,   bypass  calculation  of        fci      ,    set 

fee.      =0 

B.      If  0  <      72        (Traction) 

To  ,  for     72  ^  \/JnF&i  I 


(1)      Fc, 


(225) 
fhhL      ,  for     TL  >/jbLLFgJ 
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(2)  Is  any   Ti  > l/^L^J  ? 

(a)  If  no,  set    Fcl    ~  Fci 

(b)  If  yes,  compare  the  corresponding  product 
(  FcL  hi      )  for  the  other  wheel  of  the 
pair  (i.e.,  the  pairs  (1,  2)  and  (3,  4)  ) . 
Set  the  values  of  the  products  (  Fc^  lt[       ) 
in  the  given  pair  to  the  smaller  of  the 
individual  products.   For  example,  say 

Fcik-i       <      Fc,h,     ,    then    Fez.    =  '  Fc, 
and   Fc,        =   Fcz(^iMi)  .   (Effects 

of  differential  drive  gears.) 

C.   If  TL       <  0  (Braking), 

(1)   For  1.932  in/sec  <        JU&LI     , 


Fc  = 


4 


Tl(^nU6l)   ,  ior  ITU  ^  JjmF^  cos[arcian  j&  -  jfft^WjJJ 


(226) 


(2)      For     IU&i/      <   1.932   in/sec 

Calculate      Fc  i      as  above   (Equation  226) ,   but 
multiply  it  by 

lUeril 
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The  calculation  of  tire  side  force  includes  testing  as  to  whether 
effectively  all  of  the  available  force  (via  the  friction  circle  concept)  is 
taken  up  by  braking  or  tractive  forces,  and  whether  the  vehicle  is  in  motion. 


and 


Thus: 


If     / Pel / -  IjU-i rfii  j         -   1.0,   bypass  side  force  calculation, 

set  FsL  =   0 

If  both       I'V&il  and       I^Usol         <   0*5,   bypass   side  force 

calculation,    set         Fsl       =   0 

If  either     \V&i  \       or     \l\*i\       >  0.5,    and       I /5C  I      <  3, 


For  3<     ipLi        , 

%  -    M(Zu)^%       where       sgn  FSl  =  sSn  fc 

For  conditions  of  tire  loading  that  are  not  extreme  (i.e., 


A  ■ 


A  = 


UuMw  -  ft  J 


WVfiilft'i 


arcian 


z 


+pi  -  fynv-tOfl 


(227) 


(228) 


(229) 


(230) 


For  conditions  of  extreme  tire  overload  (i.e.,  f-g,    >  A2  Q.T      ), 


•2  ±Lr 


A 


MMr^-V 


LV^F^ 


archiri 


v*i 


\ul 


6i\ 


+P'L  -  (ynV-Jtc 


(231) 
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and 


3.7.2 


fi'i    = 


AtJW^TLri)-A0)\ 


HV0SM-VD2   Version 


*A*£  -  ~4 


re 


C&t 


Cbc 


(232) 


The  calculation  of  the  tire  loading  normal  to  the  ground,  Fg.L  , 

and  the  tire  side  force,    f^        ,  follows  the  same  assumptions  and 
derivations  given  in  Section  3.7.1.   However,  the  tire  model  employed  in 
this  program  version  makes  use  of  the  "friction  ellipse"  concept  in 
establishing  the  relationship  between  side  and  circumferential  forces. 
In  addition,  the  inclusion  of  wheelspin  as  a  degree-of-freedom  necessitates 
the  use  of  longitudinal  tire  characteristics  in  determining  the  instantaneous, 
circumferential  force  based  on  measured  tire  properties. 

The  "friction  ellipse"  is  an  extension  of  the  widely  used  "friction 
circle"  concept  (e.g.,  References  5  and  6)  that  permits  a  more  realistic 
analytical  treatment  of  interactions  between  the  circumferential  force 
(i.e.,  tractive  or  braking)  and  the  side  force  of  a  tire.   Experimental 
evidence  that  the  maximum  values  of  tire  friction  forces  are  dependent  on 
direction  relative  to  the  wheel  plane  is  presented  in  References  7  and  8. 
General  representation  of  frictional  properties  of  pneumatic  tires  requires 
independent  specification  of  lateral  and  circumferential  friction  coefficient 
and  their  variation  with  load  and  speed.   Further,  two  characteristic 
coefficients  are  required  to  represent  the  circumferential  friction,  i.e., 
the  peak  (  JLL^p     )  and  sliding  (  M-y.$   )  coefficient. 

In  the  "friction  ellipse"  form  of  treatment  of  interactions,  the 

maximum  value  of  the  resultant  tire  friction  force  in  the  tire-terrain  contact 

plane  is  assumed  to  be  bounded  by  an  ellipse  with  the  minor  axis  equal  to 
l/J-uFfi           and  major  axis  equal  to  2^/-Xf  <&          ,   where 
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effective  tire-terrain  friction  coefficient 
for  side  forces,  for  the  given  conditions 
of  vehicle  speed  and  tire  loading.   Note 
that,  in  the  absence  of  circumferential 
forces,  the  value  Li. u  F/z         constitutes  the 
maximum  achievable  side  force. 


M-Y-p  =  peak  tire  terrain  friction  coefficient  for 

circumferential  forces  for  the  given  conditions 
of  vehicle  speed  and  tire  loading. 

=  tire  loading  perpendicular  to  the  tire- 
terrain  contact  plane,  lbs. 

The  bounding  ellipse  is  depicted  in  Figure  3.23.   Note  that  a  value  of 

f^-Pp  =  fi-y         will  reduce  the  "friction  ellipse"  to  a  "friction 

circle". 

In  the  calculation  procedure  of  the  developed  analytical  treatment 

of  interactions,  the  circumferential  tire  force,  F  ,  is  given  first  priority 

in  utilization  of  the  available  friction.   The  maximum  value  of  side  force, 

(F  )   ,  corresponds  to  a  resultant  force  that  constitutes  a  radius  vector 
v  s  max       r 

of  the  bounding  ellipse,  is  then  determined  for  use  in  the  calculation  of 
side  forces  (see  Figure  3.23).  Note  that  this  sequence  of  calculations  is 
identical  with  those  of  References  4  and  5  for  the  "friction  circle" 
concept. 


In  Figure  3.23,  the  maximum  value  of  side  force,  (F  )   ,  for  a 

s  msix 

given  circumferential  force,  F  ,  is  determined  in  the  following  manner. 

Ojjnax       _i_        fc  -   (0Z  (233) 
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Solution  of   (233)   for   (F  )         yields 

s  max 


tew  -  {^f*r-*§ 


(-234) 


For  the  case  of  traction,  the  magnitude  of  the  circumferential  force  is 

limited  only  by  the  bounding  ellipse,  since  the  resultant  force  lies  in  the 

wheel  plane  under  the  condition  of  maximum  traction.  Therefore,  F  is 
r  c 


determined  directly  from 


(235) 


In  the  case  of  braking,  where  the  limiting  condition  on  the  resultant  force 
is  that  of  opposing  the  direction  of  vehicle  motion,  the  maximum  circum- 
ferential force  may  be  less  than  LLy,  Fz     •  Therefore,  the  circumferential 
component  corresponding  to  the  cited  limiting  condition  on  the  resultant 
force  is  calculated,  as  follows,  for  comparison  with  /A-y-  F&  ,  and  the 

smaller  of  the  two  values  (i.e.,  smaller  absolute  value)  is  used  in 
Equation  (234).   From  Equation  (233)  and  inspection  of  Figure  3.23, 


FcZ   +   «&,.-  /#«/+  fc 


1  - 


f± 

ML 


=  Ft  Sec2fi 


(236) 


Solution  of  (226)  for  F  yields 


&- 


Mrfl 


tan* 6  + 


7*% 


(237) 
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Thus,  the  calculation  of  F  ,  to  treat  either  traction  or  braking,  is  as 
follows: 

lAyi  tfi       ,  or 


Fc=     < 


mF/z- 


-Lanzfi  +  #-  (238) 


whichever  has  the  smaller  value.   Application  of  the  value  of  F  obtained 
from  Equation  (238)  in  Equation  (234)  yields  the  required  value  of  (F  ) 

S  IflclX 

The  determination  of  the  actual  value  of  the  circumferential  friction 
coefficient  /A./.  is  dependent  on  the  instantaneous  value  of  rotational 

slip,  SLIP.,  and  the  ^Ly.         vs.  SLIP  relationship.   This  relationship, 
as  shown  in  Figure  3.24,  is  in  turn  characterized  by  four  measurable 
quantities : 

Cf        =     the  slope  of  the  /J-y.     vs.  SLIP  curve  in 
the  small  SLIP  region 

/Ay-f>     =     the  maximum  value  of  the  circumferential 
friction 

SLIP     =  the  value  of  SLIP  at  which  f^-yp        occurs 

p-lb      =  the  value  of  /^x-      at  100%  SLIP 
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Figure  3,24   Normalized  Tractive  Force  vs.  SLIP  Model 


The  relationships  employed  to  define  the  entire 
are: 


jJ-y,     vs,  SLIP  characteristic 


=   CtJSLXP/  for         O^lSLlPl^SLlRj- 

=  CQtC,IUIPl+C1l$LXPIx        for     SLlPT^hLIPl<SLlPp 

=  ^tq/SLlP/tCslsaP/2         for     ^LTpp^lSLZPl<  1.0 


Cr 


(240) 
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where : 


SLIPT    = 


Cn       = 


Oriu 

-  .2.  /jy.P 


(*"'-  ^f 


cs 

= 

Pip  '  /^X5 

(5LXPp  -1.0)1- 

c4 

= 

-Zc5 

£3 

= 

/*■*.$  +-cs 

Cr   -       -2  cz  (SUPp) 


The  preceding  calculation  procedure  is  applied  for  all  values  of 
rotational  slip  less  than  that  which  produces  the  value  LL^p   .   For 
rotational  slips  greater  than  that  corresponding  to    U)cp   ,  Equations 
(234)  and  (238)  are  modified  by  the  substitution  of   /^-x   to  replace 

-/Aycp     .  The  effect  of  this  substitution  is  the  prevention  of  an 
excessive  recovery  of  side  force  capability  at  rotational  slips  greater  than 
that  at  which   jUx      =      M-xp    • 

The  detailed  calculational  procedure  follows. 

The  following  expression  is  used  for  the  calculation  of  tire 
loading  perpendicular  to  the  local  terrain: 


F*l    -    Fki  Sec  <f>c*L  -  FSL  tan.  <pc&L 


(239) 


The  determination  of  the  tire  side  force,  fc^  ,  for  the  current 

time  increment  requires  the  use  of  fkL        ,  and  a  value  of  F$i      .is, 

therefore,  not  available  at-  this  point  in  the  calculation.  For  this  reason, 
an  approximation  of  F$i      ,  extrapolated  from  the  previous  time  increment, 

is  initially  used.  The  subsequently  calculated  value  of  Pb£     is  compared 

with  the  extrapolated  value,  and  an  iterative  procedure  is  then  applied  to 
achieve  agreement. 
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The  following  tests  are  employed  to  determine  if  the  tire  has 
lost  contact  with  the  ground,  or  if  the  force  normal  to  the  ground,  f/^L 
is  calculated  to  be  a  negative  quantity  (this  condition  can  be  produced 
by  the  extrapolated  value  of  F$i        under  conditions  in  which  the  side 
force  reverses  signs) . 

If  F/z^  =  0,  or  (Fty—  /%£  sin  <frc&L  )  £  0,  bypass  calculati 
of  the  tire  side  and  circumferential  forces,  set  £$i  =  FcL  =  FfZ-i 
=  0. 

Rotational  Slip 

Let    ^GWL   =  U&L  COS  If/ +  V^SiyLljf- 
Test    U6>Wi, 

(1)  If   /  \XbWul   <  °-5  in/sec 
test    hL  (PP5)L 

(a)  £or\kL(&P$)L   I       <  0.5,   set      ($LtP)i     =  0, 

(b)  for  0.5   <    ll^i(HPS)iJ    set       CSLSP)L     = 
-1.00  sgn[U&Wl(&P$)L  J  • 

(2)  If  0.5   <    I  UbwJ 

(5LIp)i  =    ;  _     *£** 

For  1.000  <        IfcLXPJcl        ,    set      ($LLP)i       =  +  1.000, 
retaining  algebraic  sign. 


123 


Circumferential  Force 


(1)  For  (Pj-  Pjo)    1  0,  where   j'  =  F  for  L     =  1,2 

j  =  £  for   ^  =  3,4 

(2)  For  0  <  (Pj'Pjo) 

-psi.M'ifei(*S*'tt6VjL)        ,  or 


£s6  +  hxn" 


a  rata  n 


fe-^H 


(242) 


(245) 


whichever  has  the  smaller  absolute  value,  where 


£<;  =  i 


1 


1 


,U$lWcI^I(6lip)J 


J^  ,hr\(%I%\+\($ll?\\ 


(246) 


(/°si)im*x    =       A*ft  /,% 


^ 


iiSSfe  sgn    (SLIP.) 
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Side  Force 

The  calculation  of  tire  side  force  includes  testing  as  to  whether 
effectively  all  the  available  friction  force  is  taken  up  by  braking  or  tractive 
forces,  and  whether  the  vehicle  is  in  motion. 

Thus: 

(1)  If       eslFci     >   (     JUit  Flu      -   1.00),   bypass  side 
force  calculation,    set        fc,u       =   0. 

(2)  If  both      \Tqi\      and      \u&l\    <   0.5,   bypass   side 
force  calculation,    set        F$i       =  0. 


(3)  If  either  0.5   <_       I  V&J    or  0.5  <       \%gl\    >   test 

J>1        : 


(a)  for       //6j      <  3, 

hi  =  J /4( fid* -6*1%     \jk-  &$+  Lit? 


(247) 


(b)  for  3   <       \fiL\      , 


r*L  =     *i/A  (fkLf  -  &sL  Fcl        ,   anA   sgn  F6t  =  $qnf>L 


(248) 


For  conditions  of  tire   loading  that  are  not  extreme,    (i.e., 
Ff[L        <     AiHt      ), 


ft     * 


A%te;-^-v% 


(249) 
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and 


Pi     = 


?W  -  ^^J^J 


,(250) 


and 


For  conditions  of  extreme  tire  overload  (i.e.,    Fkc>  A?-Q-t         ) 


A  « 


#  - 


[arctav,  g_|+  ^  _  fyrL2Mi)W 


(251) 


V^t^--M*J 


.W^M-^'I^  -  *  *«*  l**l) 


(252) 


3.7.3     Components  of  Tire  Forces 

Finally,  resolution  of  the  resulting  forces  acting  on  the  tire, 
L     ,  into  components  along  the  vehicle  fixed  axis  system,  is  required 
for  application  to  the  differential  equations  for  motions  other  than  the 
wheel  rotations. 


F*HUi 

M  • 

COS  ot6g>; 

• 

^ytcc 
F*1«.i 

Kt 

cos  fen- 

k 

ces  2r62>. 

t 

rCZUc 
Fcuuc 

-•  Hi 

f 

II  **\  • 

tos<x.cc 
cos$c- 
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rcjui     J 

COS  Jet 

4 

(253) 


(254) 
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where 


cos  oc  •  ■ 

D,c 

ct 

K 

+  Dh 

+  DH 

C0S/8Ci  = 

Dt, 

H 

*»l- 

+  Dh 

COS*.    : 

I>3C 

(255) 


are  the  directional  cosines  of  a  line  perpendicular  to  the  normals  to  the 
ground  and  wheel  planes. 


^syuc 


=    fit 


COS  oCS{ 
COS#s; 

cos  *si 


(256) 


where 


COS  <%.= 


**i 


\K+*Z+<1 


cos  A 


COS     fy.    - 


K  +  K  *  *Z 


Cn 


'    K  -  K  *  *it 


(257) 


are  the  directional  cosines  of  a  line  perpendicular  to  the  normal  to  the 
ground  and   fCL,       (i.e.,   H    to  the  ground  plane). 

Further, 


cos  2rc 
cos  t&i,. 

cos/fCi 

COS  fa, 

€0 

*u  - 

cos  oce. 

COS  cC6£t. 

cos  fc . 
cos  ^ 
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COS  /Sc  .  COS  oc. 


C6s/3G#.  c°scc&£,. 


(258) 


are  the  directional  components  of  the  above  line. 

The  summation  of  the  tire  force  components  acting  along  their 
respective  vehicular  axes  are: 

fxut-    =  F/zxui  +  Fexui  +  Fszuc  (259) 


(261) 
F}*L    ~  ^tic    +   Fc%tii  +  ^u-C 

3.8       Impact  Forces 


This  subroutine,  used  in  the  HV0SM-RD2  version,  calculates  forces 
applied  directly  to  the  sprung  mass  for  simulated  collisions  of  the  vehicle 
with  horizontally  disposed  obstacles  having  flat,  vertical  faces  typified 
by  beam-type  guardrails  and  median  barriers,  or  walls.  Two  options  for 
simulating  different  types  of  barriers  are  included:   (1)  rigid  or  deformable 
barrier  option,  and  (2)  finite  or  infinite  barrier  vertical  dimensions 
option.   The  barrier  is  always  oriented  parallel  to  the  X'  space-fixed  axis 
and  is  also  treated  as  being  of  infinite  length. 

For  simulated  impacts  with  deformable  barriers,  inertial  effects 
and  curvature  of  the  barrier  are  not  considered.   The  deflection  of  the 
barrier  at  a  given  time  is  determined  by  incrementing  the  barrier  in  small 
discrete  steps  (A  u'&       ),  from  a  pseudo-deflected  position  that  results  in 
no  vehicle-barrier  interference,  until  a  barrier  deflection  is  found  for  which 
the  force  required  to  produce  that  deflection  is  balanced  with  the  force 
required  to  crush  the  vehicle  structure  in  the  region  of  mutual  interference. 
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Rigid-barrier  impact  forces  are  computed  in  a  similar  manner 
except  that,  in  this  case,  the  barrier  is  "incremented"  to  its  original, 
undeflected  position  in  determining  the  vehicle  crushed  volume  and,  thereby, 
the  forces  acting  on  the  sprung  mass. 

Provision  is  also1  made  for  the  calculation  of  point  forces  arising 
from  the  deformation  of  structural  "hard  points"  as  originally  documented 
by  the  Texas  Transportation  Institute  in  Reference  13.  The  undeformed 
location  of  three  hard  points  together  with  an  omni-directional  stiffness 
are  input  and  forces  calculated  based  on  the  amount  of  deformation  imposed 
on  these  points. 


3.8.1 


Vehicle  to  Barrier  Interface 


In  this  calculation  procedure,  the  analytic  geometry  of  the 
vehicle-barrier  interface  is  determined.  Like  the  curb  impact  subroutine, 
use  of  the  sprung  mass  impact  subroutine  is  determined  by  a  series  of  logical 
tests.   If  no  sprung  mass  impact  input  data  are  supplied,  the  subroutine  is 
not  activated  during  the  course  of  the  computations.  However,  if  input  data 
are  supplied,  a  test  for  excessive  spin-out  is  first  made  since  the  left 
rear  corner  of  the  vehicle  is  not  simulated. 

If  ty     >   135°  stop  calculations. 


barrier: 


The  corner  points  of  the  vehicle  are  tested  for  proximity  to  the 


x'c 

%Cf>K 

ycPn 

= 

y'c 

+ 

A 

• 

**» 

$CPn 

i'c 

%cpn 

(262) 


where  ~tc?n-       ,  %cp^.        »  2-cp 
the  space  fixed  axes,  and 


are  the  corner  points  with  respect  to 
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are  the  vehicle  corner  points  in  the  plane   2.  =  0  to  be  tested 


Using  the  value  of  YL     that  corresponds  to  the  largest  value  of 
^j'c^n,    ,  the  coordinates  of  the  top  and  bottom  points  on  that  corner  are 
calculated: 


ACPt 

*'< 

XCPn 

Vcrtr 

= 

y'c 

-h 

A 

• 

Vcpn 

3cpt 

rc 

%vT 

%CPb 

< 

^CPn 

y'eps 

= 

y'c 

-f- 

A 

• 

VcP-n. 

^CPb 

*c 

%V8 

and,  lAcprn.         =  the  larger  of  %'cp^ 


and 


D'cPb 
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(264) 


if    Z(Ay'e)  >  yipm-  £y'ah-i 

1.     bypass  remainder  of  sprung  mass  impact  subroutine 
for  current  time  increment, 

2,    set    Fn  ~  2  6c5  -£Fy*  -  Z^zs  ~  o 
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If  2 (A  ye    )  <  y'cPnt  -  (^  g   )^__,     ,  continue  with  calculation 

of  barrier  face  plane  constants  and  set  time  increments  to  A  tg, 

The  location  of  the  undeformed  hardpoints  with  respect  to  the 
initial  axis  system  are: 


^9rU 

xi 

!j$no 

= 

V'c 

Z'o 

■+■ 


II A II 


(265) 


The  direction  cosines  of  a  normal  to  the  barrier  face  plane  in  the 
vehicle-fixed  axes  are: 


COS   oLB 

0 

cos  /*e 

= 

Ar 

• 

1 

cos  J# 

0 

(266) 


And  the  point  of  intersection  of  the  ^   axis  and  the  barrier  face 


plane  is: 


*st. 

-*c 

yaL 

s 

\at 

• 

(y's\  -  y'c 

**i 

-t'c 

where  fa)£   -  (y^J-    n  '  fa',)  -  (i  -  /;  (&  y<B) 


(267) 
(268) 


barrier  cutting  plane,   i    =  1,  2,  ..., 
and 

YL         -   the  characteristic  of 


/I 


n.    +  1 


+  / 
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Knowing  the  direction  cosines  of  the  barrier  face  plane  normal  and 
a  point  on  the  barrier  face  plane,  (both  with  respect  to  the  vehicle  axis 
system)  a  constant  for  the  barrier  face  plane  is  determined: 


%dl  =  *&cos<*g   +  tj8l  cosfa   *■  i.5i  a>sfB 


(269) 


The  approximate  line  about  which  the  vehicle  rotates  is  determined 
by  the  intersection  of  the  previous  time  increment  barrier-vehicle 
equilibrium  plane,  and  the  present  barrier  plane  with  an  assumed  deflection 
of  (£&)t-i  A  plane  is  passed  through  this  line  perpendicular  to  the 

boundary  of  compression  of  the  vehicle  structure  are  found. 

The  coordinates  of  the  axis  of  rotation  at  the  plane   £   =(£-)&  t-> 


are: 


(cos  ocs)t  (cos  ^  (cos  *3)f 

(cos«B)^f         (cos^3)t_f       fees  *s)t_f 
0  0/ 


*8i 
If  81 


RB, 


(«*Xt 


(270) 


where,     Rgt   =  %3,  (cos  <x3)t      -h    VB,  (co* /9B)t    i-  %6,  (COS  >s)^  (271) 


and 


*S' 

Vs' 

=     \at     • 

(y'B)*_f  -  y'c 

$3' 

-*'c 

(272) 


A  constant  for  the  plane  perpendicular  to  the  present  barrier  face 
plane  and  containing  the  axis  of  rotation  is: 
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and  the  directional   cosines  of  a  normal   to  the  plane  are 


Cos   oc 


where 


and 


8/ 


COS   £3i    = 


Cos  XBf 
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<£"' 

c>« 
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a."'  =  h" (cos  >a)t    -   c"(cos/S3)t 

b  '»   =  +  c  "  (COS  cLs)t      -    *  "  (COS  *B)t 

c '"  =  +a"  (cos  /ge)t    -  Jb  "  (cos  a3)t 


<z 


'/      = 
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>'/       = 
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(cosfi8)t-t  (CosJrs)t-f 

(cos  *3)t  (c°s  <**)+ 

(cos  *a)t_f  (zoscts)t_f 

(co$«B)t  (C°s/30)t 

(cosocs)t_f  (cos/33)tf 


(273) 


(274) 


(275) 


Possible  intercepts  of  the  cutting  plane,       c      ,   on  the  vehicle 
periphery  are  then  given  as: 


COS    0£B 

A, 
A* 


cos  fo 


cos  Z3 

*"n 

*, 

c, 

• 

y*. 

= 

*i 

4 

lirt 

*2 

(276) 
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where      YL       =   1,    2,    3,    ...,    17,      See  Table  3-1  for  definitions  of     A,     , 
B,     ,     C,      ,      £,     ,    A 2      >      ^2    »     ^2-    >     /^i      corresponding  to      ^     . 

The  following  tests  are  employed  to  prevent   Equation    (276)    from 
becoming  indeterminate.      Note  that  this  condition  exists  when  the  barrier 
and  one  or  more  of  the  lines  or  planes  describing  the  vehicle  boundaries 
are  parallel. 

if  CC&cx-e,     =   °>    omit  points  4,    5,    10,    11 

if  COS  6g      =   °>   omit  points  1,    2,    7,   8 

if         CO6Y3       ~  ®*   omit  points  3,    6,    9 

if  COSO(B1COS^5-COS^61COSO(B  >   omit  points   12,    13 

if      U>$pBlC0SV&  =  a>sXB1  cosfe      >  omit  P°ints  14>   15 
if      C0SltBlai$G{a~a)S(X&1CO$X6      '   omit  Points   16,    17 

Note  that  points  to  be  omitted  refer  to  points  in  Table  3.1. 

Each  point  that  exceeds  the  following  ranges  is  rejected  as  being 
outside  the  vehicle  periphery: 

Xmr  ^   Xn-  ^  Xvf 

To  avoid  retention  of  improper  points  for  subsequent  calculation  of 
vehicle  crush  forces,  the  following  logic  is  employed.  The  plane  perpendicular 
to  the  barrier  face  plane  and  containing  the  axis  of  rotation  (the  fcgi 
plane)  is  a  compression  boundary  when  the  axis  of  rotation  is  located  within 
the  bounds  of  the  vehicle.   The  axis  of  rotation  is  inside  the  vehicle  if 
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Table  3-1 
DEFINITION  OF  ELEMENTS  IN  MATRIX  EQUATIONl(276) 


ITEM 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12* 

13 

14 

15 

16 

17 

Al 

1 

1 

1 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

1 

1 

0 

0 

Bl 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

.  1 

0 

0 

0 

0 

1 

1 

<M 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

•l 

XVF 

xVF 

xVF 

Yv 

YV 

YV 

XVR 

XVR 

XVF 

-YV 

-YV 

ZVT 

ZVB 

XVF 

XVR 

Yv 

"YV 

A2 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

(cos<*ai) 

B2 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

(cos$Bt) 

C2 

1 

1 

0 

1 

1 

0 

1 

1 

0 

1 

1 

(C0S*3() 

R2 

ZVT 

2VB 

YV 

ZVT 

ZVB 

XVR 

ZVT 

ZVB 

-YV 

ZVT 

ZVB 

(RB1) 

FRONT  o  3 


'OMIT  PTS   12-17  AT   INITIAL  CONTACT 
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The  intersection  of  the    £37  ,   Xv/r   ,  and  g       =  0  planes  is 


1.  If  the  axis  of  rotation  is  outside  the  vehicle,  omit  points 
12,  13,  14,  15,  16,  17. 

2.  If  the  axis  of  rotation  is  inside  the  vehicle  and  Y31VF  > 

Tjv      ,  reject  all  points  if  one  of  them  is  not  a  point 
12,  13,  14,  15,  16  or  17. 

If  less  than  three  points  are  retained,  let    i    -      i        +1  and 
return  to  calculation  of  the  point  of  intersection  of  the  ^j'     axis  and  the 
barrier  face  plane  (Equation  267) . 

If  option  2,  finite  vertical  barrier  dimensions,  is  desired  the 
preceding  calculations  are  carried  out,  and  it  is  necessary  to  obtain  the 
equations  of  the  top  and  bottom  barrier  planes  with  respect  to  the  vehicle 
axis  system  in  the  same  manner  as  for  the  barrier  face  plane. 


Thus 


COS  ocar 


Cos.  jr 


67 


-  Wat 


0 
0 

f 


{211) 


136 


are  the  direction  cosines  of  the  normal  to  the  barrier  top  and  bottom  planes 
with  respect  to  the  vehicle  axis  system,  and  a  point  on  the  barrier  top  plane 
is: 


"-Br 
Var 


A 


$3T  "  %C 


(278) 


A  point  on  the  barrier  bottom  plane  is: 


"X-aa 

i  -«j 

y*s 

— 

AT 

• 

-y'c 

?SB 

Ifse  -%c 

(279) 


Hence,    the  constants   for  the  barrier  top  and  bottom  planes   are 
*gr  =  *er  cos*8r  /  &8rcos/Sar  +  ^ar  cos  *sr 


(280) 


*BB   =    ^Bff  &#  *3T  +  VSB    COS  /Sgj-   i-  ?r3g  cos  Z3r 


(281) 


Determination  of  points  1,  2,  4,  5,  7,  8,  10,  11,  12,  13,  14,  15, 
16,  17  is  repeated  using  definitions  given  in  Table  3-2  except  as  follows: 


If  it       <  0,  omit  points  1,  2 


If  any  of  the  points  12-17  have  not  been  retained  from  Table  3-1, 
omit  the  corresponding. point  in  Table  3-2. 

If  both  points  1  and  2  from  Table  3-1  have  been  retained,  omit 
point  14. 

If  both  points  7  and  8  from  Table  3-1  have  been  retained,  omit 
point  15. 
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If  both  points  4  and  5  from  Table  3-1  have  been  retained,  omit 
point  16. 

If  both  points  10  and  11  from  Table  3-1  have  been  retained,  omit 
point  17. 

If  point  1  from  Table  3-1  has  been  retained,  use   ^^  =  ^gg   for 
point  14. 

If  point  8  from  Table  3-1  has  been  retained,  use   &i    =     K&r    for 
point  15. 

If  point  4  from  Table  3-1  has  been  retained,  use   K2    =     I^bb    f°r 
point  16. 

If  point  11  from  Table  3-1  has  been  retained,  use  K^    =  R.gj-    for 
point  17. 


Table  3-2 
DEFINITION  OF  ELEMENTS  IN  MATRIX  EQUATION  (276) 


ITEM 

1 

2 

4 

5 

7 

8 

10 

11 

12 

13   - 

1« 

15 

16 

17 

A1 

*1 

C1 

R1 
A2 

B2 

C2 
R2 

1 

0 
0 

XVF 

COS&  gj 

cos/f  BT 

COSt^   gj 

RBT 

1 

0 
0 

XVF 

0 
1 
0 

Yv 

0 
1 
0 

YV 

1 

0 
0 

XVR 

1 

0 
.0 

XVR 

0 

1 
0 

"Yv 

0 
1 
0 

"YV 

coso<bi 
cos^B1 

COS    Vo] 

RBI 

RBB 

RBT 

RBB 

RBT 

RBB 

RBT 

RBB 

RBT 

RBB 

RBT 

RBB 

RBT 

RBB 
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The  following  tests  are  also  employed  for  use  with  Table   3-2  to 
prevent  the  solution  to  Equation    (276)    from  becoming  indeterminate  with  the 
above  points: 

if    COS#BCO5X6T=(0SV5  COStXer  >   omit  Points  4,    5,    10,    11, 

if     COSfoCfiSfigj-^  CoSfaCDSltf  >   omit  Points   1,    2,    7,    8. 

if    coscx&[. 0)5/361  cosVqt  -  cospmcos  hli  +  tosoifrLtospb COS %-  COSfa  CDS  Xg] 

=  <2»(Xgf[a»^8064r-C^8r<S»4]    ,  omit  points  12,  13,  14,  15, 
16,  17, 

Again  reject  points  with  values  outside  the  ranges: 

Compare  values  of  \i.n,\        in  both  sets  of  points  retained,  and 
where  two  exist  for  a  given   n      ,  reject  the  point  with  the  larger  value 
of  /Z*./ 

If  less  than  three  points  are  retained,  let    L      =  L        +1  and 
return  to  the  calculation  of  the  point  of  intersection  of  the  Zj'      axis  and 
the  barrier  face  plane. 

For  each  point  retained,  calculate  its  velocity  components: 


u 


V. 


to- 


ri 
t 

Tl 
I 


n 


(282) 
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Of  the  remaining  points,  retain  only  the  four  points  with  the 
larger  values  of   ifyi      •   (Note  that  points  3,  6  and  9  are  counted  as  two 
points  in  this  selection.) 

If  a  rigid  barrier  impact  is  desired  and  the  values  of  Un.       for 
the  four  points  retained  are  negative,  bypass  the  remainder  of  the  sprung 
mass  impact  routine  for  the  current  time  increment, 


Set      ?#    ~  L^zs  ~  L>^s~^Fis  ~° 


3.8.2 


Area  Calculation 


In  this  calculation,  the  contact  area  of  the  interface  between  the 
vehicle  and  the  individual  barrier  cutting  planes,  i       ,  will  be  determined, 
as  illustrated  in  Figure  3.25. 


/  //   //  //  \ 


(*;>*;.  4 


fa  9n  4- 


(V.*;>4 


Figure  3.25  INTERSECTION  AREA 


140 


Given  the  three,  or  four,  points  that  define  the  intersection  area, 
(x"j   »  y'j    y     2.f);  >  where  j       =  1,  2,  3,  (4),  make  the  following 
selections:   (Note  that   I       identifies  the  cutting  plane.) 

Of  the  two  larger  values  of  X>j       ,  select  the  point  with  the 

"j      ,      Let  this  point  be  designated,  within  the  subroutine, 
■  i",        ).   For  the,  case  of  equal  values  of  x 


larger  value  of 
as  ( 


A 


y"i 


j 


scan 


all  points  for  the  larger  value  of   £ 
select  either  one. 


/ 


If  2  values  of  b"i      are  equal, 


Of  the  two  smaller  values  of   X/   ,  select  the  point  with  the 
smaller  value  of  £.'j     .      Let  this  point  be  designated,  within  this  subroutine, 

as  c    xi   'i   y 2 


£2.   )•   For  the  case  of  equal  values  of  X. 


scan 


all  points  for  smaller  value  of  £y 
select  either  one. 


If  2  values  of  2.-^ 


are  equal, 


Let  one  of  the  remaining  points  be  designated,  within  this 
subroutine,  as  (   ^'3  ,   y£   ,  £^   )• 

If  four  points,  identify  the  remaining  point,  within  this 
subroutine,  as  (  x'lj.     ,  *j"+      ,     2-'j-      )• 

Direction  components  of  lines  in  the  cutting  plane,    i        ,  J_ 
to  the  line  between  (  X-",      >     y'j       >     zj  .  )    and  (  x"t     >      y"z      >     i"z      ) 
are: 


d 


b'c 


/     _ 


(y"r  -  Vkh 

cos  /f3 
cos  X8 
cos  <x3 


cos  ^s 

COS  UB 

(r,  -  vi)i 

cos  /Ss 


(283) 
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H 


The  length  of  line   S?    as  shown  in  Figure  3.25  is 


& 


(284) 


The  length  of   S3   as  shown  in  Figure  3.25  is: 

If  three  points  were  retained,  set    \$>%)i      =  0. 

If  four  points  were  retained,  the  length  of   5^   as  shown  in 
Figure  3.26  is: 

(s2).  .  \«(<-**<)i+XM-Mi*<(K-rt\  (286) 

The  total  intersection  area  is  then: 


M.-(^)fV< 


(287) 
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f  if     "    i A 

[^3>y 3^3)1 


(*;.  yl  i")c 


fryltljt 


Figure  3.26        CENTROID  OF  INTERSECTION  AREA 


Centroid  Calculation 

In  this  calculation,  coordinates  will  be  determined  for  the 
intercept,  in  the  actual  vehicle-barrier  interface  plane,  of  a  normal  through 
the  centroid  of  an  individual  intersection  area. 

Calculate  the  following:   (Subletter  l      omitted  for 
convenience) 


a 


a 


'01     - 


xt  *  *; 


-  *J    ,         *ct    =    */  ******  ~cct  cos^s 
bc1    =     y/       *'  -  y%  ,        h'cf    =    eef  cos<x3  -  acr  cos  *a 


ft   +  $"*      ** 


-  >j  '       4y 


^  d*V  /f^   -  ^  COS  ois 


1    -,/f     .      //     ,,// 


t      ~" 


&,     =     <z„  %,%  +  b'cf  y'J    /  c'c/  ft 


(288) 
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XJ  *   < 


UCZ       - 


'C2       = 


-  Z 


2     » 


*;  *  y/ 


-* 


'CZ 


z"      +4." 
^ 2.2    » 


ft 


-   4*  ^     *  4 


u<rz 


4« 


cz   2  2 


bC2   cos  Js  -  ccz  cos  /S3 

Ce2    COS  CtB    -  Acz   COS  Xg 


=     a 


k  cos  /S5   -bcz  cos  «B 


(289) 


Find  the  centroid  of  the  upper  triangle  in  Figure  3.26  (i.e., 
point  P      )  by  solving  the  following: 


COS  Q£B 


a 


Cf 

/ 

CZ 


cos  0B 
Kr 


bcz 


cos  #8 

• 

= 

*L 

*p 

Gz 

(290) 


(Note  that  (£.%)±-t  corresponds  to  the  actual  equilibrium 

interface  plane  from  the  previous  time  increment.) 

If  3  points  only,  set 


(y«)£ 
(Mi 


'  (zp)c- 


If  4  points,  calculate  the  following: 


z;  /   *z 
a,,  =  -1 £  _  *£    , 


"C3    = 


CC3    - 


2 

f, 

-t- 

Vt 

*- 

6 

+ 

*; 

aC3      -    ^J   ^  ^ 


<?5 


cos  /2t 


-  y*  » 


^£3       =      ^£3     ^  °^     -  &■/•*     COS  & 


■z" 
2  +     > 


CC3      -       *«f   ***"A     -4?    ^ 


«:. 


%   -  <,  <  ■/.  ^  <  v-  ^  *; 


(291) 
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£. 


Of 


x"  +    r" 

+   -x" 


cc+  = 


z 

»; 

2 

it  * 

r< 

2    > 


AG+      =       *>€+   C0S  <*8    -  Cc+  C4S/?S 


bc4    -     <bf   cos  o£B  -  aCf  cos  Zg 


C4-      ~  ac4.    C0S/&6    ~  *>C4.    C0$    otA 


Find  the  centroid  of  the  lower  triangle  in  Figure  3,26  (i.e., 
point  <3   )  by  solving  the  following: 


COS    ois         COS  /Sa 


a 


C3 

04 


Cas  £g 


sc3 


"c+ 


H 

• 

fy 

— 

a3 

h 

$+ 

(292) 


(293) 


(Note  that  (Z$\^         corresponds  to  the  actual  equilibrium 
interface  plane  from  the  previous  time  increment.) 


The  centroid  of  the  entire  area  is  obtained  from  the  following 
relationships : 


(>A 


s,  +  sz)c 


(h  -  *r)i 


(294) 


In  this  calculation,  the  total  force  normal  to  the  barrier  will  be 
calculated  and  the  coordinates  of  the  application  point  for  an  equivalent 
concentrated  load  will  be  determined. 
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The  point  of  application  may  be  obtained  from: 


£4  = 


n+1 

t:1 


2 


%..  +  X 


*i-f 


L 


(A„r\  *  (*,„){., 


m= 


(*,„),   +  ft»r)<_, 


y*i  *  y*<-r 


w- 


L 

1*1 


n'-fl 

L 


\Attrr)i  +  (Attn-)  i-t 


\AwT)i  +   (^wr)i-/ 


f*i    +  **/„, 


L 


(Amrh    +  (Ai*r)i-t 


C295) 


(296) 


(297) 


The  total  force  normal  to  the  barrier  due  to  vehicle  crush  is 


given  by: 


n'*i 


H    -Ky(W,)l 


(Amt)l  +  (A/A/r)i-/ 


(298) 


The  force  normal  to  the  barrier  due  to  deformation  of  the  hard 
points  is  assumed  to  be  zero  if  the  lateral  position  of  the  undeformed  hard 
point  becomes  less  than  that  of  the  barrier  or  if  the  lateral  velocity  of 
the  hard  point  becomes  less  than  zero.  That  is 
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where : 


LI 


VI 


V*L 


W. 


#i 


=  UA 


(299) 


However,    if  y-^    _>  0  and         y$T6 '    >     "Ho 


(300) 


Note  that  the  position  of  the  deformed  hard  point  in  the  vehicle 


system  is: 


XfTt 


=  il/ll- 


*<fTi     ~    *C 


a. 
And  the  total  vehicle  force  is     FN   =  (Fn)+     +   -3  FnstL 


i 

i- 1 


(301) 
(302) 


If  the  barrier  is  rigid,  let    I      =    l     +   1  and  return  to  the 
calculation  of  the  point  of  intersection  of  the  y'     axis  and  the  barrier 
face  plane  until  (  y$)±        =(  M&   \  ••  Tften  proceed  to  the  calculation 
of  sprung  mass  resultant  forces  and  moments  (Section  3.8.4). 

If  the  barrier  is  deformable, 


1. 


let  the  assumed  deflection  of  the  barrier  be 


&  =  (TfsJi  ~  (y'ejo  -  £«. 
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where   Sn.        is  the  permanent  set  of  the  barrier 
(initially  zero) , 

2.     calculate  the  barrier  force,  Fq     ,  as  described 
in  the  next  section  (Section  3.8.3) 

If  the  vehicle  normal  force  and  the  barrier  normal  force  are  not  in 
balance  (within  limit,  s^      ),  i.e., 

£5  <  (fa  -Fn) 

let    L     =         l,     +   1,  and  return  to  the  calculation  of  the  point  of  inter- 
section of  the.  Tj'     axis  and  the  barrier  face  plane. 

If   *=&    >  (   Fg  -  F/4      )»  proceed  to  the  calculation  of  sprung 
mass  resultant  forces  and  moments  (Section  3,8.4). 


3.8.3     Simulation  of  Deformable  Barrier 
Load-Deflection  Characteristics 


A  previously  developed  general  subroutine  for  handling  non-linear 
load-deflection  characteristics  (Reference  12)  has  been  adapted  for 
calculating  barrier  forces  for  simulation  of  sprung  mass  impacts  with 
deformable  barriers.   Use  of  the  subroutine  requires  that  the  load-deflection 
characteristics  of  the  simulated  barrier  be  known. 

The  load-deflection  relationship  for  increasing  load  is  represented 
by  a  general  polynomial  function  of  deflection  and  the  velocity  of  deflection. 
(Velocity  terms,  included  in  the  original  application  for  which  this  subroutine 
was  developed,  have  been  retained.   However,  since  simulated  barrier 
deflections  occur  in  discrete  steps  of  A  tj'q        ,    the  coefficients  of  the 
velocity  terms  in  the  polynomial  for  the  load-deflection  relationship  should 
be  entered  as  zero  for  the  present  application  to  barriers.)   For  decreasing 
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load,  the  load-deflection  relationship  is  represented  by  a  parabolic 
function  of  the  deflection  that  is  determined  from  specified  (input  data) 
ratios  of:   (1)  conserved  energy  to  total  absorbed  energy,  and  (2)  residual 
deflection  to  maximum  deflection  (see  Figure  3.27). 

When  recurrent  loading  takes  place,  the  residual  deformation  for 
each  cycle  is  incorporated,  in  the  form  of  initial  surface  displacement,  for 
subsequent  cycles..  If  the  unloading  between  cycles  is  incomplete,  a  value 
of  residual  deformation  is  determined  such  that  the  new  loading  curve  will 
pass  through  the  final  point  on  the  preceding  unloading  curve  (see 
Figure  3.28). 

The  general  procedure  followed  in  the  nonlinear  load-deflection 
subroutine- is  described  briefly  in  the  following  paragraphs. 

1.  The  deflection  and  the  velocity  of  deflection  are 
determined  from  analytic  geometry  and  applied  in  the 
appropriate  functional  relationship  to  determine 
the  corresponding  load. 

2.  The  change  in  area  under  the  force  deflection  curve, 

(  i-t~i~-t-])(.ot~°t-\) /*-   is  calculated  and  stored 
for  each  time  increment. 

3.  When  a  maximum  load  point  is  reached  (as  indicated  by 

an  initial  negative  value  for  the  velocity  of  deflection) , 
the  total  area  under  the  loading  curve  (total  absorbed 
energy)  is  calculated,  and  the  coefficients  for  a 
parabolic  unloading  curve  are  determined  from  input 
specifications  for  (1)  the  ratio  of  conserved  energy 
to  total  absorbed  energy,  and  (2)  the  ratio  of  residual 
deflection  to  maximum  deflection. 
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DEFLECTION  AT 
WHICH      $<-0 


■    Sfinol  O, 


«w» 


DEFLECTION 

8 

Figure  3.  27  GENERAL  FORM  OF  BARRIER  LOAD-DEFLECTION  CHARACTERISTICS 
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&>«■*-  ^ 

t 

FORCE 

'  ^ 

*t 

/^-K- 

-^ —  RESIDUAL 

DEFLECTION  -»J 

FROM  FIRST 

LOADING 

CYCLE 

DEFLECTION,  S& 
(a)      SECOND  CYCLE  STARTS  PRIOR  TO  COMPLETE  UNLOADING 


FORCE 


H!_  RESIDUAL  DEFLECTION  _J 
FROM  FIRST  LOADING  CYCLE! 

DEFLECTION,  £ 


(b)     SECOND  CYCLE  STARTS  AFTER  COMPLETE  UNLOADING 


Figure    3, 28        ASSUMED  FORM  OF  LOAD-DEFLECTION  FOR  RECURRENT  LOADING 
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4.     When  the  deflection,  during  the  unloading,  becomes 

equal  to  or  less  than  the  specified  residual  deflection 
(from  the  calculated  maximum  deflection  and  the 
specified  ratio  of  residual  to  maximum  deflection) ,  the 
force  is  set  equal  to  zero. 

The  logic  and  calculation  steps  for  the  general  nonlinear  load- 
deflection  subroutine  are  presented  in  Figure  3.29.  Note  that  the 
correspondence  between  the  symbols  used  in  Figure  3.29  and  those  defined 
in  this  report  for  application  to  deformable  barriers  is  given  in  the 
following  table. 
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In  the  logic  depicted  in  Figure  3.29,  the  criteria  for  a  "proper" 
range  of  combinations  of  the  input  parameters  that  control  unloading 


characteristics  (i.e.,  SET  and  CONS)  are  based  on  an  adopted  requirement 
that  unloading  curves  (i.e.,  load  vs.  deflection  for  decreasing  deflection) 
should  be  concave  upward,  with  the  linear  case  as  an  upper  limit  and  zero 
minimum  slope  as  a  lower  limit.   During  the  original  development  of  the 
general  nonlinear  load-deflection  subroutines,  it  was  found  that  a  dynamic 
load-deflection  characteristic  for  increasing  deflection  that  was  initially 
quasi-linear  but  then  yielded  at  a  constant  force  level  required  exploratory 
variation  of  the  input  parameters  that  control  unloading  characteristics, 
since  the  "proper"  range  of  combinations  was  dependent  upon  the  maximum 
deflection.   Further  difficulty  was  encountered  in  the  case  of  multiple 
loading  cycles,  in  which  the  range  of  maximum  deflections  for  the  individual 
loading  cycles  are  generally  different. 

In  order  to  avoid  the  above  difficulties  while  still  retaining 
control  of  unloading  curve  characteristics,  the  depicted  program  logic 


automatically  adjusts  the  input  values  of  CONS,  as  required  in  initial  loading 


cycles.   Any  adjustments  made  in  the  input  values  of  CONS  are,  of  course, 
indicated  in  the  program  output.   In  addition,  whenever  a  sequence  of  loading 
cycles  leads  to  an  "improper"  unloading  curve,  either  through  the  mechanism 
of  different  ranges  of  maximum  deflections  or  through  the  accumulation  of 
"carry-over"  energy  from  previous  loading  cycles  with  incomplete  unloading, 
the  revised  program  logic  resets  the  conserved  energy  to  a  value  corresponding 
to  the  exceeded  limiting  case  of  unloading  curve  characteristics.  Again, 
the  program  output  indicates  any  adjustments  that  are  made. 
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3.8.4     Resultant  Forces  and  Moments 


The  velocity  components  of  the  points  of  application  of  the 
sprung  mass  forces  are: 

\lk     =  IIA||'   r't.tL.i-  K? 


where : 


y> " 


(303) 
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And  its  velocity  tangential  to  the  barrier  is: 


VTAA/C   =  [)/ (u'J  +  (ur$*\(sgn   <) 


(305) 
The  friction  force  acting  on  the  vehicle  due  to  contact  with  the 


barrier  is 
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(306) 
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The  resultant  force  components  along  the  vehicle  fixed  axes 
acting  on  the  sprung  mass  are: 


where, 


and 


Hfsi 

~ /*B  fetJi   cos6Bl 

EpL 

=  1M1- 

-   ^sJl 

h*L 

-  yfiFsNL  $'"£bl 

&9c 

=    ancz&rt 

*iHi 

-    (ft 

h 

hhiL 

-   Fm 

for       i 

2,    3,    4. 

(307) 


The  moments  acting  on  the  vehicle  sprung  mass  are,  therefore, 

2N^  =    Z[  $F2Sl   -  zL  Fy5  ■]  (308) 

(309) 
(310) 


3.9 


Preview-Predictor  Driver  Model 


The  driver  model  includes  several  modes  of  operation:  path- 
following,  speed  maintenance,  speed  change  and  skid  recovery  modes 
(Figure  3.30).  The  normal  steering  mode  is  path-following.  The  path- 
following  mechanism  predicts  the  future  vehicle  path  based  on  the  instan- 
taneous velocity  and  an  estimated  lateral  acceleration  due  to  the  instantaneous 
front  wheel  steer  angle.   In  the  calculation  of  the  projected  path,  the  model 
assumes  that  the  vehicle  will  maintain  its  present  velocity  except  for  the 
continuous  effect  of  the  lateral  acceleration.  The  estimated  lateral  accel- 
eration is  a  function  of  the  forward  velocity,  the  front  wheel  steer  angle 
and  the  stability  and  control  characteristics  of  the  vehicle.   Error  estimates 
are  made  between  the  predicted  path  and  the  desired  path  at  event ly  spaced 
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Figure  3,30     MODEL  OUTLJNE 
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points  along  the  predicted  path.  These  error  estimates  are  weighted  to 
account  for  the  reduction  in  lateral  acceleration  required  to  null  errors  at 
further  distances  ahead  of  the  vehicle.  The  error  estimates  are  also 
weighted  to  emphasize  the  importance  of  errors  at  particular  locations 
along  the  desired  path.   The  change  in  the  front  wheel  steer  angle  is 
proportional  to  the  average  of  the  weighted  error  estimates.  This  steer 
command  is  filtered  to  reflect  human  dynamic  capabilities,  specifically 
neuro-muscular  dynamic  characteristics. 

Operating  simultaneously  with  the  path-following  mode  is  either  the 
speed  change  .mode  or  the  speed  maintenance  mode.  The  speed  change  mode 
uses  input  data  values  of  the  desired  speeds,  their  initiation  times  and 
attainment  distances  to  determine  the  required  brake  pedal  forces  and 
accelerator  pedal  deflections.   In  calculating  the  brake  and  throttle  control 
commands,  the  model  assumes  linear  relationships  between  deceleration  and 
brake  pedal  force  and  acceleration  and  throttle  position.   Initialization 
of  the  path-following  mode  requires  that' the  initial  vehicle  position  and 
heading  be  consistent  with  the  desired  path.  To  initialize  the  speed  control 
mode  the  initial  vehicle  speed  is  set  equal  to  the  first  desired  speed. 

If,  because  of  terrain  features  or  cornering  forces,  the  vehicle 
should  gain  or  lose  speed  such  that  the  difference  between  the  desired  speed 
and  the  actual  speed  exceeds  the  threshold  values,  the  model  will  activate 
the  speed  maintenance  mode  and  attempt  to  return  the  vehicle  to  its  most 
recent  desired  speed.   These  threshold  values  are  variable  model  inputs. 

The  skid  recovery  mode  is  activated  if  the  vehicle  slip  angle, 
Oc    ,  the  angle  between  the  vehicle  heading  and  its  velocity  vector, 
exceeds  a  pre-selected  threshold  value,  Tizi       .  The  degree  of  severity  of 
the  skid  is  determined  from  comparison  of  the  vehicle  slip  angle  with  a 
second  (higher)  threshold,   7/22.  •   For  skids  of  low  severity  the  brake 
pedal  force  (  /^/e^  )  an(i  accelerator  pedal  deflection  (  AP£>    )  are  set  to 
zero;  however,  the  steering  control  remains  under  the  path-following  mode. 
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If  the  skid  is  of  high  severity,  in  which  the  higher  vehicle  slip  angle 
threshold  is  exceeded,  the  steering  commands  are  determined  by  the  skid 
recovery  routine.  Thus,  for  more  severe  skids,  brake  pedal  force  and 
accelerator  pedal  deflection  are  set  to  zero,  path-following  is  discontinued 
and  steer  commands  are  generated  strictly  to  recover  the  directional  control 
of  the  vehicle.  Under  the  skid  recovery  mode,  the  steer  angle  change  is 
proportional  to  the  average  front  wheel  slip  angle,  Ijfp     ,  and  the  sign  of 
the  steer  correlation  is  in  the  direction  of  reducing  the  average  front 
wheel  slip  angle.   Repetitive  steer  angle  correlations  are  made  until  the 
wheel  slip  angle  is  equal  to  zero.   Steer  commands  are  filtered,  as  in  the 
path-following  mode,  before  being  applied  to  the  vehicle. 

The  data  sampling  scheme  is  virtually  identical  to  the  one 
investigated  by  Kriefeld  (Reference  14).  All  of  the  model  functions 
operate  on  their  appropriate  inputs  at  discrete  time  increments,  as 
determined  by  the  input  value  of  EMDT,  the  time  between  driver  samples 
in  seconds.   Significant  changes  in  model  output  can  be  produced  by  this 
mechanism,  including  improved  correlation  with  recorded  nonlinear  responses 
of  human  operators. 

It  should  be  noted  that  the  thresholds  or  indifference  levels 
mentioned  above  are  designed  to  be  single  parameters  representing  the 
minimum  detection  level  for  that  particular  control  input  or,  if  the  driver 
chooses  not  to  act  until  a  higher  value  is  reached,  the  minimum  indifference 
level  for  that  control  input. 


3.9.1 


Path-Following  Mode 


3.9.1.1  Specification  of  Desired  Path 

The  center  line  of  the  desired  path  is  entered  into  the  program 
by  a  series  of  continuous  straight  line  segments,  i.e.,  first  order 
polynomials. 

3.9.1.2  Determination  of  Predicted  Path 


The  predicted  path  (i.e.,  driver  prediction)  is  obtained  by 
integration,  from  the  present  vehicle  position,  of  the  instantaneous  velocity 
and  estimated  lateral  acceleration.  The  estimated  lateral  acceleration  due 
to  the  front  wheel  steer  angle,  Wpmiz,       ,   can  be  calculated  by 


"r    I 


3 


idea./ 


(*'b)(l+KduTz) 


(311) 


where  0-      and  b      are  the  distances  along  the  vehicle-fixed  X     axis 
from  the  sprung  mass  center  of  gravity  to  the  center  lines  of  the  front  and 
rear  wheels,  respectively,  inches,  U-r      -   magnitude  of  the  vehicle  velocity 
and  Kck         is  a  performance  parameter  characterizing  the  understeer- 
oversteer  properties  of  the  vehicle.   Rewriting  the  lateral  acceleration  in 
the  inertial  reference  system  yields 
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where  jjAlj      is  the  transformation  matrix  from  the  vehicle  fixed  reference 
system  to  the  space  fixed  system.  The  coordinates  (  XuPi     >    H'vp/      )  at 
evenly  spaced  points  a  distance  A6>    apart  along  the  predicted  path  can 
be  found  by 


VP, 


=    K  * 


u 


OL. 


(312) 


where 


and 
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VP- 


=  Y    + 

c 


'"I'-'+ig"}  ; 


iJj-      =   magnitude  of  the  vehicle  velocity  vector 


I     =  t-tL 


(313) 


an  incremental  time  into  the  future. 


The  error,   «sc   ,  between  the  predicted  path  and  the  desired  path 
is  measured  at   £A/   different  points,  each  one  AS    times  i 
(  L     =  1,  ...,  BN       )  inches  along  the  predicted  path  ahead  of  the 
vehicle.  The  measured  errors  are  individually  weighted  and  the  results  are 
averaged  to  yield  Ej 
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3.9.1.3   Error  Det  erminat  ion 


The  manner  of  error  calculation  is  outlined  below. 
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V 


Figure  3,31     ERROR  CALCULATION 
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It  can  be  seen  in  Figure  3,31  that 


*  P;     ~    A  V  Pi  *      U, 

—  =  tan   *<■  (314) 


C."  u'p 


Since  <2.L       is  measured  normal  to  the  projected  path  at  point 
(  ^vpi     >     H  \JPi    ),  the  slope  along  which  &L       is  measured,  relative  to 
the  <j'        axis,  can  be  determined  from 

^-   -  tan   P.  (315) 

ctci 

In  the  present  case  the  indicated  derivative  is  equal  to  the  ratio  of  the  X. 
velocity  at  point  (  ^'\jpL    >     y'vPc     )  t0  tne  H         velocity. 
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Combining  Equations  (314),  (315)  and  (316),  one  obtains 
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From  the  desired  path  input 


(317) 


(318) 


where 


/I  =  1  for  yPi  -     'P, 

n  =  3  for  YP,  <   yPc  +  Y'pz 

n.  =  5  f  or  yp2  ^y'pc    ~    Y'P3 

*■  =  7  for  Yr*  *  yk  ~    Y'P+ 

^  =  9  for  Yp+  *  1j  'pL 
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Equation  (317)  can  now  be  written  as 

c,-^.     "  ~^W1W         (319) 

Solution  of  Equation  (319)  allows  determination  of  the  error  via 

1/2 


ei\  L1"^/   ^vp.J    MJ/>.    rVPyJ\  (320) 

Thev  algebraic  sign  of  Q,        is  determined  by  the  sign  of  *j  P     in'the  normal- 
tangential  reference  system  at  each  point  i 
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SGW  <e  .)   -    SGN  (gp  )  (322) 


The        u    th  error       &i.       would,   therefore,   be 
e  ■   -    \eA    SGN  (eL) 


(323) 


The  lateral  acceleration  required  to  displace  the  path  of  the  vehicle 
by  d    at  a  distance  d       ahead  of  the  vehicle  is 

2 

From  Equations  (311)  and  (324),  it  can  be  seen  that  the  steer  angle  change 
required  to  null  an  error  €.       at  a  distance  <J       ahead  of  the  vehicle  on  its 
projected  path  is 

AVf   r  71 e  (325) 

o 
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Because  the  driver  model  samples  the  error  between  the  desired  path  and  the 
projected  path  at  E.N      evenly  spaced  points  AS     apart,  the  steer  angle 
change  required  to  null  the  error  &i        of  the    i   th  point  a  distance 
from  the  vehicle  is 

A\=   (XTT? :e<-  <326> 

3.9.1.4   Importance  Weighting  Function 

If  the  importance  weighting  of  error   <£c  is   1/1/1^  ,  then  the 
resultant  steer  angle  change  is  the  weighted  average  of  the  steer  angle 
changes  .computed  for  individual  sample  points. 


or 


i=l 


EN 


AVW  =  KP  H  w£<  WI< e 


1-1 


where 


wLL       =        -z  is  the  error  weighting  function, 


WZ,        is  the  importance  weighting  function, 


l/  2(a  +  b)(l  +  Kti  ulj) 


•p  E       gt" is  the  control 

gain.   The  driver's  estimate  of  the  vehicle 
handling  parameter  (1  +  £j  \jfa  )  is 

probably  a  constant.  Therefore,  for  the 
preview-predictor  driver  model  Ke>       should 
be  determined  for  the  normal  velocity 
and  then  held  constant. 
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(328) 


3.9.1.5   Filter  Structure 


A  simplified  model  of  the  physiological  operator  is  assumed  which 
incorporates  a  time  delay   £T    a  possible  lead  Ti_  and  a  lag   7r 


This  filter  structure  corresponds  to  the  first  order  effects  of  the 
neurological  and  muscular  systems  of  a  human  driver  (Reference  17).  The 
parameters   7r  and  ~TL      determine  the  exact  shape  of  the  output  curve 

/^YFj'Ct)   ,  and   V      sets  the  time  delay  before  onset  of  the  smoothed  output 
(Figure  3. 32) . 

All  of  these  parameters  are  variable  inputs  of  the  model  and, 
therefore,  readily  permit  alterations  in  the  exact  shape  of  the  filter  output 
in  order  to  match  a  particular  set  of  driver  characteristics. 


IDEAL 


Figure  3.32        NEURO  -  MUSCULAR  FILTER  CHARACTERISTICS 
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Mathematically,  the   V  th  desired  steer  change  from  time  zh,' 
yields  an  actual  steer  change  at  time  -t       of 
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To  account  for  the  time  delay 
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(330) 


and  since  the  transients  are  assumed  to  have  disappeared  after  one  second, 
the  steady  state  response  becomes: 


A  fFU)  = A  ¥< ,  ^r  (t  -tj-r)z  1 


t 


(331) 


The  time  functional  form  of  the  actual  steer  angle  would  then  be  the  sum  of 
the  j    th  independent  inputs 


7 


(332) 


and  the  ideal  steer  angle  would  be 

t/emnr 


(333) 


where  BfADT     is  the  time  interval  between  driver  samples  (Figure  3.32), 
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If  I  frC*)  I  >     &  ,  where  .  3_f        is 

the  maximum  steer  angle,  then  (334) 

and  (335) 

3.9.2     Skid  Control  Mode 

To -effect  the  control  of  a  vehicle  skid  it  is  first  necessary  to 
detect  the  skid.  This  is  usually  done  by  the  "feel"  of  the  brake,  steering 
and/or  accelerator  systems  along  with  visual  and  inert ial  cues.  No  attempt 
has  yet  been  made  to  simulate  the  detection  mechanism  per  se,  except  for  a 
comparison  of  the  vehicle  slip  angle  with  threshold  values. 

Three  types  of  vehicle  skids  are  considered:  braking  skids,  accel- 
eration skids  and  cornering  skids.   In  all  cases,  the  first  attempt  to  control 
the  skid,  in  the  model,  is  through  a  nulling  of  all  wheel  torques.   For  a  pure 
rear  wheel  acceleration  skid,  torque  nulling  should  be  a  sufficient  action  to 
return  control  of  the  vehicle  to  the  driver,  providing  that  the  yaw  velocity 
has  not  attained  an  excessive  value.   If  a  skid  is  of  such  severity  that  the 
nulling  of  the  wheel  torques  will  not  reduce  the  vehicle  slip  angle  below  the 
threshold  value,  steer  angle  commands  are  generated  in  an  attempt  to  regain 
control.   In  such  a  case,  instead  of  following  the  desired  path  through 
steer  angle  inputs,  the  driver  model  attempts  to  orient  the  vehicle,  through 
steering  commands,  so  that  its  heading  corresponds  to  the  direction  of  its 
velocity  vector. 

From  Figure  3.33  it  can  be  seen  that 

Oc=  arcfain  (g)  (336) 
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x'i 


Figure  3.33     VEHICLE  SLIP  ANGLE 


y1 


where   Qc  is  the  vehicle  slip  angle.   The  skid  is  detected  when  the  vehicle 
slip  angle,   ©^  ,  exceeds  the  lower  skid  threshold  T/z-)        .  The  brake  and 
throttle  controls  are  set  to  zero  until  the  condition  represented  by  Equation 
(337)  is  no  longer  met,  following  which  the  model  returns  to  the  speed  con- 
trol mode. 


if  \9J>Tr.1 


then    APD  -   Ft 


3&k: 


o 


(337) 
(338) 


For  cases  of  extreme  skidding  where  the  vehicle  slip  angle  exceeds 
the  higher  skid  threshold  Trz    (^2   >  ^'^    >  Patn  following  is 
discontinued,  and  the  steer  angle  commands  are  made  proportional  to  the 
front  wheel  slip  angle. 


if     l0cj  >  Tz^ 

then    &%j    =  Zs  (-&+&) 


(339) 
(340) 
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If  the  vehicle  has  not  rotated  more  than  ninety  degrees  during  a  skid  correction 
maneuver,  the  driver  will  return  to  his  path-following  behavior.  However, 
if  the  vehicle  orientation  has  changed  more  than  ninety  degrees,  the  driver 
model  will  bring  the  vehicle  to  rest. 

3.9.3     Speed  Control  Mode 

The  speed  control  section  of  the  driver  model  is  much  less  complex 
than  the  path-following  mechanism.  The  model  determines  the  difference 
between  the  desired  and  the  actual  forward  speed  of  the  vehicle,  AV,  and 
then  attempts  to  null  the  difference  within  a  specified  distance.   It  will 
be  possible  in  the  future  to  have  the  driver  model  include  a  calculation  of 
the  desired  nulling  distance.  However,  for  the  sake  of  initial  simplicity, 
the  nulling  distance  is  prespecified  as  a  tabular  input  for  each  speed  change 
required.   For  speed  maintenance,  on  the  other  hand,  the  acceleration  is 
specified  as  a  function  of  AV.  The  speed  control  outputs  are  in  terms  of 
inches  of  accelerator  pedal  deflection  or  pounds  of  brake  pedal  effort,  The 
modeled  driver  "assumes"  that  a  vehicle  deceleration  of  0.1  g  will  result 
for  each  ten  pounds  of  applied  brake  pedal  force  and  also  that  there  is  a 
linear  relationship  between  change  in  throttle  position  and  change  in  the 
magnitude  of  the  vehicle  acceleration, 

The  order  of  calculation  for  the  speed  control  section  of  the  driver 
model  is  as  follows.  The  velocity  error  is  calculated  as 

AV  =  PS-  Ur  (341) 

A  set  of  threshold/indifference  levels  7i/   and  Tsz      are  applied 
for  positive  and  negative  AV,  respectively.  A  desired  acceleration  Pq^ 
is  calculated  using  the  desired  nulling  distance,  t?Z£T 

Pa*  =  "Sir  (342) 
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For  a  change  of  speed  task  the  new  desired  speed,  OSj^     ;  the  time 
of  the  change,  "^C£    ;  and  the  desired  distance  to  null  AV,  PX5TX^  ,   are 
among  the  data  read  into  the  program.  Dli>T        is  set  equal  to  PISTl  ^  at 
-i^_        .  At  each  time  increment  after   -tjc       ,    VIST        is  reduced  by  the 
distance  traveled  by  the  vehicle, 

CXsr -  prsr-  aT  bmdt  (344) 

until  DP3T  <^  0  when  it  is  reset  equal  to  DLSTL*    and  the  speed  maintenance 
mode  is  resumed. 

The  speed  maintenance  mode  consists  of  exactly  the  same  equations 
as  the  speed  change  mode  except  no  periodic  update  of  DTST     at  each  calcu- 
lation time  is  done,  i.e.,  the  desired  acceleration  is  proportional  to  the 

i 
vehicle  velocity  and  the  magnitude  of  the  error  times  a  gain,  £>rST" 

For   P«*   >  0,  an  accelerator  pedal  deflection  APP        is  the  model  output. 

APD  =  A?f>  +  >52  *W  (345) 

The  initial  value  of  APD       is  found  by  solving  the  net  horsepower  versus 
speed  function  at  £)5  /   for  constant  velocity.  The  vehicle,  therefore, 
starts  all  computer  runs  in  equilibrium  at  the  first  desired  speed.   Since 

APO        can  never  exceed  APOMA"/   ,  nor  be  less  than  zero,  for  cases  of 
Equation  (345)  where  this  happens,  APD       is  set  to  APD  MAX     or  zero  at 
the  exceeded  limit. 

For   iJk/   >  7^     ,  the  braking  indifference  level,  the  applied 
brake  pedal  force,  Fbr.£       >   ^s 

fa*  =  K*  Q*>  (346) 

where  K^*         is  the  assumed  response  of  the  brake  system. 
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3.9.4     Driver  Vehicle  Interface 

3.9.4.1  Braking 

The  output  from  the  driver  model  is  a  force  acting  on  the  brake 
pedal.   This  force  is  converted  to  a  master  cylinder  hydraulic  pressure 
which  is  subsequently  used  to  determine  brake  torques  as  described  in 
Section  3.6.2  by  the  following  equation 

P       -  3pFi  (F&ZIc)  ■hfyp2.(F0£k:)7~  (347) 

3.9.4.2  Engine  Torque 

Engine  output  is  controlled  by  the  accelerator  pedal  deflection 
(  /\PD     )  and  the  engine  speed.  The  accelerator  pedal  deflection  is  normalized 
with  respect  to  the  maximum  (  /\PDhAA)(   )  to  produce  a  throttle  setting  as 
defined  in  Section  3.6.2.2 

5  "  APDMA/  (348) 

The  throttle  setting  and  engine  speed  are  then  used  to  interpolate 
between  the  input  engine  torque  values.   The  engine  speed  is  determined  from 
the  average  drive  rotational  velocity  of  the  drive  wheels,  the  rear  axle 
ratio  and  the  transmission  ratio: 

RPMB  =  §  (AZ/)(TZ) ( RPSi  +  flPSu, )  (349) 

A  simplified  representation  of  an  automatic  transmission  is  used 
with  the  driver  model  to  determine  the  transmission  ratio  which  allows 
transmission  upshifts  and  downshifts  at  input  engine  speeds.   The  instantaneous 
transmission  ratio  is  determined  by  calculating  the  engine  speed  using 
Equation  (349)  with  each  input  gear  ratio  and  comparing  the  calculated  engine 
speed  with  the  input  engine  speed  shift  points.   Note  that  down  shifts  are 
assumed  to  occur  only  if  the  previous  value  of  engine  torque  is  negative. 
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Letting  XQ-     represent  the  gear  number  at  the  previous  time 
increment,  then  the  gear  ratio  is  GEAR  x&-    anc*  the  upshift  engine  speed  is 
\lGm 


[Z&  ,!<?+) 


If  R.PI\AE     >    ^&RXkZGr  j-&  ,  an  upshift  occurs,  i.e. 

T£  =  GEAR  X6rH 

Similarly,  for  downshifts  (  TGlB     <  0) 

If  &PME-    <     ^&&$T£r,x6-\       >   a  downshift  occurs,  i.e. 


3.10      Aerodynamic  and  Rolling  Resistances 
(HVOSM  Vehicle  Dynamics  Version) 


The  effects  of  aerodynamic  and  rolling  resistances  are  grossly 
approximated  by  means  of  a  force,  £Fx5         ,  acting  along  the  X  axis  of  the 
sprung  mass.  The  magnitude  of  the  force  is  determined  from  an  empirical 
relationship  which  relates  the  driving  force  that  is  required  to  overcome 
the  cited  resistances  to  the  component  of  linear  velocity  along  the  vehicle 
X  axis.   For  simplicity,  it  is  assumed  that  ^Fx<,        acts  through  the  sprung 
mass  center  of  gravity.  The  primary  objective  of  this  simplified  analytical 
treatment  is  to  approximate  the  effects  of  aerodynamic  resistance  on  stopping 
distances. 

A  number  of  empirical-analytical  relationships  can  be  found  in  the 
literature.  The  following,  representative  relationship  is  taken  from 
Reference  9: 

F«.  -  0.0011Z8A\/z-h  Vj(0.O0OI395M-h  0.016JG)  (35°) 
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where 

PfZ       =  force  required  to  over  come  aero- 
dynamic and  rolling  resistances,  lbs. 

A   =  front  area,  square  feet. 

V   =  vehicle  speed,  miles  per  hour. 

W      =  total  vehicle  weight,  lbs. 

Converting  Equation  (350)  to  a  general  format  that  is  suitable  for  the  present 
simulation  program,  the  following  expression  is  obtained: 


ZF«  = 


m 


'  0     ,    fer  \U\  *  1  iVt/sec 

.  -  Cj  ujvl 'CtU-C5^nV    ,  -for  Ivl  1 1  in/s&c 


(351) 


where 


Li        =  scalar  component  of  linear  velocity 
along  the  X  axis,  in/sec. 

CfjCz  C$  =  coefficients  to  approximate  the  total 
force  required  to  overcome  aero- 
dynamic resistance. 

In  applications  where  measured  resistance  data  are  available,  the 
coefficients  C. ,  C  ,  C  can  be  directly  selected  to  achieve  an  empirical 
fit.   In  the  absence  of  such  measured  data,  Equation  (350)  can  provide  a 
basis  for  approximating  C  ,  C_  and  C  ,  by  means  of  the  following 
conversions: 

C   =   (3.6415X10"6)A 
C2  =   (7.926X10"6)W 


C   =   0.01676W 
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where 


A  =  frontal  area,  square  feet,  and 
W  =  total  vehicle  weight,  lbs. 

For  A  ^       26 .4  square  feet  and  W   ~   3600  pounds,  the  following, 
representative  values  are  obtained  for  the  coefficients: 

C   =  9.6136X10"5 
C   =  0.02853 
C   =  60.336 

A  rolling  resistance  torque  is  approximated  by  a  linear  function 
of  tire  normal  load  and  is  applied  directly  to  the  appropriate  wheel  spin 
equation  of  motion: 

ZrmL    =   "  Cf&hA,  ?*[  S$n(ZPSL)  (352) 

3.11  Resultant  Forces  and  Moments 

The  resultant  forces  on  the  vehicle  that  act  on  the  unsprung  masses 
in  the        X        and        U        directions  are  determined  from  the  following 
relationships: 

ZFXU    =    ifixuL+FcvuL  +  fswi)  (353) 

Zfvu  =  2  (fa/*i  rfcju  +6s»d  (354) 


-1 
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Forces  acting  on  the  unsprung  masses  in  the   2   direction  must 
be  treated  individually  in  view  of  the  unsprung  mass  degrees  of  freedom: 

FlUt  =    FztUL+fciUL   +F611M  (355) 

The  resultant  force  that   is  transmitted  through  the  suspensions  in 
the     z       direction  is  obtained  from 

ZFt^St  +  St+Sz+Si  (356) 

The  resultant  moments  acting  on  the  sprung  and  unsprung  masses  arising 
from  tire  and  suspension  forces  are  given  below  for  the  three  suspension 
options. 

Independent  Front/Solid  Axle  Rear  Suspension 

-  (Fyus  *  Fyu4 )  (*  **)  +  5* (T^  fAr>"^ 

-SCTrfr+fify)    -h   (%-%)%Jl  (357) 

\  (358) 

f  F/U4  (h  +si  V-Tefa/Z  i-  h+Q*  YK4) 
2^Vw       ~      Ftju,(a  +  li,&$(yKl)    4  Fy-U2  (a^hzcoscxhz) 

+  FXuz  (Tftt+A  fe  -  k2  &*fiki  )-FXu,  (%  H +A  Tm  i-k,  a>*fa , ) 

(359) 
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2A/*t      =      Fiu3  (1k/2-/xfc  +k<**fa)  -  hu4  fattyfa  -h*  a>*fa) 

f  (%-%)%/%  (360) 

Independent  Front  and  Rear  Suspension 

ZNfa      =    $jTF/l+&THpz)  ~$,(TF/2i-A%Fl)  +$4yi+bT^)  -%(TElt+m&) 

-  Fijui(  %f  r*/+k,  cos  XKl )  -  ^«2  ^zf  *£  *  ^2  ^&2  ) 

+  Fxa+CzztSi  +  h+costirt)  (362) 

2  AJ<a<     =     F^ ,  (a  -fk,  cos  <xKi )  t  Fyul  (a.  +  kz  cosaru ) 

-  fyuz  (b  ~  ^i co^m)  -  Fyw  Cb  -  h^CD^a^ ) 

<■  Fml  (  TFM+ATHF2-kzco$fa2)  -  FXUI(TF/Z  tATHF/  +ht  cosfa, ) 

(363) 

■f  Fxw  CTnk  +&Th&  -  K4 &*  fa )  "  Fxm  (Tfi/2  tAThR3+  h*>  coffin*,) 


Solid  Front  and  Rear  Axles 


2A/#i      -    -(fyui  +-Fyu.1)[2Fi-<fl)  -CF^u^tFiiuA.)(ig.tSa>) 

+($r*,)%Ffr    +fy-%)Ts/2  (364) 
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+-f*uz  (zF+<*i  +pF  -TF<f>p/z  +k2  coeXnz)  +  EfUl  (%+&,+?  +Tzfa/ltkdcc6K^ 
+  hiM&e^P  -Tifa/Z  +IH &***+)  (365) 

*  ^xu2 ft/2  V/=  *  "^  ^Au )  " FXU3  (T* /2  7^  ^  ^fe) 

*  ^>m  ft  /2  +/  ^  -  k+  a>e^>K4)  ,-  366 ) 

-  fyui  (fci-Tfhfo  ^,cob\)  -  fm  (pp  -TF<j>F/2  +k2co$Xkl) 

+  (5r62)T^/2  (367) 

2Mte     =      l^uiC'bte'/fa  +  ha&ffa)   '^4^^^-^^J 

-  ^3  0>  * Ze &  /z  f/l/  ^^1/ )  -/^ Lp-1kfal£+  h+  cot  Xkf) 

+  (S3-S4)T$/2  (370) 

Note  that  driveline  torque  reactions  are  included  in  the  roll  moments  acting 
on  the  sprung  and  unsprung  masses  in  the  HVOSM-VD2  version  if  the  drive  axle 
is  solid.   For  the  case  of  a  solid  rear  axle,  with  rear  drive,  the  following 
terms  are  added  to  the  equations  for  sprung  mass  roll  moment  and  rear  axle 
roll  moment: 

ZN<j>u  =    ...  +n(m)DR 
£%  =    ...  -12 era),* 
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For  the  case  of  a  solid  front  axle  with  front  drive  the  following  terms 
are  added  to  the  equations  for  sprung  mass  roll  moment  and  front  axle 
roll  moment: 

Zttju    =  ...  +12  (TQ)df 
Zhlfc    =    ...  -12  (TQ)pf 
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